J. Meteor. Res.  2018, Vol. 32 Issue (1): 124-134 PDF
http://dx.doi.org/10.1007/s13351-018-7023-z
The Chinese Meteorological Society
0

#### Article Information

GUAN, Liang, and Xuyang GE, 2018.
How Does Tropical Cyclone Size Affect the Onset Timing of Secondary Eyewall Formation?. 2018.
J. Meteor. Res., 32(1): 124-134
http://dx.doi.org/10.1007/s13351-018-7023-z

### Article History

in final form October 18, 2017
How Does Tropical Cyclone Size Affect the Onset Timing of Secondary Eyewall Formation?
Liang GUAN1,2, Xuyang GE1
1. Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Key Laboratory of Meteorological Disaster, Nanjing University of Information Science & Technology, Nanjing 210044;
2. Shanghai Meteorological Science and Technology Service Center, Shanghai 200030
ABSTRACT: By using idealized numerical simulations, the impact of tropical cyclone size on secondary eyewall formation (SEF) is examined. Both unbalanced boundary layer and balanced processes are examined to reveal the underlying mechanism. The results show that a tropical cyclone (TC) with a larger initial size favors a quicker SEF and a larger outer eyewall. For a TC with a larger initial size, it will lead to a stronger surface entropy flux, and thus more active outer convection. Meanwhile, a greater inertial stability helps the conversion from diabatic heating to kinetic energy. Furthermore, the progressively broadening of the tangential wind field will induce significant boundary layer imbalances. This unbalanced boundary layer process results in a supergradient wind zone that acts as an important mechanism for triggering and maintaining deep convection. In short, different behaviors of balanced and unbalanced processes associated with the initial wind profile lead to different development rates of the secondary eyewall.
Key words: secondary eyewall formation     tropical cyclone size     surface entropy flux
1 Introduction

Tropical cyclone (TC) intensity forecast is one of the most challenging tasks for operational forecasting. TC intensity change is closely related to its structure change. For instance, a strong TC (e.g., the maximum wind speed greater than 60 m s–1) often has a secondary eyewall (SE) and then experiences eyewall replacement cycle (ERC). The SE, also called concentric eyewall (CE), is generally accompanied by ring-like convective clouds with a local maximum wind speed. The weak-echo area between the outer and primary inner eyewalls is called the moat region. As the SE is formed, the primary one starts to weaken and will be replaced by the outer one (i.e., the SE) eventually. During this ERC, the TC may experience dramatic intensity fluctuation and structure change (Willoughby et al., 1982; Black and Willoughby, 1992; Houze et al., 2006, 2007; Kuo et al., 2009; Sitkowski et al., 2011; Zhou and Wang, 2011, 2013; Zhou et al., 2011; Kossin and Sitkowski, 2012; Yang et al., 2013, 2014; Ge, 2015). Meanwhile, the convective SE formation results in damaging winds broadening and thus affects the extent of heavy rainfall and storm surge.

Recently, Ge et al. (2016) found that the onset timing of the SE has large variability. In their study, 45 typhoons over the western North Pacific were identified as CE events based on the Morphed Integrated Microwave Imagery (MIMIC) of 15-min intervals from the Cooperative Institute for Meteorological Satellite Studies (CIMSS)/University of Wisconsin–Madison. The SE formation (SEF) timing is defined as the period from the official warning time to the onset of SEF when the deep convection of the outer region covers at least two-thirds of a circle. Based on these TC samples, the mean timing of SEF is about 84 h, with standard deviation of 33 h. It is speculated that both internal dynamics and external forcing play vital roles in SEF (Willoughby et al., 1984; Nong and Emanuel, 2003; Kuo et al., 2004, 2008; Zhou and Wang, 2011; Huang et al., 2012; Rozoff et al., 2012; Abarca and Montgomery, 2013; Kepert, 2013; Qiu and Tan, 2013; Wang et al., 2013; Menelaou et al., 2014; Zhu and Zhu, 2014; Wang et al., 2016). Montgomery and Kallenbach (1997) addressed the roles of the vortex Rossby waves (VRWs) in the SEF. Nong and Emanuel (2003) suggested that SEF requires an external forcing. Terwey and Montgomery (2008) pointed out that the so-called “β-skirt” is necessary for SEF. Huang et al. (2012) suggested that the planetary boundary layer (PBL) imbalance serves as an important mechanism for triggering and supporting a ring of deep convection in a narrow supergradient wind zone around the ensuing outer eyewall. Qiu and Tan (2013) emphasized the role of asymmetric active outer spiral rainbands. Some studies focused on the impacts of environmental conditions on SEF (Kossin and Sitkowski, 2009; Zhou and Wang, 2013; Ge, 2015). For instance, moist environment favors a quicker SEF (Ge, 2015).

Ge et al. (2016) examined the sensitivity of SEF to TC initial size represented by the initial radius of maximum wind (RMW) and the outer radius of the gale-force wind (R34), and found that a TC with a larger initial size is apt to form SE (see their Fig. 1). They showed a relationship between the SEF timing and the initial radius of maxi-mum wind (RMW). The initial RMW is obtained from the best-track data of the Joint Typhoon Warning Center (JTWC). Specifically, there exists a significant negative correlation between the SEF timing and the initial RMW of TC; namely, the larger initial RMW the TC has, the quicker the SEF is. Nevertheless, no underlying mechanism was provided. Hence, the main purpose of this study is to reveal the possible physical processes involved.

The organization of this paper is as follows. In Section 2, the model and experimental design are introduced. An overview of the SEF simulation results is given in Section 3. Detailed analyses are conducted to give physical interpretations in Section 4. Conclusions and discussion are given in Section 5.

2 Model and experimental design

In this study, we use the ARW-WRF model (version 3.3.1) that has four nested domains with horizontal grids of 54, 18, 6, and 2 km, respectively. The Kain–Frisch convective parameterization scheme (Kain and Fritsch, 1993) is applied to the two outermost domains, and the microphysics scheme from Lin et al. (1983) is use for all the domains. All the simulations are conducted on an f-plane centered at 10°N in a quiescent environment with a constant sea surface temperature of 29°C. The longwave radiation scheme is the RRTM (A Rapid Radiative Transfer Model, Mlawer et al. 1997), and the shortwave radiation scheme is from Dudhia (1989). The YSU (Yonsei University) scheme (Hong et al., 2006) is used to calculate the PBL process.

The initial vortex is specified as a modified Rankine vortex. Specifically,

 $V t = { V max r RMW if r ≤ RMW V max ( RMW r ) b if r > RMW ,$ (1)

where r is the radius, b = 0.75 is the shape factor, and the initial maximum wind speed (Vmax) is set to 25 m s–1. By changing the RMW, five idealized vortices are tested with the RMW of 50, 70, 90, 110, and 130 km, respectively. The initial radial profiles of near-surface winds are presented in Fig. 1. For simplicity, these experiments are named as RMW50, RMW70, RMW90, RMW110, and RMW130, respectively. These vortices are simulated with a spin-up period of 10 days.

 Figure 1 Initial radial profiles of near-surface tangential wind (m s–1) in the five experiments.
3 Simulation results

To obtain an overall comparison of the evolution features, Fig. 2 displays the time–radius cross-section of the azimuthal-mean tangential wind at the height ofz = 1 km for all the experiments. We can see that the tangential wind gradually enhances, and the outer size broadens with time. A second maximum wind appears in certain outer region in each experiment. Once the outer entity forms, the inner one demises progressively, and is eventually replaced by the outer one. Of particular interest is that the SEF occurs earlier in the cases with larger initial RMWs. For example, the onset of SEF occurs the earliest at t = 107 h in RMW130. In contrast, the onset of SEF is much later (about t = 202 h) in RMW50, since its initial size is the smallest; that is, the development rate of SE is highly sensitive to TC initial size.

 Figure 2 The radius–time cross-section of azimuthally-averaged tangential velocities (m s–1) at the height of z = 1 km in the five experiments.

To further illustrate the salient differences in SEF, Fig. 3 presents the time–radius cross-section of azimuthally-averaged vertical velocity atz = 1 km. Initially, there is a maximum vertical velocity associated with the primary eyewall in the inner-core area. With time, an outer active convective ring emanates, which is collocated with the SE. Once the outer convective ring is established, the inner one weakens quickly and the outer one becomes dominated. In this regard, the simulations capture the major evolution characteristics of convection during an ERC. The timing of the establishment of the active outer convection reflects the development rate of SE. Hence, an earlier prolific convective updraft indicates a quicker SEF in the cases with a larger initial size.

 Figure 3 As in Fig. 2, except for vertical velocity (m s–1) at the height of z = 1 km in the five experiments.

Figure 4 presents the sequences of simulated radar reflectivity at the height of z = 1 km during the SEF, using RMW50 and RMW130 as examples. In this study, the inner and outer eyewall radii are referred to as the mean radius of convective ring, and the moat width is represented by the distance between the two eyewalls. For RMW130, the inner eyewall is located at the radius of about 40 km around t = 102 h, surrounded by outer ring-like convection. Thereafter, the outer convection contracts gradually and evolves into the SE at the radius of 120 km. A clear moat region can be identified during this period of interest. Once SE starts to emerge, TC intensity begins to decrease (figure omitted). The SEF in RMW50 bears many similarities as that in RMW130, except that the SEF occurs much later (i.e., t = 202 h). In addition, there are some remarkable differences in the structure of SE. For example, the moat size shows a positive correlation with the initial TC size. More specifically, the moat is much wider (about 80 km) in RMW130 than in RMW50 (about 40 km). Furthermore, the intensity fluctuation is greater in RMW130 during the ERC, which is consistent with Zhou and Wang (2013). They found that a TC with a larger SE may experience more significant intensity fluctuation. Since our focus is to examine the development speed of SE, the structure and intensity changes during the ERC are not covered in the current work.

 Figure 4 Evolution of simulated radar reflectivity in (a, b, c) RMW50 and (d, e, f) RMW130 during the SEF at a time interval of 4 h.

To summarize, the onset timing of SEF exhibits remarkable sensitivity to initial TC size, indicating that the development rate of the SE is highly sensitive to TC structure. A TC with a larger initial size is apt to create a SE more quickly. This supports the observational results of Ge et al. (2016). In addition, the SE is apt to emanate farther away from the center. Naturally, a question arises—How does the initial TC size influence the sequential SEF? We will address this question in the next section.

4 Possible mechanisms

Physically, TC size reflects the radial distribution of surface wind speed, which directly determines surface heat fluxes. Surface heat fluxes can modulate atmosphe-ric moisture supplies, and thus convective activities. In this regard, the differences in surface heat fluxes between RMW50 and RMW130 should be compared. Figure 5 presents the time–radius cross-sections of azimuthally-averaged surface heat fluxes. As anticipated, after 96 h of integration, pronounced surface heat fluxes are evident in the region outside of the primary eyewall (i.e., around the radius of 120 km) in RMW130. On the contrary, the surface heat fluxes in the outer region in RMW50 increase much more slowly untilt = 192 h. One remaining issue here is the role of the heat fluxes from the lateral boundary. As shown in Fig. 5, beyond the radius of 180 km, the differences in heat fluxes in both cases are smaller than those in the inner-core area. It is likely that the local surface heat fluxes play a dominant role. However, Gao et al. (2017) pointed out that total precipitable water in the outer environment, where moisture is mainly provided by surface evaporation, is also vital to storm’s rapid intensification (RI) because RI is favored when there is less dry air intruded into the storm circulation. Ge (2015) also found that a TC is apt to create SEF more quickly under a moist environment. Both studies suggested the important impacts of lateral boundary or environmental condition, which awaits further studies.

 Figure 5 The radius–time cross-section of azimuthally-averaged surface heat fluxes (W m–2) in (a) RMW50 and (b) RMW130.

Accordingly, stronger surface entropy fluxes in the outer region favor greater convective available potential energy (CAPE) and thus active convection. Toward this end, Fig. 6 displays the evolution of simulated diabatic heating averaged between 4 and 8 km in RMW50 and RMW130. Obviously, the onset timing of strong diaba-tic heating is consistent with that of surface heat flux. Hence, different initial wind profiles can likely account for distinct radial distributions of surface heat fluxes and convection.

 Figure 6 The radius–time cross-section of azimuthally-averaged diabatic heating (× 10–3 K s–1) averaged between the heights of 4 and 8 km in (a) RMW50 and (b) RMW130.

In general, the latent heating released in the outer region helps enhance the low-level radial inflows and accelerate tangential winds outside the eyewall, leading to a broadening of tangential wind field. To demonstrate this, the radius–time cross-section of azimuthally-averaged radial wind field at the low level is given inFig. 7. We can see that marked enhancement of radial inflows occurs in the outer region. In turn, the radial inflow helps accelerate the tangential wind.

 Figure 7 The radius–time cross-section of azimuthally-averaged radial winds (m s–1) at the height of z = 1 km in (a) RMW50 and (b) RMW130.

A greater tangential wind reflects a larger inertial parameter. From the balanced dynamics viewpoint, a greater inertial stability helps the conversion from diaba-tic heating to kinetic energy, and favors local secondary maximum wind (Shapiro and Willoughby, 1982; Rozoff et al., 2012; Sun et al., 2013). In this regard, the evolution of inertial stability should be compared. Figure 8 presents the radial distribution of inertial stability averaged during 104–116 h at the height of z = 4 km. Essentially, the inertial stability is larger in the inner-core area. As the radius increases, the stability decreases rapidly. This tendency is more pronounced in RMW50 than in RMW130. As such, the inertial stability becomes larger beyond the radius of 40 km in RMW130.

 Figure 8 Time-averaged azimuthal-mean inertial stability (× 10–4 s–1) at the height of z = 4 km in RMW50 (solid line) and RMW130 (line with circles) during 104–116 h.

A budget analysis of the azimuthal-mean tangential wind tendency (Wang et al., 2016) is conducted. The budget equation can be approximated by:

 $∂ v t ¯ ∂ t = − v r ¯ ζ a ¯ − w ¯ ∂ v t ¯ ∂ z − v ′ r ζ ′ a ¯ − w ′ ∂ v ′ t ∂ z ¯ + F ¯ ,$ (2)

where vt, vr, and w are the tangential, radial, and vertical velocities, respectively; and ζa is the absolute vorticity. The overbar denotes the azimuthal mean and the prime is the deviation from the azimuthal mean. The terms on the right-hand side are the mean radial advection, mean vertical advection, eddy radial advection, eddy vertical advection, and friction, respectively. The analysis shows that the mean radial advection term ( $− v r ¯ ζ a ¯$ ) is the key contributor to the enhancement of tangential wind in the boundary layer. As shown in Fig. 9, this effect is pronounced in the outer region (i.e., beyond the radius of 100 km) in RMW130, which agrees well with a faster enhancement of tangential wind therein.

 Figure 9 Radial distribution of the mean radial advection term (m s–2) averaged during 104–116 h at the height of z = 1 km.