Characteristics of Atmospheric Kinetic Energy Spectra during the Intensiﬁcation of Typhoon Lekima (2019)

: The intensiﬁcation of Typhoon Lekima (2019) is simulated with the Weather Research and Forecasting model to study the atmospheric horizontal kinetic energy (HKE) spectra and corresponding spectral HKE budgets under the control of real tropical cyclone (TC). The results show that the TC has the ability to modify the canonical atmospheric energy spectrum during its evolution, which is dominated by its rotational mode. With the intensiﬁcation of Lekima, the HKE spectrum in the troposphere swells over the central mesoscale and develops an arc-like shape. The stronger the TC, the more pronounced the arc-like shape is and the smaller scale it extends to. The roles various physical processes play at di ﬀ erent heights and horizontal scales during the intensiﬁcation of Lekima are investigated and the dependence of the e ﬀ ect of physical processes on scale and height is revealed. Meanwhile, the potential relationship between the intensiﬁcation of TC, the activation of energy activity at smaller scales, and the downscale extension of the arc-like spectral shape is found.


Introduction
For the past few decades, despite large improvement in the track forecast of tropical cyclone (TC), there is almost no substantial improvement in the intensity forecast for all lead times [1]. Model skill in predicting TC formation and rapid intensification is still very limited [2]. One reason why TC intensity is so poorly predicted lies in that, compared to track, which is more dependent on the environmental conditions, TC intensity can be more strongly dependent on internal dynamics and moist convection, which occur on smaller scales, and are more chaotic, less well understood, and intrinsically less predictable [3]. In this context, a better understanding of TC internal dynamics associated with its intensity variation is essential for improving intensity forecasts [4,5].
Among various methods, energy spectrum analysis is known to reveal signals at different scales. It is thus particularly suitable for investigating physical processes associated with TC, which is known as a typical multiscale system. Atmospheric energy spectrum refers to the distribution (i.e., variance) of atmospheric energy with wavenumber and/or frequency, which can intuitively reflect the scale characteristics of energy. Since the pioneering work of Nastrom and Gage [6], based on data from aircraft flights, many researchers have examined the horizontal wavenumber (k h ) spectrum of the atmospheric energy, especially for the mesoscale range. Numerous observational studies have reached a consensus-the atmospheric horizontal kinetic energy (HKE) spectrum as a function of k h possesses robust canonical shape, which follows a k −3 h power law at the synoptic scale and a k −5/3 h power law at

Typhoon Lekima
Typhoon Lekima formed in the Northwestern Pacific at around 0900 UTC 4 August 2019, moved northwestward and subsequently experienced an intensification at sea. After developing into a super Typhoon at around 1500 UTC 7 August, Lekima entered the mature stage, which persisted to around 1800 UTC 9 August. It then made landfall in Zhejiang, China, moved northward and began to weaken continuously. According to the best-track data of Japan Meteorological Agency (JMA), Lekima reached its maximum intensity at around 1200 UTC 8 August, with a lowest minimum sea level pressure of~925 hPa and a highest 10-m maximum wind speed of~54 m s −1 .

Experimental Design
In this study, we perform a 48-h numerical simulation of Lekima from 0000 UTC 6 August to 0000 UTC 8 August, which involved the period of primary intensification and initial maturity of the Typhoon. The model used here is the Advanced Research version of WRF (ARW-WRF) model, version 3.6. The model domain is duplicately nested through two-way nesting, and the outer domain is centered at (20. For both domains, the simulation uses the Thompson microphysics scheme [25], the Yonsei University (YSU) boundary layer scheme [26], and the Rapid Radiative Transfer Model for general circulation models (RRTMG) longwave and shortwave radiation scheme [27]. Rayleigh damping is applied to the vertical velocity in the upper 5 km of the model domains to prevent artificial reflection of gravity waves from the model top [28]. In addition, the Kain-Fritsch cumulus parameterization scheme [29] is used in the outer domain. To provide the initial fields and lateral boundaries for the simulation, the European Center for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis data at a horizontal resolution of 0.25 • × 0.25 • and a time interval of 6 h are used. All analyses in this study are conducted based on the results from the inner domain, which are output at 1-h interval.

Overview of the Simulation
The simulated track, minimum sea level pressure, and maximum 10-m wind speed of Lekima are compared with the best-track data of JMA at 6-h resolution in Figure 1. It can be found that the simulated track of Lekima agrees well with the observation throughout the simulation (Figure 1a). The model reproduces well the intensification of Lekima during the simulated period, and achieves a lowest minimum sea level pressure of 938.8 hPa and highest 10-m maximum wind speed of 48.5 m s −1 (Figure 1b).    Figure 2a shows the vertical profile of the basic-state Brunt-Väisälä frequency during the simulation, which is defined as the inner domain-averaged value and reflects the static stability of the background atmosphere. The frequency is almost constant below z = 5 km, with a mean value of 1.26 × 10 −2 s −1 . It then decreases slowly from 1.28 × 10 −2 s −1 at z = 5 km to 6.48 × 10 −3 s −1 at z = 14 km, while experiences a sharp growth from z = 14 km and achieves a maximum value of 2.91 × 10 −2 s −1 at z = 18 km. Thus, the heights from 14 to 17 km can be roughly treated as a transition zone from troposphere to stratosphere.   The diabatic heating released by the microphysical processes can be output by the WRF model directly. The time-height distribution of the heating rate averaged over a 200 × 200 m 2 square horizontal area being centered in the Typhoon is shown in Figure 2b. During the simulation, the powerful heating always penetrates the whole troposphere, implying the vigorous deep convection throughout the intensification of Lekima. In addition, the heating does not distribute uniformly in the troposphere, but exhibits two relatively separated vertical subranges, with heights of 2-5 and 6-12 km, respectively. The lower subrange stems mainly from the release of heating by the condensation process; the upper subrange stems mainly from the release of heating by the condensation, freezing, and sublimation processes; and the gap of heating between these two subranges comes from the absorption of heating induced by the melting process (not shown).

Methodology
In WZ2018, they calculated the energy spectra using the two-dimensional (2D) discrete Fourier transform (DFT). However, using the DFT on regional domains with an aperiodic structure will result in the aliasing of large-scale variance into shorter scales; thus, destroying all usefulness of spectra at large wavenumbers; while the method of discrete cosine transform (DCT) can solve this problem effectively [30]. Thus, in this study, HKE spectra and various spectral budget terms are calculated based on the 2D DCT of fields over the inner domain. Letφ(k) be the DCT of field ϕ(x, y), where k = (k x , k y ) is the horizontal wave vector. The HKE spectrum per unit volume is defined as where ρ is the height-variant basic-state density and u = (u, v) is the horizontal velocity vector.
The HKE spectrum can be further decomposed into horizontally rotational and divergent components [E R (k) and E D (k)], respectively [18], which are defined as E R (k) = ρζ(k)ζ(k)/ 2|k| 2 and E D (k) = ρδ(k)δ(k)/ 2|k| 2 (2) where ζ = ∂v/∂x − ∂u/∂y is the vertical vorticity and δ = ∂u/∂x + ∂v/∂y is the horizontal divergence. Following Peng et al. [31], the spectral HKE budget formulation is given by where ∂E h /∂t is net tendency term of HKE; t h is spectral transfer term due to nonlinear interactions; F h↑ and F p↑ are HKE vertical flux and pressure vertical flux, respectively, and ∂ z F h↑ and ∂ z F p↑ are the vertical divergence of the corresponding vertical flux; C A→h represents the spectral conversion of available potential energy (APE) to HKE; Div h is the spectral tendency due to the three-dimensional (3D) divergence; H h is the spectral tendency due to diabatic processes; J h is the adiabatic nonconservation term; and D h is the dissipative term. The detailed expressions for the terms in Equation (3) are written as Div h (k) = ρû · DCT[u(∂ z w + ∇ · u)]/2 (8) H h (k) = c p ρĤ mπ (9) where DCT() is equivalent to ∧ (), ∇ is the horizontal gradient operator, w is the vertical velocity, c p is the specific heat of dry air at constant pressure, θ is the height-variant basic-state potential temperature; π = π − π = (p/p 0 ) R d /c p − (p/p 0 ) R d /c p is the perturbed Exner pressure, with p, p 0 , p and R d denoting pressure, reference surface pressure, basic-state pressure and gas constant for dry air, respectively; H m = (1 + 1.61q v )S θ + 1.61θS qv is the combined diabatic contributions, with q v denoting water vapor mixing ratio, θ denoting potential temperature, and S θ and S qv denoting the diabatic contributions to dθ/dt and dq v /dt, respectively; and D u represents the dissipation of u.
For the convenience of displaying the results, one-dimensional (1D) wavenumber spectra as a function of total horizontal wavenumber k h (= |k| = k 2 x + k 2 y ) are constructed by angular averaging over the wavenumber band k h − ∆k/2 ≤ |k| < k h + ∆k/2 as where ε represents any one of the spectral terms in this study; ∆k = π/(∆N), ∆ is the horizontal grid spacing, N = min (N i , N j ), and N i and N j are the number of grid points along the zonal and meridional directions, respectively.

Evolutions of the HKE Spectra
Time evolutions of the vertically averaged HKE spectra over the troposphere (z = 0-12 km) and lower stratosphere (z = 17-19 km) are shown in respectively. In the troposphere, it can be found from the distribution of HKE spectrum as a function of horizontal wavelength λ h (wavenumber k h ), which, with the decrease (increase) of λ h (k h ), the magnitude of HKE in general reduces sharply at the initial time (0000 UTC 6 August; Figure 3a), especially for λ h < 200 km (k h > π × 10 −5 rad m −1 ). It suggests that although the ECMWF ERA5 reanalysis data reach a horizontal resolution of 0.25 • , the atmospheric information they contain at mesoscale is still limited. The model then takes about 6 h to complete the spinup adjustment, which reflects in the variation of the HKE spectrum. The spectrum experiences a significant lift for λ h < 400 km (k h > π/2 × 10 −5 rad m −1 ) during the first 3 h and a continued lift for λ h < 80 km (k h > π/4 × 10 −4 rad m −1 ) during the second 3 h. Obviously, compared with larger scales, the spinup adjustment exerts more influence on λ h < 100 km Accompanied with the intensification of Lekima, the HKE in the troposphere then increases at various horizontal scales in the following time (Figure 3b-d). Note that the HKE spectrum has already reproduced well the −3 slope at the synoptic scale and the −5/3 slope at the mesoscale of the real atmosphere just after the completion of spinup (0600 UTC 6 August; Figure 3a). Thus, it would have kept this canonical spectral shape if energy at various scales could increase comparatively in the log coordinate, whereas it is not the case here. The HKE spectrum can be roughly divided into three wavelength regions: , the spectrum in general appears as lifting consistently; thus, the spectral shape over this two wavelength regions remains stable. Figure 4a-d for a zoomed-in view); however, the significant growth of energy firstly occurs at larger scales, and then extends to smaller scales. Consequently, the spectrum appears to be swelling during this wavelength region and shows an arc-like spectral shape.  The arc-like shape of HKE spectrum in the troposphere over the central mesoscale is an interesting phenomenon, which has also been revealed by the idealized TC simulation of WZ2018 and the observational hurricane study of Vonich and Hakim [32]. Beyond that, this phenomenon has not been identified by any other studies related to atmospheric energy spectrum. Therefore, it can be treated as a unique feature for TC. For the case in this study, it can be found from the evolution of the spectrum that before the intensification of Lekima, the HKE already possesses canonical spectrum, which is flat over the central mesoscale with a slope of approximately −5/3; while, with the rapid intensification of Lekima, the HKE spectrum gradually swells. In other words, the formation and development of this unique spectral feature is associated with the intensification of TC, which can be demonstrated by the observations. In Vonich and Hakim [32], they calculated the HKE spectra based on numbers of in situ observations and then classified the spectra according to storm intensity. Although they did not describe the hurricane spectra with expression similar to "arc-like shape" or "swelling", their conclusion is in essence consistent with ours. From their results (e.g., their Figure 7), the spectrum for tropical storms has a slope close to −5/3 over the central mesoscale. As the storm intensifies from category-1 to category-5, the swelling of the spectrum becomes increasingly obvious, which extends to increasingly smaller wavelengths. This leads them to the conclusion that "as storm intensity increases, spectral slope becomes steeper and the transition zone where spectral slope begins to steepen shifts downscale". The variation of their hurricane spectrum with storm intensity is in perfect agreement with the evolution of our TC spectrum during the intensification of Lekima. Note that in our simulation, Lekima reaches a maximum intensity of category-2 and the swelling of the spectrum extends to wavelengths (wavenumbers) around 100 km (2π × 10 −5 rad m −1 ); while in WZ2018, their idealized TC reaches a maximum intensity of category-4 and the swelling extends to wavelengths (wavenumbers) around 50 km (2π/5 × 10 −4 rad m −1 ). This also proves the relationship between TC intensification and the arc-like spectral shape.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 18 Note that in our simulation, Lekima reaches a maximum intensity of category-2 and the swelling of the spectrum extends to wavelengths (wavenumbers) around 100 km (2π × 10 −5 rad m −1 ); while in WZ2018, their idealized TC reaches a maximum intensity of category-4 and the swelling extends to wavelengths (wavenumbers) around 50 km (2π/5 × 10 −4 rad m −1 ). This also proves the relationship between TC intensification and the arc-like spectral shape.  Figure 3, but zoomed in for 100 < λh < 300 km (2π/3 × 10 −5 < kh < 2π × 10 −5 rad m −1 ).  Compared to that in the troposphere, the growth of HKE spectrum in the lower stratosphere shows different characteristics (Figure 5a-d). After the spinup adjustment of the model, the spectral shape keeps relatively stable throughout the intensification of Lekima. Differing from the canonical atmospheric spectrum, the lower-stratospheric HKE spectrum here has consistent −5/3 slope for λ h < 1000 km (k h > 2π × 10 −6 rad m −1 ), but not −3 slope at the large-scale end. Meanwhile, the spectrum shows flat shape throughout the mesoscale but not arc-like shape over the central mesoscale as in the troposphere. These features are also identical to those found in WZ2018.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 18 To further understand the formation of the unique HKE spectra for TC, Figure 6 shows the rotational kinetic energy (RKE) and divergent kinetic energy (DKE) spectra at selected times. The RKE has spatiotemporal variable spectral shape. At 0600 UTC 6 August, the RKE spectrum in the troposphere is steeper at large scales and shallower at smaller scales, and the visible swelling of the spectrum appears for 200 < λh < 400 km (π/2 × 10 −5 < kh < π × 10 −5 rad m −1 ; Figure 6a); while the RKE spectrum in the lower stratosphere is flat for λh < 1000 km (kh > 2π × 10 −6 rad m −1 ) with a slope approaching −5/3 (Figure 6b). After the intensification of Lekima at 0000 UTC 8 August, the swelling of the RKE in the troposphere is more significant, which extends downscale to λh = 100 km (kh = 2π × 10 −5 rad m −1 ; Figure 6c); while the lower-stratospheric RKE spectrum suffers much less change ( Figure  6d). The DKE spectrum, in contrast, always shows a flat shape, which changes little with time and height. 12   To further understand the formation of the unique HKE spectra for TC, Figure 6 shows the rotational kinetic energy (RKE) and divergent kinetic energy (DKE) spectra at selected times. The RKE has spatiotemporal variable spectral shape. At 0600 UTC 6 August, the RKE spectrum in the troposphere is steeper at large scales and shallower at smaller scales, and the visible swelling of the spectrum appears for 200 < λ h < 400 km (π/2 × 10 −5 < k h < π × 10 −5 rad m −1 ; Figure 6a); while the RKE spectrum in the lower stratosphere is flat for λ h < 1000 km (k h > 2π × 10 −6 rad m −1 ) with a slope approaching −5/3 (Figure 6b). After the intensification of Lekima at 0000 UTC 8 August, the swelling of the RKE in the troposphere is more significant, which extends downscale to λ h = 100 km (k h = 2π × 10 −5 rad m −1 ; Figure 6c); while the lower-stratospheric RKE spectrum suffers much less change (Figure 6d). The DKE spectrum, in contrast, always shows a flat shape, which changes little with time and height. From Figure 6, it is clear that in the troposphere, the HKE is dominated by the RKE over almost all wavelengths. Considering that the TC is characterized by intensive rotation, this is to be expected. Therefore, the tropospheric spectral shape of HKE is determined by the RKE spectrum to a large extent, and the gradual swelling of the HKE spectrum during the intensification of Lekima is induced by the variation of the RKE spectrum. This means that the unique arc-like spectrum for TC in the troposphere stems from the rotational mode of TC motion. When it comes to the lower stratosphere, the RKE spectrum becomes as flat as the DKE spectrum and the HKE tends to be dominated by the divergent mode, probably because that the vortex mode of the TC which is remarkable in the troposphere does not extend upward to higher heights, and the remnant of the TC outflow layer we examine here is strongly divergent. As both the RKE and the DKE have flat spectral shape with a slope of −5/3, the resultant lower-stratospheric HKE spectrum also inherits this feature.

Spectral HKE Budget Diagnosis
To determine the roles of different physical processes play during the intensification of Lekima, here we investigate the spectral budget of HKE using Equation (3). The 12-h period from 0000 UTC 7 August to 1200 UTC 7 is chosen for analysis. The earlier period is not chosen to avoid the impact of   From Figure 6, it is clear that in the troposphere, the HKE is dominated by the RKE over almost all wavelengths. Considering that the TC is characterized by intensive rotation, this is to be expected. Therefore, the tropospheric spectral shape of HKE is determined by the RKE spectrum to a large extent, and the gradual swelling of the HKE spectrum during the intensification of Lekima is induced by the variation of the RKE spectrum. This means that the unique arc-like spectrum for TC in the troposphere stems from the rotational mode of TC motion. When it comes to the lower stratosphere, the RKE spectrum becomes as flat as the DKE spectrum and the HKE tends to be dominated by the divergent mode, probably because that the vortex mode of the TC which is remarkable in the troposphere does not extend upward to higher heights, and the remnant of the TC outflow layer we examine here is strongly divergent. As both the RKE and the DKE have flat spectral shape with a slope of −5/3, the resultant lower-stratospheric HKE spectrum also inherits this feature.

Spectral HKE Budget Diagnosis
To determine the roles of different physical processes play during the intensification of Lekima, here we investigate the spectral budget of HKE using Equation (3). The 12-h period from 0000 UTC 7 August to 1200 UTC 7 is chosen for analysis. The earlier period is not chosen to avoid the impact of initial spinup adjustment of the model on the results. The later period is not chosen mainly for two reasons: (1) this period involves the mature stage of Lekima, and (2) Lekima in this period is too close to some islands, which leads to the local absence of physical quantities and affects the implementation of 2D DCT over the lower layers. Considering the particularity of the boundary layer and the vertical distribution of diabatic heating, the details of dynamics and energetics of HKE should vary with height in the troposphere. Thus, the troposphere is further divided into the lower troposphere (z = 0-1.5 km), midtroposphere (z = 2.5-5 km), and upper troposphere (z = 6-12 km). Figure 7 shows the various spectral budget terms averaged in time over the intensification of Lekima and vertically over different layers. For the energy conversion term C A→h (black solid lines in the left panels of Figure 7), a positive value means a conversion induced by the buoyancy effect from APE to HKE and vice versa. In the lower troposphere, the magnitude of C A→h is relatively weak at various horizontal scales ( Figure 7a); while with the increase of height, the positive buoyancy dominates the troposphere and the importance of C A→h becomes prominent. In the midtroposphere, except for sporadic scales such as wavelengths (wavenumbers) around 100 and 400 km (2π × 10 −5 and π/2 × 10 −5 rad m −1 ), it mainly makes a positive contribution, implying a dominant conversion from APE to HKE (Figure 7c). The sporadic negative contribution comes from the local-scale negative buoyancy forcing induced by the increase of air density, which originates from the evaporation and melting processes. In the upper troposphere, it becomes entirely positive, and develops two peaks at wavelengths (wavenumbers) of 300 and 1000 km (2π/3 × 10 −5 and 2π × 10 −6 rad m −1 ), respectively (Figure 7e). In general, a level of neutral buoyancy which is characterized as a transition zone from positive to negative buoyancy always exists near the tropopause. Thus, differing from the troposphere, the lower stratosphere is dominated by negative buoyancy forcing, and the C A→h term always makes a negative contribution at all scales, implying a conversion from HKE to APE (Figure 7g).
The pressure vertical flux F p↑ corresponds to the flux of IGW energy at lowest order; thus, the divergence term ∂ z F p↑ can reflect the effect of vertical transportation of IGWs on the budget of HKE. The wavelength-height distribution of F p↑ is shown in Figure 8a, and note that a positive (negative) F p↑ means a downward (upward) flux. On the whole, F p↑ is upward in the lower stratosphere and downward in the troposphere at various scales, implying the domination of IGWs which occur in the upper troposphere and then propagate to higher and lower layers.
There are also nonnegligible differences at different scales. For 400 < λ h < 1000 km (2π × 10 −6 < k h < π/2 × 10 −5 rad m −1 ), the boundary between upward and downward F p↑ is around z = 15 km; namely, the highest source region of the waves can reach the tropopause. With the reduction of height from 15 to 6 km, the magnitude of F p↑ increases, which suggests that the downward-propagating waves are not only generated at the tropopause, but also generated throughout the upper troposphere. The F p↑ achieves its local maximum around z = 6 km and decreases from 6 to 4 km, indicating that the IGW energy are deposited into this height range during the downward propagation of the waves. It then increases again and achieves another local maximum around z = 2.5 km. Obviously, a secondary wave source exists in the midtroposphere for z = 2.5-4 km. At lower heights, the F p↑ gradually decreases, and the IGW energy are powerfully deposited into the lower troposphere.
The vertical distribution of F p↑ for λ h < 100 km (k h > 2π × 10 −5 rad m −1 ) is similar to that for 400 < λ h < 1000 km (2π × 10 −6 < k h < π/2 × 10 −5 rad m −1 ), except that the boundary height between positive and negative F p↑ reduces to z = 13 km; namely, the highest source region of the waves drops in height. For wavelengths (wavenumbers) slightly longer (smaller) than 100 km (2π × 10 −5 rad m −1 ); however, the distribution is not quite the same. Firstly, weak upward F p↑ appears for heights below 4 km, which can be attributed to upward-propagating IGWs excited in the boundary layer. Secondly, the boundary between positive and negative F p↑ in the upper troposphere drops rapidly to z = 6 km. In other words, the upper troposphere is not an effective source region for downward-propagating waves at these scales.        , and (g,h) lower stratosphere, and in time over the intensification period. The legend for the left panels is given in (a) and that for the right panels in (b). The asterisk * represents any one of the budget terms. The spectra are multiplied by k h to preserve area in log-linear coordinates. The vertical distribution of p F ↑ for λh < 100 km (kh > 2π × 10 −5 rad m −1 ) is similar to that for 400 < λh < 1000 km (2π × 10 −6 < kh < π/2 × 10 −5 rad m −1 ), except that the boundary height between positive and negative p F ↑ reduces to z = 13 km; namely, the highest source region of the waves drops in height. For wavelengths (wavenumbers) slightly longer (smaller) than 100 km (2π × 10 −5 rad m −1 );  Although there is weak upward F p↑ at wavelengths (wavenumbers) slightly greater (smaller) than 100 km (2π × 10 −5 rad m −1 ), the lower troposphere is primarily the net inflow layer of downward F p↑ . Meanwhile, the upper troposphere is always the net outflow layer of F p↑ and the lower stratosphere is the net inflow layer of upward F p↑ . Consequently, the ∂ z F p↑ term (black shaded lines in the left panels of  (Figure 7g), respectively. The midtroposphere, however, is both the inflow and the outflow layer of F p↑ , depending on the wavelength; thus the contribution of ∂ z F p↑ with wavelength over this layer is complicated (Figure 7c).
The spectral transfer term t h (red solid lines in the left panels of Figure 7) represents the transfer of HKE among different horizontal scales, which is caused by nonlinear interactions. It actually reflects the role of the so called "energy cascade" on budget of HKE at different scales. By constructing nonlinear spectral flux of HKE using t h , one can identify the direction of energy cascade, which is a key topic of research on atmospheric energy (WZ2018). In this paper, however, as the limited-area domain we analyze has open lateral boundaries, the spectral flux integrated over the whole domain is not equal to zero, and the direction of energy cascade cannot be identified accurately; thus, we only discuss the distribution of t h . In the lower troposphere, the energy cascade makes a negative contribution at wavelengths (wavenumbers) around 300 km (2π/3 × 10 −5 rad m −1 ) and makes a positive contribution for the rest scales (Figure 7a). In the midtroposphere, the effect of energy cascade is complicated and t h crosses zero several times with the variation of wavelength. However, it makes a positive contribution at most scales ( Figure 7c). When it comes to the upper troposphere, the role of t h is weak, with negative contribution for 400 < λ h < 800 km (π/4 × 10 −5 < k h < π/2 × 10 −5 rad m −1 ) and 200 < λ h < 300 km (2π/3 × 10 −5 < k h < 2π × 10 −5 rad m −1 ), and positive contribution for the remainder (Figure 7e). The distribution of t h is relatively clear in the lower stratosphere, which contributes positively for λ h > 500 km (k h < 2π/5 × 10 −5 rad m −1 ) and negatively at smaller scales ( Figure 7g).
The F h↑ term is the HKE vertical flux, and the corresponding divergence term ∂ z F h↑ reflects the direct contribution of vertical transportation of convection to HKE. It can be found from Figure 8b that, for most scales, the F h↑ is consistently upward from the lower troposphere to the lower stratosphere, which suggests that convection within the troposphere transports the energy upward. Two prominent subranges of F h↑ are collocated with those of the diabatic heating (Figure 2b), implying the close connection between the convection intensity and release of heating. There is also sporadic weak downward F h↑ in the lower troposphere and midtroposphere for wavelengths (wavenumbers) around 100 km (2π × 10 −5 rad m −1 ) and below 60 km (above π/3 × 10 −4 rad m −1 ), which is likely to be associated with the evaporation and the drag effect of precipitation particles. Reflected in the budget of HKE, the ∂ z F h↑ term (red shaded lines in the left panels of Figure 7) in general makes a negative contribution in the lower troposphere ( Figure 7a) and a positive contribution over the higher layers (Figure 7c,e,g).
The 3D divergence term Div h and adiabatic nonconservation term J h (blue solid and dashed lines in the right panels of Figure 7, respectively) stem from the compressibility and adiabatic nonconservation of the atmosphere, respectively. They are both significant in the midtroposphere and upper troposphere, and always play the balanced positive and negative contributions (Figure 7d,f).
The diabatic term H h (green solid lines in the right panels of Figure 7) represents the direct effect of diabatic processes on HKE, which is only relatively remarkable in the midtroposphere with negative contribution (Figure 7d).
The dissipative term D h (green dashed lines in the right panels of Figure 7) represents the effect of dissipative processes. Because of boundary-layer friction, its role is significant in the lower troposphere, with large negative contribution for λ h > 300 km (k h < 2π/3 × 10 −5 rad m −1 ; Figure 7b). However, for higher layers, it plays weak role (Figure 7d,f,h).
Overall, despite differences in magnitude and specific scales the effect are exerted on, the contributions of various budget terms at different heights during the intensification of Lekima are qualitatively consistent with those of the idealized TC in WZ2018, which suggests that the roles of various physical processes play are relatively steady. The spectral transfer term is an exception, which shows more complicated distribution here than that in WZ2018, implying the complexity of energy cascade associated with TC in real atmosphere.
During the intensification of Lekima, according to the C A→h term, the system is in general baroclinically unstable in the troposphere and stable in the lower stratosphere. Due to the indetermination of the direction of energy cascade, the barotropic instability associated with the t h term cannot be determined. However, by comparing the magnitude of C A→h and t h , the barotropic process dominates the baroclinic process in the lower troposphere, and the opposite happens in the upper troposphere; while these two processes are comparable in the midtroposphere. In other words, with the increase of height in the troposphere, the relative importance of barotropic process weakens, while the importance of baroclinic process enhances. Based on a forced shallow-water model and primitive equation model, respectively, Hendricks et al. [33] and Williams [34] discussed the potential vorticity dynamics associated with TC, and found the importance of barotropic instability during the evolution of TC inner region. Although the focuses of their research and the models and methods they used were all different from us in this paper, their conclusions on barotropic and baroclinic processes were essentially consistent with the results here, which supports our finding from a more theoretical perspective.
In addition, note that the various terms in the troposphere are stagnant for λ h < 100 km (k h > 2π × 10 −5 rad m −1 ) compared with those at larger scales. On the one hand, this matches the fact that the swelling of HKE spectrum only extends downscale to wavelengths (wavenumbers) around 100 km (2π × 10 −5 rad m −1 ). On the other hand, it implies that the activation of energy activity at smaller scales is the key to the downscale extension of the arc-like spectral shape with the intensification of TC.

Discussion and Conclusions
In this paper, the intensification of Typhoon Lekima is simulated for the purpose of investigating the evolution and characteristics of HKE spectra and diagnosing the corresponding spectral HKE budgets for real TC.
Accompanied with the intensification of Lekima, the HKE spectra at various heights vary remarkably. After the completion of the intensification, the spectra reproduce the canonical atmospheric spectrum to some extent, but also exhibit significant differences from the canonical spectrum. In the troposphere, the spectrum is in perfect agreement with the canonical spectrum before the intensification of Lekima; while accompanied with the intensification, the spectrum is gradually modified, and an arc-like spectral shape forms over the central mesoscale and extends to smaller scales. The investigation on RKE and DKE shows that the appearance of the arc-like shape is a direct result of the increasing vorticity during the intensification of Lekima. Combined with the idealized TC simulation of WZ2018 and the observational hurricane study of Vonich and Hakim [32], this arc-like shape is a unique and robust feature for atmospheric HKE spectrum associated with TC. The stronger the TC, the more pronounced this feature is and the smaller scale it extends to. In the lower stratosphere, the spectrum is in perfect agreement with the canonical spectrum at the mesoscale with a slope of −5/3, but the −3 slope at the synoptic scale does not appear.
The canonical atmospheric spectrum is obtained by averaging numerous of observations, and is considered representative and universal to the earth's atmosphere. However, the results in this study demonstrate that the canonical spectrum does not always apply to all cases. In fact, the evolution of HKE spectrum here in the troposphere directly reflects the modification of the spectrum by the development of Lekima. It is thus a reasonable inference that weather systems with unique structure and vibrant energy activity have the potential of modifying the canonical atmospheric spectrum during their development and evolution. Future investigations on the atmospheric energy spectrum should focus upon the spectral feature associated with specific weather systems, and seek to determine the formation mechanism behind the potential difference from the canonical spectrum. For example, why does the enhancement of rotational mode of TC lead to the swelling of HKE spectrum over the central mesoscale?
The predictability of TC indicated by the energy spectra, which has not been discussed here, also deserves to be explored. According to Lorenz [35], atmospheric system governed by a k −p power law has limited predictability when p < 3 and unlimited predictability when p ≥ 3. The arc-like shape of HKE spectrum over the central mesoscale means that p for TC system changes dramatically over a relatively narrow scale range, which may imply that the predictability of TC is complicated at the mesoscale.
By diagnosing the spectral HKE budget, the roles various physical processes play during the intensification of Lekima are investigated. Some processes, such as the buoyancy forcing, vertical transportation of IGWs, nonlinear interactions, and vertical transportation of convection, make remarkable contributions at almost all heights from lower troposphere to lower stratosphere. Other processes, by contrast, only make significant contributions at specific heights. The results for various budget terms here are qualitatively consistent with those in WZ2018, which both verifies the reasonability and authenticity of results from the idealized TC simulation and also suggests that the roles of various physical processes play during the intensification of TC are relatively steady.
The dependence of budget terms on scale and height is intriguing. For example, from the pressure vertical flux and its divergence, it can be found that upward-and downward-propagating IGWs with broad horizontal-wavelength range are excited near the tropopause. This phenomenon is consistent with the description of mechanical oscillator mechanism, which has been proposed to account for the generation of gravity waves excited by convective systems [36,37]. Meanwhile, strong generation sources of IGWs also exist in the upper troposphere and midtroposphere. Note that the two maximum centers of pressure vertical flux (Figure 8a) are right at the bottom of the two diabatic heating subranges ( Figure 2b); thus, it is reasonable to relate the generation of IGWs in the upper troposphere and midtroposphere with the forcing of heating [38]. Other interesting phenomena about IGWs are also worth noting: the specific height where the waves are generated near the tropopause is different for waves with different scales; for waves with wavelengths (wavenumbers) between 90 and 200 km (2π/9 × 10 −4 and π × 10 −5 rad m −1 ), the tropopause and upper troposphere are even not effective emission layers; and the upward-propagating waves in the lower troposphere are only remarkable for wavelengths (wavenumbers) slightly larger (smaller) than 100 km (2π × 10 −5 rad m −1 ). The selectivity of the generated IGWs on height and horizontal scale deserves further confirmation and study.
In addition, why does sporadic negative buoyancy forcing in the midtroposphere only appear at wavelengths (wavenumbers) around 100 and 400 km (2π × 10 −5 and π/2 × 10 −5 rad m −1 )? What makes the effect of nonlinear interactions so complicated in the troposphere for TC in the real atmosphere? And what determines the downward HKE vertical flux associated with the drag effect of precipitation particles at wavelengths (wavenumbers) around 100 km (2π × 10 −5 rad m −1 ) and below 60 km (above π/3 × 10 −4 rad m −1 )? These problems referring to selectivity on scale are also under explored.
For our simulation, the swelling of the HKE spectrum in the troposphere just extends to wavelengths (wavenumbers) around 100 km (2π × 10 −5 rad m −1 ); meanwhile, the vibrant spectral HKE budgets are restricted to wavelengths larger than 100 km. It is thus reasonable to connect the intensification of TC, the activation of energy activity at smaller scales, and the downscale extension of the arc-like spectral shape. The details and mechanisms are still unclear, pending further investigation.

Conflicts of Interest:
The authors declare no conflict of interest.