Advanced Search
Article Contents

Kinetic Energy Budgets during the Rapid Intensification of Typhoon Rammasun (2014)


doi: 10.1007/s00376-022-2060-z

  • In this study, Typhoon Rammasun (2014) was simulated using the Weather Research and Forecasting model to examine the kinetic energy during rapid intensification (RI). Budget analyses revealed that in the inner area of the typhoon, the conversion from symmetric divergent kinetic energy associated with the collocation of strong cyclonic circulation and inward flow led to an increase in the symmetric rotational kinetic energy in the lower troposphere. The increase in the symmetric rotational kinetic energy in the mid and upper troposphere resulted from the upward transport of symmetric rotational kinetic energy from the lower troposphere. In the outer area, both typhoon and Earth’s rotation played equally important roles in the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy in the lower troposphere. The decrease in the symmetric rotational kinetic energy in the upper troposphere was caused by the conversion to asymmetric rotational kinetic energy through the collocation of symmetric tangential rotational winds and the radial advection of asymmetric tangential rotational winds by radial environmental winds.
    摘要: 以往对于热带气旋快速增强机制的研究通常集中在最大风速所在的内核区域进行。本研究对2014年的台风个例威马逊进行了分析,除内核区域外,还分析了台风外部区域快速增强期间的动能收支。对低层旋转动能收支的分析显示,在台风内区,气旋式切向风和向内的辐散流相互作用造成的对称辐散动能向对称旋转动能的转换导致了对流层低层对称旋转动能的增加。对流层中高层对称旋转动能的增加是由于对称旋转动能从对流层低层向中高层输送引起的。在台风外区,台风和地球自转在对流层低层的对称辐散动能向对称旋转动能的转换中发挥了同样重要的作用。对流层高层对称旋转动能的减弱是通过对称切向旋转风和径向环境风对非对称切向旋转风的径向平流的相互作用引起的低层旋转动能向非对称旋转动能的转换所致。本研究的结果表明,低层入流在快速增强期间发挥了以下重要作用:低层入流与向外位势梯度方向相反,向外位势梯度将可将位能转换为对称辐散动能;低层入流与气旋切向流相互作用,将对称辐散动能转换为对称旋转动能;低空入流导致台风内区向上运动,将动能从对流层低层输送到对流层中高层。
  • 加载中
  • Figure 1.  Map of the three nested model domains. The color bar indicates the height (m) of the topography.

    Figure 2.  (a) Maximum wind speed (solid, m s−1) and minimum sea level pressure (dashed, hPa) and (b) track of Typhoon Rammasun (2014) from simulation (blue) and observation (orange). Time is in LST.

    Figure 3.  Vertical profiles of symmetric rotational kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

    Figure 4.  Vertical profiles of symmetric divergent kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

    Figure 5.  Vertical profiles of asymmetric rotational kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

    Figure 6.  Vertical profiles of asymmetric divergent kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

    Figure 7.  Radial variation of tangentially averaged partitioned components of (a) the conversion between environmental kinetic energy and symmetric rotational kinetic energy $\left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right)$ and the flux divergence of symmetric rotational kinetic energy $\left({{\rm{FK}}}_{{\psi }_{\rm s}}\right)$, (b) the conversion between asymmetric rotational kinetic energy and symmetric rotational kinetic energy $\left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right)$, and (c) $ ({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}) $ at 300 hPa, (d) ${{\rm{FK}}}_{\psi s}$,(e) $ \left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right) $, $ ({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}) $, the conversion between symmetric divergent kinetic energy and symmetric rotational kinetic energy $\left({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}\right)$ and ${{\rm{FK}}}_{{\psi }_{\rm s}}$, and (f) $ ({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}) $ at 500 hPa and (g–i) $ ({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}) $ at 850 hPa. Units: m2 s−3.

    Figure 8.  Radial variation of tangentially averaged partitioned components of (a) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and (b) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $. Except for $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ in (a) and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},\mathrm{C}} $ in (b), the other partitioned components of $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ are negligibly small. Units: m2 s−3.

    Figure 9.  Radial variation of tangentially averaged partitioned components of (a) $ {({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2} $ (m2 s−3) at 300 hPa, and (b) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and (c) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ at 850 hPa. Except for $ {({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2\theta } $ in (a), $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ in (b), and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},\mathrm{C}} $ in (c), the other partitioned components of $ {({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2} $, $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $, and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ are negligibly small. Units: m2 s−3.

    Figure 10.  Radial variation of tangentially averaged partitioned components of (a1, a2) $ {({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1} $ and (b1, b2) ${({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}}}$ at (a1, b1) 300 hPa, and (a2, b2) 500 hPa, and (a3) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and (b3) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ at 850 hPa. Units: m2 s−3.

    Figure 11.  Radial variation of tangentially averaged partitioned components of the conversion from potential energy to symmetric rotational kinetic energy $\left[{(P,K}_{{\chi }_{{\rm{s}}}}\right)]$ at 850 hPa. Units: m2 s−3.

    Figure 12.  Radial variation of tangentially averaged partitioned components of the flux divergence of asymmetric rotational kinetic energy $\left({{\rm{FK}}}_{{\psi }_{{\rm{a}}}}\right)$ at (a) 300 hPa and (b) 500 hPa. Units: m2 s−3.

    Figure 13.  Schematic diagram that summarizes dominant dynamic processes associated with the change in symmetric rotational kinetic energy during the rapid intensification of Typhoon Rammasun (2014). $ {K}_{{\psi }_{\rm s}} $ is the symmetric rotational kinetic energy; $ \left({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}\right) $ is the conversion between symmetric divergent kinetic energy ($ {K}_{{\chi }_{\rm s}} $) and symmetric rotational kinetic energy; $ \left({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}\right) $ is the conversion between asymmetric rotational kinetic energy ($ {K}_{{\psi }_{\rm a}} $) and symmetric rotational kinetic energy; $\mathrm{FK}_{{\psi }_{\rm s}}$ is the flux divergence of symmetric rotational kinetic energy; $ {v}_{{\psi }_{\rm s}\theta } $ is the tangential component of symmetric rotational wind; r is the radial distance from typhoon center; f is the Coriolis parameter due to the Earth’s rotation.

  • Bougeault, P., and P. Lacarrere, 1989: Parameterization of orography-induced turbulence in a mesobeta-scale model. Mon. Wea. Rev., 117(8), 1872−1890, https://doi.org/10.1175/1520-0493(1989)117<1872:POOITI>2.0.CO;2.
    Chang, C.-C., and C.-C. Wu, 2017: On the processes leading to the rapid intensification of Typhoon Megi (2010). J. Atmos. Sci., 74(4), 1169−1200, https://doi.org/10.1175/JAS-D-16-0075.1.
    Chen, G. H., 2016: Determination of the effect of initial inner-core structure on tropical cyclone intensification and track on a beta plane. Adv. Atmos. Sci., 33(8), 945−954, https://doi.org/10.1007/s00376-016-5241-9.
    Chen, X. M., J. A. Zhang, and F. D. Marks, 2019: A thermodynamic pathway leading to rapid intensification of tropical cyclones in shear. Geophys. Res. Lett., 46(15), 9241−9251, https://doi.org/10.1029/2019GL083667.
    Chen, X. M., Y. Q. Wang, K. Zhao, and D. Wu, 2017: A numerical study on rapid intensification of Typhoon Vicente (2012) in the South China Sea. Part I: Verification of simulation, storm-scale evolution, and environmental contribution. Mon. Wea. Rev., 145(3), 877−898, https://doi.org/10.1175/MWR-D-16-0147.1.
    Chen, X. M., Y. Q. Wang, J. Fang, and M. Xue, 2018: A numerical study on rapid intensification of Typhoon Vicente (2012) in the South China Sea. Part II: Roles of inner-core processes. J. Atmos. Sci., 75(1), 235−255, https://doi.org/10.1175/JAS-D-17-0129.1.
    Cheng, C.-J., and C.-C. Wu, 2020: The role of WISHE in the rapid intensification of tropical cyclones. J. Atmos. Sci., 77(9), 3139−3160, https://doi.org/10.1175/JAS-D-20-0006.1.
    Cheng, X. P., J. F. Fei, X. G. Huang, and J. Zheng, 2012: Effects of sea spray evaporation and dissipative heating on intensity and structure of tropical cyclone. Adv. Atmos. Sci., 29(4), 810−822, https://doi.org/10.1007/s00376-012-1082-3.
    DeMaria, M., C. R. Sampson, J. A. Knaff, and K. D. Musgrave, 2014: Is tropical cyclone intensity guidance improving. Bull. Amer. Meteor. Soc., 95(3), 387−398, https://doi.org/10.1175/BAMS-D-12-00240.1.
    DeMaria, M., J. Kaplan, and J.-J. Baik, 1993: Upper-level eddy angular momentum fluxes and tropical cyclone intensity change. J. Atmos. Sci., 50(8), 1133−1147, https://doi.org/10.1175/1520-0469(1993)050<1133:ULEAMF>2.0.CO;2.
    Ding, Y. H., and Y. Z. Liu, 1986: Study on Kinetic Energy Budget of a Typhoon-The conversion between Kinetic energies of divergent winds and non-divergent winds. Scientia Sinica (Series B), 29(4), 397−410.
    Dudhia, J., 1989: Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model. J. Atmos. Sci., 46(20), 3077−3107, https://doi.org/10.1175/1520-0469(1989)046<3077:NSOCOD>2.0.CO;2.
    Ek, M. B., K. E. Mitchell, Y. Lin, E. Rogers, P. Grunmann, V. Koren, G. Gayno, and J. D. Tarpley, 2003: Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model. J. Geophys. Res., 108(D22), 8851, https://doi.org/10.1029/2002JD003296.
    Elsberry, R. L., T. D. B. Lambert, and M. A. Boothe, 2007: Accuracy of Atlantic and eastern North Pacific tropical cyclone intensity forecast guidance. Wea. Forecasting, 22(4), 747−762, https://doi.org/10.1175/WAF1015.1.
    Feng, T., G.-H. Chen, R.-H. Huang, and X.-Y. Shen, 2014: Large-scale circulation patterns favourable to tropical cyclogenesis over the western North Pacific and associated barotropic energy conversions. International Journal of Climatology, 34(1), 216−227, https://doi.org/10.1002/joc.3680.
    Hong, S. Y., and J.-O. Lim, 2006: The WRF single-moment 6-class microphysics scheme (WSM6). Journal of the Korean Meteorological Society, 42, 129−151.
    Kain, J. S., 2004: The kain- fritsch convective parameterization: An update. J. Appl. Meteor., 43(1), 170−181, https://doi.org/10.1175/1520-0450(2004)043<0170:TKCPAU>2.0.CO;2.
    Kanada, S., and A. Wada, 2015: Numerical study on the extremely rapid intensification of an intense tropical cyclone: Typhoon Ida (1958). J. Atmos. Sci., 72(11), 4194−4217, https://doi.org/10.1175/JAS-D-14-0247.1.
    Kaplan, J., and M. Demaria, 2003: Large-scale characteristics of rapidly intensifying tropical cyclones in the North Atlantic basin. Wea. Forecasting, 18(6), 1093−1108, https://doi.org/10.1175/1520-0434(2003)018<1093:LCORIT>2.0.CO;2.
    Kornegay, F. C., and D. G. Vincent, 1976: Kinetic energy budget analysis during interaction of Tropical Storm Candy (1968) with an extratropical frontal system. Mon. Wea. Rev., 104(7), 849−859, https://doi.org/10.1175/1520-0493(1976)104<0849:KEBADI>2.0.CO;2.
    Kurihara, Y., M. A. Bender, and R. J. Ross, 1993: An initialization scheme of hurricane models by vortex specification. Mon. Wea. Rev., 121(7), 2030−2045, https://doi.org/10.1175/1520-0493(1993)121<2030:AISOHM>2.0.CO;2.
    Kurihara, Y., M. A. Bender, R. E. Tuleya, and R. J. Ross, 1995: Improvements in the GFDL hurricane prediction system. Mon. Wea. Rev., 123(9), 2791−2801, https://doi.org/10.1175/1520-0493(1995)123<2791:IITGHP>2.0.CO;2.
    Leroux, M.-D., M. Plu, D. Barbary, F. Roux, and P. Arbogast, 2013: Dynamical and physical processes leading to tropical cyclone intensification under upper-level trough forcing. J. Atmos. Sci., 70(8), 2547−2565, https://doi.org/10.1175/JAS-D-12-0293.1.
    Leroux, M.-D., M. Plu, and F. Roux, 2016: On the sensitivity of tropical cyclone intensification under upper-level trough forcing. Mon. Wea. Rev., 144(3), 1179−1202, https://doi.org/10.1175/MWR-D-15-0224.1.
    Li, M. X., F. Ping, J. Chen, and L. R. Xu, 2016: A simulation study on the rapid intensification of Typhoon Megi (2010) in vertical wind shear. Atmospheric Science Letters, 17(12), 630−638, https://doi.org/10.1002/asl.713.
    Li, X. F., 1993: On the dynamics of vortex propagation and associated asymmetric gyres. PhD dissertation, University of Hawaii at Manoa, 102 pp.
    Liang, X. S., and D. G. M. Anderson, 2007: Multiscale window transform. Multiscale Modeling & Simulation, 6(2), 437−467, https://doi.org/10.1137/06066895X.
    Miyamoto, Y., and T. Takemi, 2015: A triggering mechanism for rapid intensification of tropical cyclones. J. Atmos. Sci., 72(7), 2666−2681, https://doi.org/10.1175/JAS-D-14-0193.1.
    Mlawer, E. J., S. J. Taubman, P. D. Brown, M. J. Iacono, and S. A. Clough, 1997: Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave. J. Geophys. Res., 102(D14), 16 663−16 682,
    Montgomery, M. T., and R. J. Kallenbach, 1997: A theory for vortex rossby-waves and its application to spiral bands and intensity changes in hurricanes. Quart. J. Roy. Meteor. Soc., 123(538), 435−465, https://doi.org/10.1002/qj.49712353810.
    Moon, I.-J., I. Ginis, T. Hara, and B. Thomas, 2007: A physics-based parameterization of air–sea momentum flux at high wind speeds and its impact on hurricane intensity predictions. Mon. Wea. Rev., 135(8), 2869−2878, https://doi.org/10.1175/MWR3432.1.
    Persing, J., M. T. Montgomery, and R. E. Tuleya, 2002: Environmental interactions in the GFDL hurricane model for hurricane opal. Mon. Wea. Rev., 130(2), 298−317, https://doi.org/10.1175/1520-0493(2002)130<0298:EIITGH>2.0.CO;2.
    Prasanth, S., D. R. Chavas, F. D. Jr. Marks, S. Dubey, A. Shreevastava, and T. N. Krishnamurti, 2020: Characterizing the energetics of vortex-scale and sub-vortex-scale asymmetries during tropical cyclone rapid intensity changes. J. Atmos. Sci., 77(1), 315−336, https://doi.org/10.1175/JAS-D-19-0067.1.
    Rappaport, E. N., and Coauthors, 2009: Advances and challenges at the national hurricane center. Wea. Forecasting, 24(2), 395−419, https://doi.org/10.1175/2008WAF2222128.1.
    Shieh, O. H., M. Fiorino, M. E. Kucas, and B. Wang, 2013: Extreme rapid intensification of Typhoon Vicente (2012) in the South China Sea. Wea. Forecasting, 28(6), 1578−1587, https://doi.org/10.1175/WAF-D-13-00076.1.
    Sun, Y., Z. Zhong, and Y. Wang, 2012: Kinetic energy budget of Typhoon Yagi (2006) during its extratropical transition. Meteorol. Atmos. Phys., 118(1−2), 65−78, https://doi.org/10.1007/s00703-012-0200-1.
    Wada, A., 2015: Unusually rapid intensification of Typhoon Man-yi in 2013 under preexisting warm-water conditions near the Kuroshio front south of Japan. Journal of Oceanography, 71(5), 597−622, https://doi.org/10.1007/s10872-015-0273-9.
    Wang, Y. P., X. P. Cui, X. F. Li, W. L. Zhang, and Y. J. Huang, 2016: Kinetic energy budget during the genesis period of Tropical Cyclone Durian (2001) in the South China Sea. Mon. Wea. Rev., 144(8), 2831−2854, https://doi.org/10.1175/MWR-D-15-0042.1.
    Wang, Y. Q., 2002: Vortex rossby waves in a numerically simulated tropical cyclone. Part I: Overall structure, potential vorticity, and kinetic energy budgets. J. Atmos. Sci., 59(7), 1213−1238, https://doi.org/10.1175/1520-0469(2002)059<1213:VRWIAN>2.0.CO;2.
    Yamada, H., Q. Moteki, and M. Yoshizaki, 2010: The unusual track and rapid intensification of Cyclone Nargis in 2008 under a characteristic environmental flow over the Bay of Bengal. J. Meteor. Soc. Japan, 88(3), 437−453, https://doi.org/10.2151/jmsj.2010-311.
    Yang, M.-J., D.-L. Zhang, and H.-L. Huang, 2008: A modeling study of Typhoon Nari (2001) at landfall. Part I: Topographic effects. J. Atmos. Sci., 65(10), 3095−3115, https://doi.org/10.1175/2008JAS2453.1.
    Zeng, Z. H., Y. Q. Wang, and L. S. Chen, 2010: A statistical analysis of vertical shear effect on tropical cyclone intensity change in the North Atlantic. Geophys. Res. Lett., 37(2), L02802, https://doi.org/10.1029/2009GL041788.
    Zhang, R. W., J. L. Huangfu, and T. Hu, 2019: Dynamic mechanism for the evolution and rapid intensification of Typhoon Hato (2017). Atmospheric Science Letters, 20(8), e930, https://doi.org/10.1002/asl.930.
    Zhao, H. R., L. G. Wu, and X. S. Liang, 2015: Preliminary application of the MWT method to separate tropical cyclone circulation. Advances in Meteorological Science and Technology, 5(6), 31−36, https://doi.org/10.3969/j.issn.2095-1973.2015.06.005. (in Chinese with English abstract
  • [1] ZHOU Lingli, DU Huiliang, ZHAI Guoqing, WANG Donghai, 2013: Numerical Simulation of the Sudden Rainstorm Associated with the Remnants of Typhoon Meranti (2010), ADVANCES IN ATMOSPHERIC SCIENCES, 30, 1353-1372.  doi: 10.1007/s00376-012-2127-3
    [2] Xun LI, Noel E. DAVIDSON, Yihong DUAN, Kevin J. TORY, Zhian SUN, Qinbo CAI, 2020: Analysis of an Ensemble of High-Resolution WRF Simulations for the Rapid Intensification of Super Typhoon Rammasun (2014), ADVANCES IN ATMOSPHERIC SCIENCES, 37, 187-210.  doi: 10.1007/s00376-019-8274-z
    [3] Tong XIE, Liguang WU, Yecheng FENG, Jinghua YU, 2024: Alignment of Track Oscillations during Tropical Cyclone Rapid Intensification, ADVANCES IN ATMOSPHERIC SCIENCES, 41, 655-670.  doi: 10.1007/s00376-023-3073-y
    [4] Xinghai ZHANG, Yihong DUAN, Yuqing WANG, Na WEI, Hao HU, 2017: A High-resolution Simulation of Supertyphoon Rammasun (2014) —— Part I: Model Verification and Surface Energetics Analysis, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 757-770.  doi: 10.1007/s00376-017-6255-7
    [5] Shuai YANG, Xiba TANG, Shuixin ZHONG, Bin CHEN, Yushu ZHOU, Shouting GAO, Chengxin WANG, 2019: Convective Bursts Episode of the Rapidly Intensified Typhoon Mujigae (2015), ADVANCES IN ATMOSPHERIC SCIENCES, 36, 541-556.  doi: 10.1007/s00376-019-8142-x
    [6] Qian WANG, Dajun ZHAO, Yihong DUAN, Shoude GUAN, Lin DONG, Hongxiong XU, Hui WANG, 2023: Super Typhoon Hinnamnor (2022) with a Record-Breaking Lifespan over the Western North Pacific, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 1558-1566.  doi: 10.1007/s00376-023-2336-y
    [7] P. L. S. Rao, 2000: The Influence of Systematic Errors on the Asian Summer Monsoon Circulation, ADVANCES IN ATMOSPHERIC SCIENCES, 17, 576-586.  doi: 10.1007/s00376-000-0021-x
    [8] Kin Sik LIU, Johnny C. L. CHAN, 2022: Growing Threat of Rapidly-Intensifying Tropical Cyclones in East Asia, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 222-234.  doi: 10.1007/s00376-021-1126-7
    [9] Donglei SHI, Guanghua CHEN, Ke WANG, Xinxin BI, Kexin CHEN, 2020: Evaluation of Two Initialization Schemes for Simulating the Rapid Intensification of Typhoon Lekima (2019), ADVANCES IN ATMOSPHERIC SCIENCES, 37, 987-1006.  doi: 10.1007/s00376-020-2038-7
    [10] Xiaohong BAO, Xiuping YAO, 2022: Intensity Evolution of Zonal Shear Line over the Tibetan Plateau in Summer: A Perspective of Divergent and Rotational Kinetic Energies, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 1021-1033.  doi: 10.1007/s00376-021-1302-9
    [11] Shou Shaowen, Li Shenshen, 1991: Diagnosis of Kinetic Energy Balance of a Decaying Onland Typhoon, ADVANCES IN ATMOSPHERIC SCIENCES, 8, 479-488.  doi: 10.1007/BF02919270
    [12] Wang Zuoshu, Zhang Ruojun, Peng Zhengyi, 1989: The Kinetic Energy Budget and Circulation Characteristics of the Tropical Storm Irma during AMEX Phase II, ADVANCES IN ATMOSPHERIC SCIENCES, 6, 414-423.  doi: 10.1007/BF02659076
    [13] Jianjun LIU, Feimin ZHANG, Zhaoxia PU, 2017: Numerical Simulation of the Rapid Intensification of Hurricane Katrina (2005): Sensitivity to Boundary Layer Parameterization Schemes, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 482-496.  doi: 10.1007/s00376-016-6209-5
    [14] Fei Jianfang, Lu Hancheng, 1996: Study on Instability in Baroclinic Vortex Symmetric Disturbance under Effect of Nonuniform Environmental Parameters, ADVANCES IN ATMOSPHERIC SCIENCES, 13, 461-470.  doi: 10.1007/BF03342037
    [15] Bo SUN, 2018: Asymmetric Variations in the Tropical Ascending Branches of Hadley Circulations and the Associated Mechanisms and Effects, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 317-333.  doi: 10.1007/s00376-017-7089-z
    [16] Runkua Yang, J. Shukla, P.J. Sellers, 1994: The Influence of Changes in Vegetation Type on the Surface Energy Budget, ADVANCES IN ATMOSPHERIC SCIENCES, 11, 139-161.  doi: 10.1007/BF02666542
    [17] ZUO Qunjie, GAO Shouting, and LÜ Daren, 2014: Eddy Kinetic Energy Study of the Snowstorm over Southern China in January 2008, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 972-984.  doi: 10.1007/s00376-013-3122-z
    [18] Yaoming MA, Weiqiang MA, Huaguang DAI, Lei ZHANG, Fanglin SUN, Jinqiang ZHANG, Nan YAO, Jianan HE, Zhixuan BAI, Yuejian XUAN, Yunshuai ZHANG, Yuan YUAN, Chenyi YANG, Weijun SUN, Ping ZHAO, Minghu DING, Kongju ZHU, Jie HU, Bian Bazhuga, Bai Juepingcuo, Zhuo Ma, Ren Qingnima, Suo Langwangdui, Yang Zong, Haikun WEN, 2023: Earth Summit Mission 2022: Scientific Expedition and Research on Mt. Qomolangma Helps Reveal the Synergy between Westerly Winds and Monsoon and the Resulting Climatic and Environmental Effects, ADVANCES IN ATMOSPHERIC SCIENCES, 40, 187-193.  doi: 10.1007/s00376-022-2166-3
    [19] Jianpu BIAN, Juan FANG, Guanghua CHEN, Chengji LIU, 2018: Circulation Features Associated with the Record-breaking Typhoon Silence in August 2014, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 1321-1336.  doi: 10.1007/s00376-018-7294-4
    [20] Lei LIU, Guihua WANG, Ze ZHANG, Huizan WANG, 2022: Effects of Drag Coefficients on Surface Heat Flux during Typhoon Kalmaegi (2014), ADVANCES IN ATMOSPHERIC SCIENCES, 39, 1501-1518.  doi: 10.1007/s00376-022-1285-1

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 08 March 2022
Manuscript revised: 23 May 2022
Manuscript accepted: 30 May 2022
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Kinetic Energy Budgets during the Rapid Intensification of Typhoon Rammasun (2014)

    Corresponding author: Xiaofan LI, xiaofanli@zju.edu
  • Key Laboratory of Geoscience Big Data and Deep Resource of Zhejiang Province, School of Earth Sciences, Zhejiang University, Hangzhou, Zhejiang 310027, China

Abstract: In this study, Typhoon Rammasun (2014) was simulated using the Weather Research and Forecasting model to examine the kinetic energy during rapid intensification (RI). Budget analyses revealed that in the inner area of the typhoon, the conversion from symmetric divergent kinetic energy associated with the collocation of strong cyclonic circulation and inward flow led to an increase in the symmetric rotational kinetic energy in the lower troposphere. The increase in the symmetric rotational kinetic energy in the mid and upper troposphere resulted from the upward transport of symmetric rotational kinetic energy from the lower troposphere. In the outer area, both typhoon and Earth’s rotation played equally important roles in the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy in the lower troposphere. The decrease in the symmetric rotational kinetic energy in the upper troposphere was caused by the conversion to asymmetric rotational kinetic energy through the collocation of symmetric tangential rotational winds and the radial advection of asymmetric tangential rotational winds by radial environmental winds.

摘要: 以往对于热带气旋快速增强机制的研究通常集中在最大风速所在的内核区域进行。本研究对2014年的台风个例威马逊进行了分析,除内核区域外,还分析了台风外部区域快速增强期间的动能收支。对低层旋转动能收支的分析显示,在台风内区,气旋式切向风和向内的辐散流相互作用造成的对称辐散动能向对称旋转动能的转换导致了对流层低层对称旋转动能的增加。对流层中高层对称旋转动能的增加是由于对称旋转动能从对流层低层向中高层输送引起的。在台风外区,台风和地球自转在对流层低层的对称辐散动能向对称旋转动能的转换中发挥了同样重要的作用。对流层高层对称旋转动能的减弱是通过对称切向旋转风和径向环境风对非对称切向旋转风的径向平流的相互作用引起的低层旋转动能向非对称旋转动能的转换所致。本研究的结果表明,低层入流在快速增强期间发挥了以下重要作用:低层入流与向外位势梯度方向相反,向外位势梯度将可将位能转换为对称辐散动能;低层入流与气旋切向流相互作用,将对称辐散动能转换为对称旋转动能;低空入流导致台风内区向上运动,将动能从对流层低层输送到对流层中高层。

    • Typhoons may experience a rapid intensification (RI) period before landfall, which often leads to large economic losses and deaths in the coastal areas. For example, the maximum wind speed in the super Typhoon Rammasun (2014), which is analyzed in this study, increased from 40 m s−1 to 60 m s−1 within 24 h before landfall at Hainan, China on 18 July 2014, causing a huge economic loss of over 26 billion Chinese Yuan and approximately 30 deaths. To effectively reduce the damage caused by typhoons, accurate typhoon forecasting is required to facilitate informed governmental decisions before landfall. When compared to the significant improvement in the track forecasts of tropical cyclones (TC) in recent decades, improvement in TC intensity forecasts have shown relatively slower progress (Elsberry et al., 2007; Rappaport et al., 2009; DeMaria et al., 2014). The intensity forecasts mainly rely on numerical model guidance; however, improvement of numerical predictions requires an in-depth understanding of the physical processes associated with the intensity change. Primary physical processes and factors associated with TC intensity change include environmental effects (DeMaria et al., 1993; Zeng et al., 2010; Feng et al., 2014), inner core dynamics (Montgomery and Kallenbach, 1997; Wang, 2002; Miyamoto and Takemi, 2015; Chen, 2016; Chen et al., 2019), underlying surface forcing (Yang et al., 2008; Cheng et al., 2012; Cheng and Wu, 2020), and convection and latent heat release (Kanada and Wada, 2015; Li et al., 2016).

      Many studies have contributed to enhancing the understanding of dynamic processes associated with the RI of TCs. The RI of Hurricane Opal (1995) was caused by the mean vertical advection and mean vorticity flux terms (Persing et al., 2002). The RI of TC Dora (2007) was related to the superposition of potential vorticity structure of the trough with strong deformation (Leroux et al., 2013). The RI of Typhoon Man-yi (2013) was associated with the enhanced mesovortex under the condition of reduced static stability (Wada, 2015). The RI of Typhoon Megi (2010) resulted from strengthened rotational circulations in the mid and upper troposphere through the transport of vorticity (Chang and Wu, 2017). The RI of Typhoon Vicente (2012) was largely affected by an upper-tropospheric “inverted” trough (Shieh et al., 2013). The analysis of the asymmetric impacts at multiple scales during RI found that the changes in asymmetric circulations were associated with the baroclinic conversion from available potential to kinetic energies, and the conversion between asymmetric kinetic energies (Prasanth et al., 2020).

      Budgets of tangential wind, vorticity, and kinetic energy have been analyzed to examine dynamic processes associated with the RI in previous studies. The RI of Vicente was related to mean radial flux of absolute mean vertical vorticity and vertical advection of the azimuthal mean tangential wind by the azimuthal mean vertical motion in the azimuthal mean tangential wind budget (Chen et al., 2017), and to tilting of the horizontal vorticity and upward vorticity advection in the mean vorticity budget of circular area (Chen et al., 2018). The positive tendency of the TC kinetic energy associated with the RI of Typhoon Hato (2017) resulted mainly from the conversion of kinetic energy of large-scale circulations in the barotropic eddy kinetic energy equation (Zhang et al., 2019).

      Previous studies have focused on the analyses of RI mechanisms using analyses that have been carried out in the inner core of the tropical cyclone where maximum winds are located. In the current case study of Typhoon Rammasun (2014), kinetic energy budgets in outer area during RI are also examined in addition to the analyses in the inner region. The questions that are discussed in this study are: How do kinetic energies during RI change vertically and horizontally? What dynamic processes control the kinetic energy changes in inner and outer areas and in the lower, middle, and upper troposphere during RI?

      The paper is organized as follows. An overview of super Typhoon Rammasun (2014) is presented, the model simulation design is introduced, and a comparison between observations and simulations is presented in section 2. The vertical structures of kinetic energy budgets are analyzed in section 3, and the dominant dynamic terms in kinetic energy budgets are discussed in-depth through partitioned analysis in section 4. Lastly, a summary is given in section 5.

    2.   Overview of Typhoon Rammasun (2014), model design, and kinetic energy budgets
    • Typhoon Rammasun (2014) originated in the northwest Pacific Ocean at 1400 Local Standard Time (LST; LST = UTC + 8) on 12 July 2014 and moved westward, developing into a severe typhoon. The first landfall was made on the southern Luzon Island of the Philippines at 1700 LST on 15 July. Rammasun (2014) weakened to a typhoon following landing, continued to move northwestward after 1400 LST on 17 July, and subsequently began to strengthen rapidly. It developed into a super typhoon at 1000 LST on 18 July and landed at Wenchang, Hainan, China at 1400 LST. During landing, the central pressure dropped to 918 hPa and the maximum wind speed increased to 72 m s−1, making it the strongest typhoon among those that made landfall in China since 1973.

      Kaplan and Demaria (2003) defined RI when change in the maximum wind speed of the typhoon in 24 h is more than 15 m s−1. This definition has been adopted by previous studies (Yamada et al., 2010; Leroux et al., 2016). Based on this definition of RI, the 13-h period from 0200 LST to 1400 LST on 18 July was considered as the RI of Rammasun (2014) in this study.

      Version 4 of the Weather Research and Forecasting (WRF) Model was used to simulate Typhoon Rammasun (2014), and Global Forecast System analysis data (GFS; 0.5° × 0.5°) were used to construct the initial and lateral boundary conditions. The model used two-way nesting and three nested domains (Fig. 1) with horizontal grid spacings of 18, 6, and 2 km and domain sizes of 299 (lon.) × 235 (lat.), 634 (lon.) × 538 (lat.), and 1015 (lon.) × 904 (lat.) grid points, respectively. The model had 50 vertical layers, with the top at 10 hPa. The model physics packages included the WRF single-moment 6-class cloud microphysics scheme (WSM6; Hong and Lim, 2006), BouLac planetary boundary layer (Bougeault and Lacarrere, 1989), Kain–Fritsch cumulus parameterization scheme (Kain, 2004), which was only used in the outermost domain, longwave radiation parameterization scheme (Mlawer et al., 1997), shortwave radiation parameterization scheme (Dudhia, 1989), Monin-Obukhov surface layer model (Moon et al., 2007), and Noah Land surface model (Ek et al., 2003). The model was integrated from 1400 LST on 16 July to 1400 LST on 19 July, which covered the 13-h RI period of Rammasun (2014) in the South China Sea. An average of the 13-h simulation data was used to conduct the following analysis of the kinetic energy budget.

      Figure 1.  Map of the three nested model domains. The color bar indicates the height (m) of the topography.

      Comparing the simulation of Rammasun (2014) with the Joint Typhoon Warning Center best track data (Fig. 2), it can be seen that the simulated maximum wind speed of Rammasun (2014) is consistent with the observation, and the difference in peak intensity between the simulation and observation only was 0.74 m s−1 at 1400 LST on 18 July. The simulated minimum sea-level pressure generally followed observations and the simulated typhoon track was slightly north of the observed track. The largest difference between the simulated and observed centers was less than 90 km at 1400 LST on 17 July. Thus, the simulation reasonably reproduced the intensity and track of Rammasun (2014).

      Figure 2.  (a) Maximum wind speed (solid, m s−1) and minimum sea level pressure (dashed, hPa) and (b) track of Typhoon Rammasun (2014) from simulation (blue) and observation (orange). Time is in LST.

      The change in TC intensity is accompanied by complex energy conversion processes, including the conversion between the kinetic energies of the TC and the energy conversion between the large-scale environmental and TC components, as well as the conversion between the potential and kinetic energies. When a cylindrical coordinate system is established at the center of the typhoon and moves with the typhoon, the kinetic energy associated with typhoon circulation can be decomposed into four components: symmetric rotational, symmetric divergent, asymmetric rotational, and asymmetric divergent components. The associated four kinetic energy budgets can then be derived accordingly (Li, 1993; Wang et al., 2016). The derivation of the four kinetic energy equations is provided in the appendix. Note that the tendencies of kinetic energy at a given time were calculated taking the differences in kinetic energy between the simulation data a time step ahead and a time step behind to guarantee the accuracy of kinetic energy budgets.

      The conversion between the environmental kinetic energy and TC kinetic energy in Eq. (A3) can be used to study the impacts of environmental wind, particularly environmental wind shear, on the change in typhoon intensity. We separated environmental and typhoon circulations using Kurihara’s filtering method (Kurihara et al., 1993, 1995), multiscale window transform (Liang and Anderson, 2007), and a combination of these two methods, respectively. We found that the combination of these two methods resulted in the better separation than the individual methods as shown in a study by Zhao et al. (2015). Therefore, a combination of these two methods was used in this study.

    3.   Analysis of area-averaged kinetic energy budgets
    • Since winds in the inner area of the typhoon are considerably larger than those in the outer area, the kinetic energy budgets in the inner area may be different from those in the outer area during RI. Thus, kinetic energy budgets were calculated in a circular area with a radius of 200 km (inner area) and annular areas with radii of 200−400 km and 400−600 km (outer area), respectively. Because the maximum wind of a typhoon is mainly tangential rotational wind, we first analyzed the symmetric rotational kinetic energy budget (Fig. 3). Analysis of symmetric rotational kinetic energy budget averaged over a circular area with a radius of 200 km showed that the symmetric rotational kinetic energy increased in the troposphere; however, the physical processes that were responsible for the increase varied with increasing height (Fig. 3a). Below 800 hPa, the increase in the symmetric rotational kinetic energy resulted from the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy $\left[\left({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}\right) > 0\right]$. From 800 hPa to 600 hPa, the increase in symmetric rotational kinetic energy was equally contributed from the convergence of the flux of the symmetric rotational kinetic energy and geopotential $\left({\mathrm{F}\mathrm{K}}_{{\psi }_{\rm s}} > 0\right)$ and the conversion from asymmetric rotational kinetic energy to symmetric rotational kinetic energy $ \left[\left({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}\right) > 0\right] $. Above 600 hPa, the increase in symmetric rotational kinetic energy originated from the convergence of the flux of the symmetric rotational kinetic energy and geopotential $\left({\mathrm{F}\mathrm{K}}_{{\psi }_{\rm s}} > 0\right)$. Near 300 hPa, the conversion from environmental kinetic energy to symmetric kinetic energy $ \left[\left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right) > 0\right] $ and $ {\mathrm{F}\mathrm{K}}_{{\psi }_{\rm s}} $ contributed equally to the increase in the symmetric rotational kinetic energy.

      Figure 3.  Vertical profiles of symmetric rotational kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

      Similarly, the analysis of symmetric rotational kinetic energy budget averaged over the annular area with a radius of 200–400 km revealed an increase in the symmetric rotational kinetic energy in the troposphere (Fig. 3b). However, the increase only appeared below 350 hPa in the calculation of the symmetric rotational kinetic energy budget averaged in the annular area with a radius of 400–600 km (Fig. 3c). The magnitude of the increase in symmetric rotational kinetic energy and the associated conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy significantly decreased as the radius increased (Fig. 3). The conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy made a dominant contribution to the increase in the symmetric rotational kinetic energy below 800 hPa within a radius of less than 200 km (Fig. 3a) and below 550 hPa within a radius of 200-600 km (Figs. 3b and 3c). Unlike the analysis of the 200 km circular area average, the analysis of the annular area average with the radius of 200–400 km showed that $\left({K}_{{\rm{e}}},{K}_{{\psi }_{\rm s}}\right)$ contributed more to the increase in the symmetric rotational kinetic energy than $ {{\rm{FK}}}_{{\psi }_{\rm s}} $ near 300 hPa (Fig. 3b). Analysis of the annular area average of a radius of 400–600 km revealed that the conversion from asymmetric rotational kinetic energy to symmetric rotational kinetic energy $ [\left({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}\right) > 0] $ played a dominant role in the increase in symmetric rotational kinetic energy in the mid-troposphere (Fig. 3c).

      Since the increase in symmetric rotational kinetic energy in the mid and lower troposphere mainly originated from the conversion from symmetric divergent kinetic energy, we analyzed the symmetric divergent kinetic energy budget (Fig. 4). In the mid and lower troposphere, the loss of symmetric divergent kinetic energy was compensated for by the conversion from potential energy $ \left[{(P,K}_{{\chi }_{\rm s}}) > 0\right] $, which is consistent with the results of Ding and Liu (1986). The balance between $ -\left({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}\right) $ and $ {(P,K}_{{\chi }_{\rm s}}) $ occurred below 800 hPa, 500 hPa, and 300 hPa in the budget averaged over the circular areas with a radius of 200 km (Fig. 4a), 200−400 km (Fig. 4b), and 400−600 km (Fig. 4c). The conversion from symmetric divergent kinetic energy to potential energy appeared above 800 hPa in the budget averaged over a circular area with a radius of 200 km (Fig. 4a). The budget averaged over the annual area with a radius of 200–400 km showed that the symmetric divergent kinetic energy gained from the conversion of symmetric rotational kinetic energy above 500 hPa was converted to potential energy from 500 hPa to 200 hPa, whereas it led to an increase in symmetric divergent kinetic energy above 200 hPa (Fig. 4b). The budget averaged over the annual area with a radius of 400–600 km revealed that in the upper troposphere, symmetric rotational kinetic energy and potential energy were converted to symmetric divergent kinetic energy, which led to an increase in symmetric divergent kinetic energy (Fig. 4c).

      Figure 4.  Vertical profiles of symmetric divergent kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

      The asymmetric rotational kinetic energy decreased with time (Fig. 5). The decreasing rate was one order of magnitude lower than the increasing rate of the symmetric rotational kinetic energy. The budget averaged in the circular area with a radius of 200 km showed that the decrease in asymmetric rotational kinetic energy below 300 hPa mainly resulted from the divergence of the flux of asymmetric rotational kinetic energy $\left({{\rm{FK}}}_{{\psi }_{{\rm{a}}}} < 0\right)$ below 800 hPa and from 300–400 hPa and from the conversion to symmetric rotational kinetic energy $\left[-\left({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}}\right) < 0\right]$ from 400–800 hPa (Fig. 5a). The budgets averaged in the annual areas with radii of 200-400 km and 400-600 km showed that the asymmetric rotational kinetic energy was converted to environmental kinetic energy $\left[\left({K}_{{\rm{e}}},{K}_{{\psi }_{{\rm{s}}}}\right) < 0\right]$, which caused a decrease in the asymmetric rotational kinetic energy below 400 hPa (Figs. 5b and 5c).

      Figure 5.  Vertical profiles of asymmetric rotational kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

      The change in asymmetric divergent kinetic energy was small (Fig. 6). The budget averaged in the circular area with a radius of 200 km revealed that below 600 hPa, the conversion from asymmetric divergent kinetic energy to symmetric rotational kinetic energy $ \left[-\left({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}\right) > 0\right] $ was balanced by the divergence of the flux of asymmetric divergent kinetic energy $\left({{\rm{FK}}}_{{\psi }_{\rm a}} < 0\right)$ (Fig. 6a). At approximately 500 hPa, the conversion from symmetric rotational kinetic energy to asymmetric divergent kinetic energy was offset by the conversion from asymmetric divergent kinetic energy to asymmetric rotational kinetic energy $ \left[-\left({K}_{{\chi }_{\rm a}},{K}_{{\psi }_{\rm a}}\right) < 0\right] $. Above 300 hPa, the conversion from environmental kinetic energy to asymmetric divergent kinetic energy $ \left[\left({K}_{\rm e},{K}_{{\chi }_{\rm a}}\right) > 0\right] $ was balanced by the divergence of the flux of asymmetric divergent kinetic energy as well as the conversion from asymmetric divergent kinetic energy to potential energy $ \left[\left(P,{K}_{{\chi }_{\rm a}}\right) < 0\right] $. The terms of the asymmetric divergent kinetic energy budgets averaged in the annular areas with radii of 200–400 km and 400-600 km (Figs. 6b and 6c) are considerably smaller than those averaged in the circular area with a radius of 200 km, which is not discussed in this study.

      Figure 6.  Vertical profiles of asymmetric divergent kinetic energy budgets (m2 s−3) averaged in (a) the circular area with a radius of 200 km and the annular area with radii of (b) 200–400 km and (c) 400–600 km, respectively.

    4.   Partitioned analysis of major physical terms
    • The analysis of kinetic energy budgets in the previous section showed that the important terms depend on the radius and vertical levels. To further explore the dominant processes associated with these major terms, partitioned analyses were conducted in the lower (850 hPa), middle (500 hPa) and upper (300 hPa) troposphere. The results from the previous section revealed that in a circular area with a radius of 200 km, the increases in symmetric rotational kinetic energy were related to the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy at 850 hPa, the convergence of the flux of symmetric rotational kinetic energy at 500 hPa, the flux convergence and the conversion from environmental kinetic energy to symmetric rotational kinetic energy at 300 hPa. Thus, the terms are partitioned as follows: following (A3a), (A4a) (A4b) and (A5a) in the Appendix, $\left({K}_{{\rm{e}}},{K}_{{\psi }_{{\rm{s}}}}\right)$, $\left({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}}\right)$, $\left({K}_{{\chi }_{{\rm{s}}}},{K}_{{\psi }_{{\rm{s}}}}\right)$ and ${{\rm{FK}}}_{{\psi }_{{\rm{s}}}}$ can be written as:

      where the first (subscript “1”) and second (subscript “2”) terms on the right-hand side of Eqs. (1)–(3) are from the first and second terms on the right-hand side of Eqs. (A3a), (A4a), and (A4b), which are the conversions related to the horizontal and vertical advection processes of environmental winds, asymmetric rotational winds and symmetric divergent winds, respectively; the third (subscript “3”) terms on the right-hand side of (1)–(3) are from the third terms on the right-hand side of (A3a), (A4a), and (A4b), which are the Coriolis force of environmental winds, asymmetric rotational winds and symmetric divergent winds, respectively; the first (subscript “1”) and second (subscript “2”) terms on the right-hand side of (4) are from the first and second terms on the right-hand side of (A5a), which are the convergence of horizontal and vertical fluxes of symmetric rotational kinetic energy, respectively; and the third (subscript “3”) term in the right-hand side of (4) is from the third term on the right-hand side of (A5a), which is the convergence of horizontal geopotential flux by symmetric rotational winds.

      Within the radius of 200 km from the typhoon center, the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy was determined by the horizontal advection processes $ {\left[\right({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1}] $ at 850 hPa (Fig. 7g). The flux convergence was due to the convergence of vertical flux of symmetric rotational kinetic energy $\left(\text{F}{\mathrm{K}}_{{\psi }_{\rm s}2}\right)$ at 500 hPa (Fig. 7d) and 300 hPa (Fig. 7a). The strong upward motions resulted from prevailing low-level convergent inflow of air that transported the symmetric rotational kinetic energy associated with the maximum wind from the lower troposphere to the mid and upper troposphere, strengthening the symmetric rotational kinetic energy in the mid and upper troposphere. The conversion from environmental kinetic energy to symmetric rotational kinetic energy was related to vertical advection processes $ {\left[\right({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2}] $ (Fig. 7a). The conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy occurred within 100 km from the typhoon center, with a maximum at approximately 30 km. The flux convergence at 500 hPa and 300 hPa and the conversion from environmental kinetic energy to symmetric rotational kinetic energy at 300 hPa occurred within 20−40 km with their maxima at approximately 30 km.

      Figure 7.  Radial variation of tangentially averaged partitioned components of (a) the conversion between environmental kinetic energy and symmetric rotational kinetic energy $\left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right)$ and the flux divergence of symmetric rotational kinetic energy $\left({{\rm{FK}}}_{{\psi }_{\rm s}}\right)$, (b) the conversion between asymmetric rotational kinetic energy and symmetric rotational kinetic energy $\left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right)$, and (c) $ ({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}) $ at 300 hPa, (d) ${{\rm{FK}}}_{\psi s}$,(e) $ \left({K}_{\rm e},{K}_{{\psi }_{\rm s}}\right) $, $ ({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}) $, the conversion between symmetric divergent kinetic energy and symmetric rotational kinetic energy $\left({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}\right)$ and ${{\rm{FK}}}_{{\psi }_{\rm s}}$, and (f) $ ({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}) $ at 500 hPa and (g–i) $ ({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}) $ at 850 hPa. Units: m2 s−3.

      Analysis of the tangential average from 200 km to 400 km showed that the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy was related to the horizontal advection processes $ \left[{({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1}\right] $ from 200 km to 400 km and the Coriolis force $ \left[{({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{3}\right] $ from 300 km to 400 km (Fig. 7h). Figure 3b shows that the conversion from both asymmetric rotational $({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})$and symmetric divergent kinetic energy to symmetric rotational kinetic energy and flux convergence contributed equally to the increase in the symmetric rotational kinetic energy at 500 hPa. The comparison between 500 hPa (Fig. 7e), 300 hPa (Fig. 7b), and 850 hPa (Fig. 7h) revealed that the four terms that contributed to the increase in symmetric rotational kinetic energy at 500 hPa had negligibly small magnitudes, which will not be discussed in this study. The conversion from environmental kinetic energy to symmetric rotational kinetic energy at 300 hPa was determined by the vertical advection processes $ \left[{({K}_{{\rm{e}}},{K}_{{\psi }_{\rm s}})}_{2}\right] $ (Fig. 7b).

      Analysis of the tangential average from 400 to 600 km revealed that the conversion from symmetric divergent kinetic energy to symmetric rotational kinetic energy at 850 hPa was due to both $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{3} $ (Fig. 7i). The conversion from asymmetric rotational kinetic energy to symmetric rotational kinetic energy at 500 hPa was related to $ \left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1} > 0\right] $ (Fig. 7f). The conversion from symmetric rotational kinetic energy to asymmetric rotational kinetic energy at 300 hPa is related to $ \left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1} < 0\right] $ (Fig. 7c).

      Substituting Eq. (A1a) into the first term on the right-hand side of Eq. (A4b), ${({K}_{{\chi }_{{\rm{s}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1}$ can be separated as:

      Within a radius of 600 km, $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ at 850 hPa was controlled by $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ (Figs. 8a, 9b and 10a3). In cylindrical coordinates, $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ can be further partitioned into terms related to the horizontal advection $\left[{({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},{\rm{A}}}+ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},{\rm{B}}}\right]$, and the terms ${({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},{\rm{C}}}$ and $\left[{({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},{\rm{D}}}\right]$ that are generated due to the transformation from Cartesian coordinates to cylindrical coordinates, that is,

      Figure 8.  Radial variation of tangentially averaged partitioned components of (a) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and (b) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $. Except for $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ in (a) and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},\mathrm{C}} $ in (b), the other partitioned components of $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ are negligibly small. Units: m2 s−3.

      Within a radius of 600 km, $ {({K}_{{\chi }_{{\rm{s}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\psi }_{{\rm{s}}}} $ at 850 hPa is determined by the collocation of tangential component of symmetric rotational winds and radial component of symmetric divergent winds$\left[{({K}_{{\chi }_{{\rm{s}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\psi }_{{\rm{s}}},{\rm{C}}}\right]$ (Figs. 8b, 9c and 10b3). Maximum cyclonic winds and inward convergent flows occurred in the lower troposphere during RI. The interaction between strong symmetric cyclonic tangential winds and inward convergent radial flows led to the conversion of symmetric divergent kinetic energy to symmetric rotational kinetic energy in the lower troposphere, which further increased symmetric rotational kinetic energy and strengthened cyclonic circulations. This indicates that the low-level inflow played an important role during the RI.

      ${({K}_{{\rm{e}}},{K}_{{\psi }_{{\rm{s}}}})}_{2}$ at 300 hPa can be further partitioned into:

      Figure 9a shows that ${({K}_{{\rm{e}}},{K}_{{\psi }_{{\rm{s}}}})}_{2}$ is determined by ${({K}_{{\rm{e}}},{K}_{{\psi }_{{\rm{s}}}})}_{2\theta }$, suggesting that the conversion from environmental kinetic energy to symmetric rotational kinetic energy results from the collocation of vertical advection of tangential environmental winds and tangential component of symmetric rotational winds.

      Figure 9.  Radial variation of tangentially averaged partitioned components of (a) $ {({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2} $ (m2 s−3) at 300 hPa, and (b) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and (c) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ at 850 hPa. Except for $ {({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2\theta } $ in (a), $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ in (b), and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s},\mathrm{C}} $ in (c), the other partitioned components of $ {({K}_{\rm e},{K}_{{\psi }_{\rm s}})}_{2} $, $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $, and $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ are negligibly small. Units: m2 s−3.

      Substituting Eq. (A1a) into the first term on the right-hand side of Eq. (A4a), ${({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1}$ can be separated into:

      The tangential average of 400−600 km revealed that the increase in symmetric rotational kinetic energy owing to the conversion from asymmetric rotational kinetic energy at 500 hPa resulted from the horizontal advection by the environmental winds $\left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}}} > 0\right]$ (Fig. 10a2). The decrease in the symmetric rotational kinetic energy due to the conversion to asymmetric rotational kinetic energy at 300 hPa was related to the horizontal advection by the environmental winds $\left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}}} < 0\right]$ (Fig. 10a1). In cylindrical coordinates, $\left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}}}\right]$ can be further separated into terms related to the horizontal advection of asymmetric rotational circulation by environmental winds $\left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}},\mathrm{A}}+ ({K}_{{\psi }_{\rm a}}, {K}_{{\psi }_{\rm s}})_{1{\rm{e}},\mathrm{B}}+{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}},\mathrm{C}}+{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}},\mathrm{D}}\right]$, and the terms $\left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}},\mathrm{E}}\right]$, and $\left[{({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}},\mathrm{F}}\right]$ that are generated due to the transformation from Cartesian coordinates to cylindrical coordinates, that is,

      Figure 10.  Radial variation of tangentially averaged partitioned components of (a1, a2) $ {({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1} $ and (b1, b2) ${({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}})}_{1{\rm{e}}}$ at (a1, b1) 300 hPa, and (a2, b2) 500 hPa, and (a3) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1} $ and (b3) $ {({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}})}_{1{\psi }_{\rm s}} $ at 850 hPa. Units: m2 s−3.

      The conversion from asymmetric rotational kinetic energy to symmetric rotational kinetic energy due to the horizontal advection of asymmetric rotational circulations by environmental winds $\left[{({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\rm{e}}} > 0\right]$ at 500 hPa was caused by the collocation of the tangential component of symmetric rotational winds and the radial advection of asymmetric tangential rotational winds by the radial environmental winds $\left[{({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\rm{e}},\mathrm{C}} > 0\right]$ and the interaction between radial component of asymmetric rotational winds and tangential components of environmental winds and symmetric rotational winds $\left[{({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\rm{e}},\mathrm{E}} > 0\right]$ (Fig. 10b2). The conversion from symmetric rotational kinetic energy to asymmetric rotational kinetic energy due to the horizontal advection of environmental winds $\left[{({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\rm{e}}} < 0\right]$ at 300 hPa was mainly attributed to the collocation of tangential component of symmetric rotational winds and the radial advection of asymmetric tangential rotational winds by radial environmental winds $\left[{({K}_{{\psi }_{{\rm{a}}}},{K}_{{\psi }_{{\rm{s}}}})}_{1{\rm{e}},\mathrm{C}} < 0\right]$ (Fig. 10b1).

      ${(P,K}_{{\chi }_{{\rm{s}}}})$ can be partitioned into:

      Within a radius of 200 km, the conversion from potential energy to symmetric rotational kinetic energy at 850 hPa resulted from the work done through radial geopotential advection by symmetric convergent inflows $\left[{{(P,K}_{{\chi }_{{\rm{s}}}})}_{{\rm{r}}}\right]$ (Fig. 11). In the lower troposphere, the geopotential increased dramatically as the radius increased, while convergent inflows prevailed. The collocation of an outward geopotential gradient and prevailing inflows transferred the potential energy to symmetric rotational kinetic energy in the lower troposphere. Similar results were found in the tangential average within 200–400 km and 400–600 km (not shown).

      Figure 11.  Radial variation of tangentially averaged partitioned components of the conversion from potential energy to symmetric rotational kinetic energy $\left[{(P,K}_{{\chi }_{{\rm{s}}}}\right)]$ at 850 hPa. Units: m2 s−3.

      Following Eq. (A5b), ${\mathrm{F}\mathrm{K}}_{{\psi }_{{\rm{a}}}}$ can be separated into:

      where the first ($ {\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}1} $) and second ($ {\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}2} $) terms on the right-hand side of (11) are from those of (A5b), which are the convergence of horizontal and vertical fluxes of asymmetric rotation kinetic energy, respectively; the third (${\mathrm{F}\mathrm{K}}_{{\psi }_{{\rm{a}}}3}$) and fourth ($ {\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}4} $) terms on the right-hand side of (11) are from those of (A5b), which are the convergence of horizontal and vertical fluxes of dot product of symmetric and asymmetric rotational winds, respectively; the fifth term ($ {\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}5} $) on the right-hand side of (11) is from that of (A5b), which is the convergence of horizontal geopotential flux by asymmetric rotational winds.

      Within a radius of 400–600 km, the conversion from asymmetric rotational kinetic energy to symmetric rotational kinetic energy at 500 hPa was supported by both flux convergence of the dot product of symmetric and asymmetric rotational winds $ \left[{\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}3} > 0\right] $ and flux convergence of the geopotential field $ \left[{\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}5} > 0\right] $ (Fig. 12b). The conversion from symmetric rotational kinetic energy to asymmetric rotational kinetic energy at 300 hPa was balanced by the flux divergence of the dot product of the symmetric and asymmetric rotational winds $ \left[{\mathrm{F}\mathrm{K}}_{{\psi }_{\rm a}3} < 0\right] $ (Fig. 12a).

      Figure 12.  Radial variation of tangentially averaged partitioned components of the flux divergence of asymmetric rotational kinetic energy $\left({{\rm{FK}}}_{{\psi }_{{\rm{a}}}}\right)$ at (a) 300 hPa and (b) 500 hPa. Units: m2 s−3.

    5.   Summary
    • In this study, the RI of Typhoon Rammasun (2014) that occurred in 2014 was simulated using the WRF model. The simulation was compared with observations in terms of maximum winds, minimum sea level pressures and typhoon movement. The simulation data were used to analyze the kinetic energy budget during RI. The winds were first partitioned into environmental winds and typhoon circulations, and the typhoon circulations were further separated into symmetric and asymmetric rotational winds, and symmetric and asymmetric divergent winds in cylindrical coordinates. The symmetric rotational kinetic energy generally showed an increase within a radius of 600 km from the typhoon center during RI.

      The dominant dynamic processes associated with the increase in the symmetric rotational kinetic energy during the RI of Typhoon Rammasun (2014) are summarized in Fig. 13. In the lower troposphere, the increase in the symmetric rotational kinetic energy results from the conversion of symmetric divergent kinetic energy. The conversion is caused by the collocation of strong symmetric cyclonic tangential flows and symmetric inward radial flows. The loss of symmetric divergent kinetic energy was compensated for by the transfer of potential energy through the collocation of inward decreased geopotential and prevailing inward flows. In the inner area of the typhoon (e.g., within a radius of 200 km from the typhoon center), where the maximum winds were located, the typhoon rotation was considerably larger than Earth’s rotation. As a result, strong symmetric cyclonic tangential rotational winds with symmetric inward convergent radial flows under typhoon rotation were responsible for the increase in symmetric rotational kinetic energy. In the outer area (e.g., within the radii of 400–600 km), the typhoon rotation became weak, having the same order of magnitude as the Earth’s rotation. Thus, symmetric cyclonic tangential rotational winds with symmetric inward radial flows due to typhoon and the Earth’s rotation accounted for the increase in the symmetric rotational kinetic energy.

      Figure 13.  Schematic diagram that summarizes dominant dynamic processes associated with the change in symmetric rotational kinetic energy during the rapid intensification of Typhoon Rammasun (2014). $ {K}_{{\psi }_{\rm s}} $ is the symmetric rotational kinetic energy; $ \left({K}_{{\chi }_{\rm s}},{K}_{{\psi }_{\rm s}}\right) $ is the conversion between symmetric divergent kinetic energy ($ {K}_{{\chi }_{\rm s}} $) and symmetric rotational kinetic energy; $ \left({K}_{{\psi }_{\rm a}},{K}_{{\psi }_{\rm s}}\right) $ is the conversion between asymmetric rotational kinetic energy ($ {K}_{{\psi }_{\rm a}} $) and symmetric rotational kinetic energy; $\mathrm{FK}_{{\psi }_{\rm s}}$ is the flux divergence of symmetric rotational kinetic energy; $ {v}_{{\psi }_{\rm s}\theta } $ is the tangential component of symmetric rotational wind; r is the radial distance from typhoon center; f is the Coriolis parameter due to the Earth’s rotation.

      In the upper troposphere, the upward advection of the symmetric rotational kinetic energy in the mid and upper troposphere led to an increase in the symmetric rotational kinetic energy in the inner area. The radial advection of asymmetric tangential rotational winds by radial environmental winds generates a reduction in the symmetric rotational kinetic energy in the outer area.

      The results from this study showed that the low-level inflow played the following important roles during the RI of Rammasun (2014): (1) the low-level inflow did work against the outward geopotential gradient, which converted available potential energy to symmetric divergent kinetic energy; (2) the low-level inflow interacted with cyclonic flow, which converted symmetric divergent kinetic energy to symmetric rotational kinetic energy; and (3) the low-level inflow resulted in upward motion in the inner area, which transported kinetic energy from the lower troposphere to the mid and upper troposphere.

      Several differences between this study and previous studies are noted. First, the analysis of symmetric divergent kinetic energy budget in this study revealed that the increase in symmetric rotational kinetic energy associated with the RI eventually resulted from the conversion from available potential energy, but the tangential momentum budget did not have any term related to the geopotential gradient. Also, the analysis of the vorticity budget did not show that the solenoidal term associated with geopotential gradient was important in previous studies. Secondly, the winds were separated into environmental and typhoon components in this study, whereas they have not been separated as such in previous studies. Finally, this study presents a case analysis of kinetic energy budgets in inner and outer areas during RI, whereas previous studies focused on the analysis of RI mechanisms only in inner area.

      It should be understood that, while the kinetic energy budget analysis of rapidly intensifying TC has been conducted in this study, the analyses of slowly intensifying TCs are required to identify the dominate dynamic processes for RI in future studies.

      Acknowledgements. The authors thank the Training Center of Atmospheric Sciences of Zhejiang University and two anonymous reviewers for their constructive comments and suggestions. We also thank Huiyan XU for deriving kinetic energy equations, Liguang WU and Huarui ZHAO for providing the MWT filter codes, Xinyong SHEN for improving the simulation of Typhoon Rammasun (2014), and Liangliang LI, and Chi ZHANG for valuable discussions. This study was supported by the National Natural Science Foundation of China (Grant No. 41930967).

    APPENDIX
    • The governing momentum equation in a p-coordinate system can be written as:

      Here, V [= (u, v)] is the horizontal wind vector, f is the Coriolis parameter due to the Earth’s rotation, and $\phi$ is the geopotential, $ \boldsymbol{F} $ is frictional dissipation, $ \nabla $ is horizontal differential operator, and $ \boldsymbol{k} $ is a unit vector in vertical direction.

      The winds $ \left(\boldsymbol{V}\right) $ can be partitioned into environmental $ \left({\boldsymbol{V}}_{\rm e}\right) $ and tropical cyclone (TC) $ \left({\boldsymbol{V}}_{{\rm{ty}}}\right) $ components. TC winds can be separated into rotational $ \left({\boldsymbol{V}}_{\psi }\right) $ and divergent $ \left({\boldsymbol{V}}_{\chi }\right) $ components, which can be further separated into symmetric $ ({\boldsymbol{V}}_{{\psi }_{\rm s}},{\boldsymbol{V}}_{{\chi }_{\rm s}}) $ and asymmetric $ ({\boldsymbol{V}}_{{\psi }_{\rm a}},{\boldsymbol{V}}_{{\chi }_{\rm a}}) $ components in cylindrical coordinates. Thus, the wind vectors can be expressed as:

      Here $ {\psi }_{\rm s} $ and $ {\chi }_{\rm s} $ are symmetric componential of stream function and velocity potential, and $ {\psi }_{\rm a} $ and $ {\chi }_{\rm a} $ are asymmetric componential of stream function and velocity potential. ${\boldsymbol{V}}_{\rm{e}}$, $ {\boldsymbol{V}}_{{\psi }_{\rm s}} $, $ {\boldsymbol{V}}_{{\chi }_{\rm s}} $, $ {\boldsymbol{V}}_{{\psi }_{\rm a}} $, and $ {\boldsymbol{V}}_{{\chi }_{\rm a}} $ are environmental wind, symmetric components of rotational and divergent wind, and asymmetric components of rotational and divergent wind, respectively.

      By multiplying momentum equation (A1) with the symmetric (${\boldsymbol{V}}_{{\psi }_{\rm s}}$) and asymmetric ($ {\boldsymbol{V}}_{{\psi }_{\rm a}} $) rotational wind vectors and symmetric (${\boldsymbol{V}}_{{\chi }_{\rm s}}$) and asymmetric ($ {\boldsymbol{V}}_{{\chi }_{\rm a}} $) divergent wind vectors, respectively, the corresponding kinetic energy equations can be written as:

      where

      Here, $ \mathrm{D}\mathrm{K} $ is the kinetic energy tendency of every wind component, K and P are kinetic and potential energies, respectively, r and $ \theta $ are radial distance and tangential angle, and $ {{r}}_{0} $ is radius of circular area for average. Equations (2a–d) are budgets of symmetric and asymmetric rotational kinetic energies and symmetric and asymmetric divergent kinetic energies, respectively. Equations (3a–d) are the conversion terms from environmental kinetic energy to TC kinetic energy, Eqs. (4a–f) are the conversion terms between TC kinetic energy components, Eqs. (5a–d) are flux divergence terms, Eqs. (6a–d) are the residual terms mainly due to dissipation (Kornegay and Vincent, 1976; Sun et al., 2012; Wang et al., 2016), Eqs. (7a–d) are the correlation terms of the wind components with their tendencies. Equations (8a–b) are the conversion terms between the potential energy and divergent kinetic energy. The angle bracket represents the area average that is defined as Eq. (9).

      Wang et al. (2016) derived similar kinetic energy equations for cylindrical coordinates. The differences between the kinetic energy equations in this study and those derived by Wang et al. include the following: (1) In this study, the winds are separated into environmental and typhoon components; thus, the environmental effects on typhoon intensity can be examined; however, Wang et al. (2016) did not perform this separation. (2) Our study indicates the presence of the divergence of the flux of dot product of rotational wind and symmetric divergent wind in symmetric divergent kinetic energy budget and the divergence of the flux of dot product of rotational wind and asymmetric divergent wind in asymmetric divergent kinetic energy budget [Eqs. (A5c)–(A5d)]; however, they are absent in the study by Wang et al. [Eqs. (1c)–(1d)]. (3) The term R in the study of Wang et al. [Eqs. (1a)–(1d)] can be simplified to the term R in this study [Eqs. (A7a)–(A7d)].

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return