Dynamic Mechanisms Associated with the Structure and Evolution of Roll Vortices and Coherent Turbulence in the Hurricane Boundary Layer: A Large Eddy Simulation During the Landfall of Hurricane Harvey

Roll vortices are a series of large-scale turbulent eddies that nearly align with the mean wind direction and prevail in the hurricane boundary layer. In this study, the one-way nested WRF-LES model simulation results from Li et al. (J Atmos Sci 78(6):1847–1867, https://doi.org/10.1175/JAS-D-20-0270.1, 2021) are used to examine the structure and generation mechanism of roll vortices and associated coherent turbulence in the hurricane boundary layer during the landfall of Hurricane Harvey from 00 UTC 25 to 18 UTC 27 August 2017. Results indicate that roll vortices prevail in the hurricane boundary layer. The intense roll vortices and associated large turbulent eddies above them (at a height of ~ 200 to 3000 m) accumulate within a hurricane radius of 20–40 km. Their intensity is proportional to hurricane intensity during the simulation period. Before and during hurricane landfall, strong inflow convergence leads to horizontal advection of roll vortices throughout the entire hurricane boundary layer. Combined with the strong wind shear, the strongest roll vortices and associated large turbulent eddies are generated near the eyewall with suitable thermodynamic (Richardson number at around − 0.2 to 0.2) and dynamic conditions (strong negative inflow wind shear). After landfall, the decayed inflow weakens the inflow convergence and quickly reduces the strong roll vortices and associated large turbulent eddies. Diagnosis of vertical turbulent kinetic energy indicates that atmospheric pressure perturbation, caused by horizontal convergence, transfers the horizontal component of turbulence to the vertical component with a mean wavelength of about 1 km. The buoyancy term is weak and negative, and the large turbulent eddies are suppressed.


Introduction
Roll vortices are a series of large-scale turbulent eddies (Etling and Brown 1993;Mourad and Walter 1996;Weckwerth et al. 1997;Young et al. 2002). They can be resident in the hurricane boundary layer (Wurman and Winslow 1998;Katsaros et al. 2000Katsaros et al. , 2002Huang et al. 2018). They can also enhance momentum, energy, and composition exchange in the atmospheric boundary layer (Fujita 1992;Kepert 1996;Makin 1998;Lorsolo et al. 2008;Zhang et al. 2008;Ellis and Businger 2010;Gao et al. 2017). Improved understanding of the mechanism of roll vortices can enhance comprehension of the planetary boundary layer (PBL) and consequently atmospheric predictions, especially for intensely damaging hurricanes (Ernst et al. 2019).
Previous studies have found that roll vortices in the hurricane boundary layer can nearly align with or be significantly different from the mean wind, with wavelength varying from 200 m to over 3000 m (Gall et al. 1998;Wurman and Winslow 1998;Lorsolo et al. 2008;Morrison et al. 2005;Nolan 2005). Morrison et al. (2005) verified from radar observations that roll vortices are oriented approximately 4°to the left of the mean wind and extend from 200 to 800 m above the ground in the hurricane boundary layer. Through radar and satellite observations, Huang et al. (2018) found that the wavelength of a roll in the marine atmospheric boundary layer is shortest near the storm centre and increases and then decreases with distance from the storm centre. Moreover, compared to numerical flux simulations with the PBL scheme in mesoscale models (e.g., Hong and Pan 1996;Hong et al. 2006;Hong 2010), some of the fluxes associated with roll vortices can be about 2-3 times greater (Morrison et al. 2005;Zhu 2008b; Gao and Ginis 2016). Strong fluxes induced by roll vortices change the hurricane structure and finally impact hurricane damage through intense surface winds (Saffir 1973;Wakimoto and Black 1994).
While many previous studies have emphasized the structure of roll vortices and their influence on hurricane structure and evolution, the dynamic mechanism that causes the genesis of hurricane roll vortices has remained under continuing investigation. LeMone (1973) found that rolls in the PBL are maintained primarily by the energy of the buoyancy and/or smaller turbulent eddies that reinforce the rolls. A theoretical study (Foster 2005) and several idealized numerical simulations (Gao and Ginis 2014) have found that hurricane roll vortices are caused by inflow inflection point instability, namely dynamic instability near the height with the sharpest inflow change of the profiles, similar to roll vortices in the Ekman layer (Lilly 1966;Brown 1980). Young et al. (2002) found that the different instability of wind components can lead to a different kind of roll. Foster (2005) used a theoretical study of roll vortices to verify that aligned with the mean wind, roll vortices are dominated by the inflection point of radial wind. Moreover, in the hurricane boundary layer, the dynamic instability with wind shear generates a quasi-equilibrium roll near neutral conditions, once the divergence of the roll fluxes transfers kinetic and/or potential energy from the perturbation back into the mean flow at the same rate at which it is extracted (LeMone 1973;Brown 1980). Specifically, strong buoyancy forcing inhibits the rolls, and they can be generated only when the flux Richardson number (Rf) is less than 0.25. With an idealized simulation, Gao and Ginis (2014) further confirmed that the rate of generation of roll vortices (w ) is positively correlated with the strength of the inflow shear there. However, their finding has not yet been confirmed with a real hurricane due to the lack of high-resolution observations.
Fortunately, numerical simulations can provide insights into the structure of roll vortices inside a hurricane (Nakanishi and Niino 2012;Gao and Ginis 2014. Large eddy simulation (LES), which is commonly used to simulate turbulence in the atmospheric boundary layer, can also be used to simulate roll vortices in the hurricane boundary layer (Nakanishi and Niino 2012;Wang and Jiang 2017). Since LES was introduced into the Weather Research and Forecasting Model (WRF) (Skamarock et al. 2019), large turbulent eddies, including roll vortices, have been simulated in real hurricanes by LES through nested grids within an advanced research version of the community Weather Research and Forecasting (WRF-LES) model (Zhu 2008a;Zhu et al. 2016). With LES, the strong fluxes associated with large eddies can be simulated without using a PBL parameterization scheme (Zhu 2008b). In a recent study,  used the WRF-LES model to simulate Hurricane Harvey (2017) around its landfall and analyzed the impact of large turbulent eddies on the evolution of the hurricane before, during, and after landfall. They found that roll vortices were reasonably simulated by WRF-LES. Their diagnoses also found that large turbulent eddies, including roll vortices, could generate a strong flux near the hurricane eyewall and influence the wind structure and precipitation of the hurricane during its landfall.
Following , in this study we further evaluate the dynamic mechanisms associated with the structure and evolution of roll vortices and other large turbulent eddies in the hurricane boundary layer, with the goal of enhancing understanding of mechanisms that lead to the generation and decay of large eddies in a real hurricane. The WRF-LES model setup and turbulence data processing are briefly described in Sect. 2. Validation of coherent eddies, including roll vortices and other large turbulent eddies, is discussed in Sect. 3. A structural analysis and genesis mechanism are provided in Sect. 4. Diagnoses of the maintenance mechanism of large turbulent eddies are provided in Sect. 5. Conclusions are presented in Sect. 6.

The WRF-LES Model
As described in , the one-way nested WRF model (version 3.9.1.1) with LES was used to simulate Hurricane Harvey (2017) during its landfall from 06 UTC 25 to 18 UTC 27 August 2017. In the model, the US National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) final analysis (FNL) 0.25°× 0.25°data were used to derive the initial and boundary conditions. The YSU PBL scheme (Hong et al. 2006) was used in the first two domains, while the inner two domains used LES directly with the TKE 1.5 closure sub-scale turbulence scheme. These four domains were configured with grid meshes (horizontal resolution) of 150 × 150 (12.5 km), 251 × 281 (2.5 km), 951 × 1131 (0.5 km), and 1401 × 1401 (0.1 km), respectively. For the 71 vertical levels, over 39 levels were below 1500 m, with a minimum height of approximately 3 m. Domain 4 was moved with time to follow the hurricane. Configuration details can be found in .

Derivation of the Turbulence Field
As described in , the turbulence in LES is derived from the simulation result by subtracting the smoothed fields produced from the horizontal two-dimensional Gaussian filter:  Figure 1a and d shows the 10 m wind speed from the original simulation conducted with LES and the simulation results after smoothing by the Gaussian filter, respectively, at 01 UTC 26 August 2017. The smoothed field (Fig. 1d) retains the hurricane structure and smooths out the perturbations. Therefore, the turbulence field can be computed by subtracting the smoothed field (e.g., Fig. 1d) from the original simulated field (e.g., Fig. 1a). In addition, the sensitivity of filtering results to the filter region was examined by varying the filter region to 50 × 50 grid. The results (not shown) are similar to those in Fig. 1, confirming the usefulness of the filtering method described in Eq. (1). Hereafter, we use these turbulence fields produced by subtracting the smoothed field generated by Eq. (1) from the original simulations to analyze the turbulence structure in LES.  found that turbulence was reasonably resolved with roll vortices in domain 4 of the WRF-LES simulation. The turbulence in LES leads to strong wind perturbations, resulting in a higher surface wind that is consistent with the observations in Wurman and Kosiba (2018) and Fernández-Cabán et al. (2019) for Hurricane Harvey. In this study, we focus on domain 4 to analyze the roll vortices and other eddies.

Structure of Turbulence Fields
Figures 1b, c and 2a-d show the 10 m turbulence of w (Fig. 1b, c), u r (inflow, Fig. 2a, c), and u t (tangential wind, Fig. 2b, d) at 01 UTC 26 August 2017. Compared to Fig. 1d, many perturbations throughout the entire storm nearly align with the mean wind. A partially enlarged view in Figs. 1c and 2c, d shows that these perturbations are associated with updraft and downdraft turbulence that aligns along the mean wind direction and rotationally downstream of the storm center. The coupled positive and negative inflow and tangential turbulence along the mean wind direction combined with the coupled updraft and downdraft turbulence organize a spiraling turbulent eddy extending along the mean wind direction to the storm center.
To reveal the turbulence vertical structure, Fig. 1c, e, and f show the w' fields at different heights from 10 to 1000 m at 01 UTC 26 August 2017. With increasing height from 10 to 100 m, w is gradually enhanced in the hurricane boundary layer and the turbulent eddy always nearly aligns with the mean wind direction. When the radius is reduced to the storm center, spiraling turbulent eddies throughout the storm boundary layer gradually converge along the storm flow field, and the strong convergence makes the well-structured spiraling turbulent eddy relatively chaotic near the eyewall. At a height of 1000 m, the well-structured spiraling turbulent eddy nearly disappears, but there is a strong chaotic turbulent eddy near the eyewall. There is no clear relationship between the rainband (dashed contour line) and the turbulent eddy. The distribution of rainband did not affect the distribution of strong turbulence in the hurricane. Hurricane observations from Guimond et al. (2018) reveal that strong turbulence exists near the eyewall and can extend upward to over 1000 m. The simulation here produces  Morrison et al. (2005), indicating quasi-steady-state turbulent eddies, namely, roll vortices (rolls hereafter). These rolls can extend from near the surface to 200-400 m. Further, near the eyewall (radius of 20-40 km), there are strong perturbations above the rolls extending from approximately 200 m to over 1000 m, even to around 3000 m. These large turbulent eddies (eddies hereafter) are associated with rolls but are different from previous roll observations ( Morrison et al. 2005) and simulations Ginis 2014, 2016), which focused on rolls relatively far away from the storm center. However, these eddies are very similar to the turbulent eddies observed near the hurricane eyewall (Guimond et al. 2018). We found that these strong eddies can connect to the regular rolls near the top of the roll. This connection can also be found away from the storm center but usually has weaker eddy than that near the eyewall. LeMone (1973) suggested that mesoscale turbulence or structure can interact with rolls. Therefore, in the LES simulation presented in this study, eddies (200-1000 m) are associated with rolls (from near the surface to 200 m) and form coherent eddies (including rolls and eddies) near the eyewall from near the surface to 1000 m and even beyond (3 000 m).

Wavelength
The wavelength varies in different turbulent structures (Young et al. 2002;Foster 2005). To distinguish the differences and examine the connection between eddies and rolls, the wavelengths of strong eddies and rolls are compared. As shown in Fig. 1f and 2e-g, strong eddies occur mainly within a 20-40 km radius of the storm center and extend upward from approximately 200 m to over 1000 m. We examined four cross-section lines from the hurricane center to the north, east, south, and west at heights of 100 m and 1000 m during the simulation period to obtain roll and eddy perturbations, respectively. The w is considered an eddy at a height of 1000 m and an roll at a height of 100 m over an area within a radius of 20-40 km.
Since rolls and eddies both have coupled upward and download turbulence motion, the wavelength is detected by the distance between two adjacent local maxima or minima of the cross-section data. Figure 3a shows an example of the probability distribution of the wavelength of rolls within a radius of 20-40 km and 40-100 km at a height of 100 m during hurricane landfall at 00-02 UTC on 26 August. The wavelengths based on local maxima and minima are both shown in the figure. The probability distributions of wavelengths at different radii from the local maxima and minima are close to each other, with wavelength maxima at around 0.5-1.5 km and a mean wavelength of 0.8 km. This wavelength is slightly different from the mean wavelength of about 1.6 km found by Morrison et al. (2005) but is the same as that in Huang et al. (2018), who observed an approximately 1 km wavelength roll in the same region through radar and satellite observations. The similar wavelength probability distributions of rolls at different radii, i.e., 20-40 km and 40-100 km, indicate that the structure of rolls is the same near the eyewall and away from the eyewall. Furthermore, the distribution of wavelength of rolls did not change (change less than 1%) for different filter region choices (not shown), indicating that the scale of the rolls is conserved with respect to filter region. With a depth of approximately 200-400 m ( Fig. 2e-g) and a mean wavelength of The global mean, median, and number (n) are also shown 0.8 km, the roll aspect of 2-4 confirms that the LES-simulated rolls are close to the summary by Young et al. (2002), who found that roll aspect varied from 2 to 6 with a PBL height of around 200 m.
For eddies, Fig. 3b shows an example of the probability distribution of wavelengths within a radius of 20-40 km at a height of 1000 m during hurricane landfall at 00-02 UTC 26 August. The probability distributions of wavelengths from local maxima and minima are close to each other, with global maxima at around 0.5-1.5 km and a mean wavelength of 0.9-1.1 km, which is very close to the observations of Guimond et al. (2018). Compared to the rolls, the wavelength probability distributions of eddies have a slightly lower percentage in the 0.5-1 km wavelength range (10% to eddies and 15% to rolls). Apart from this difference, the similar probability distributions of wavelengths between eddies and rolls indicate the connection between eddies and rolls near the eyewall.
To understand the generation and evolution of these coherent eddies, including both rolls and eddies, in the following sections, we examine their structure and mechanism.

Evolution of Roll Vortices and Turbulence
To examine the evolution of rolls and eddies, Fig. 4 shows the azimuthally averaged vertical turbulence intensity (namely, Iw w 2 /2, contour line, w is the turbulence vertical component) before Hurricane Harvey's landfall at 12 UTC (  During the simulation period, although rolls prevail over the entire hurricane boundary layer (Figs. 1 and 2), rolls are weak away from the eyewall, which makes azimuthally averaged Iw usually less than 0.5 m 2 s −2 . The Iw near the eyewall is proportional to hurricane intensity (MWS), while the Iw distribution along the radius is not proportional to wind speed. Along the radius, wind speed can be reduced to only half the maximum at a radius of 70-80 km, while Iw can be reduced to less than one-eighth of the maximum Iw.
Horizontal turbulence is also analyzed in Fig. 4, which shows the azimuthally averaged horizontal turbulence intensity (namely, Iu (u r 2 + u t 2 )/2, shaded; u r is the turbulence radial component; u t is the turbulence tangential component) before (Fig. 4a-c), during ( Fig. 4d-f), and after hurricane landfall (Fig. 4g-i). Before landfall, under strong Iw, there is always a strong Iu associated with rolls at heights below 200 m, with a maximum of over 15 m 2 s −2 below 200 m. Within the region of strong Iw, eddy-associated Iu is relatively weak, with a maximum of over 3 m 2 s −2 . During landfall, Iu below 200 m gradually reduces from a maximum of over 33-27 m 2 s −2 with reduced MWS. Above 200 m, eddy-associated Iu also gradually decreases, with a maximum of over 21-5 m 2 s −2 . After landfall, below 200 m, there is still horizontal turbulence, but it is also weakened, with maximum Iu decreasing from 13 to 9 m 2 s −2 . Above 200 m, horizontal turbulence is the same as vertical turbulence and almost disappears, with maximum Iu less than 3 m 2 s −2 .
Similar to vertical turbulence, horizontal turbulence near the eyewall, including turbulence contributed by eddies and rolls, is also proportional to hurricane intensity, while the decrease in Iu along the radius is significantly quicker than the decrease in wind speed. Although rolls usually extract energy from the mean flow (LeMone 1973;Foster 2005), the reduction of Iu and Iw along the radius is always quicker than that of wind speed in the simulation result. Therefore, the extremely intense coherent eddies (rolls and eddies) near the hurricane eyewall are not caused only by the strong wind there. The specific stream structure of the hurricane-namely, accumulation by horizontal convergence-may also be critical for extremely intense turbulent kinetic energy (TKE) near the eyewall. Meanwhile, roll-related Iu (below 200 m) is usually larger than that related to eddies, while roll-related Iw (below 200 m) is usually weaker than that related to eddies, indicating a different turbulent kinetic energy distribution of eddies and rolls.

Stability Conditions
Previous studies have found that rolls are often generated with a flux Richardson number (Rf) less than 0.25 (Gao and Ginis 2014). Near-neutral conditions can suppress a roll enhanced by convection (Foster 2005). To clarify the stability of thermodynamic conditions associated with rolls and eddies, azimuthally averaged Rf ( g θ v w θ v /( u w ∂ u ∂z + v w ∂ v ∂z )) at the same time as in Fig. 4 is shown in Fig. 5 before (Fig. 5a-c), during ( Fig. 5d-f), and after hurricane landfall (Fig. 5g-i). The height of the thermodynamic hurricane boundary layer is determined by the height where the potential temperature gradient is over 3 K km −1 with a temperature inversion above to create a potential temperature difference of over 2 K (Heffter 1980). This calculated hurricane boundary layer height (HBLH) is also azimuthally averaged and added in Fig. 5. For convenience, azimuthally averaged Iw, representing coherent eddy rolls and eddies, and Iu (dashed line), which majorly representing rolls, are also displayed in Fig. 5 as contour lines.
From Fig. 5, we found that before landfall, rolls and eddies are generated mainly with Rf at − 0.2 to 0.2. Only a small region has rolls and eddies with Rf exceeding this threshold. The HBLH always has a lower value of about 200 m in the region where strong rolls and eddies prevail. The main part of eddies is always above the HBLH. During landfall, rolls and eddies gradually decay where Rf gradually exceeds the threshold of − 0.2 to 0.2. The HBLH also has a lower value in the roll-prevailing region except at 09 UTC 26 August, when a lower HBLH is located near the storm center. After landfall, although the strong rolls and eddies almost disappear, some regions still have Rf within the threshold of − 0.2 to 0.2. The HBLH is almost unchanged at about 200 m, except in the eye center. During the simulation period, there are suitable thermodynamic conditions for rolls with Rf within − 0.2 to 0.2. After landfall, these suitable conditions still exist, but strong rolls and eddies disappear. Results here indicate that thermodynamic conditions are necessary but not sufficient for the generation of strong rolls and eddies near the hurricane eyewall.  Foster (2005) pointed out that a dynamic condition-namely, radial wind inflection point instability in the hurricane boundary layer-is fundamental for the generation of rolls. To investigate this dynamic condition, the azimuthally averaged radial wind shear (positive for inflow) at the same time as in Fig. 4 is shown in Fig. 6 before (Fig. 6a-c), during (Fig. 6d-f), and after hurricane landfall (Fig. 6g-i). Azimuthally averaged Iw and Iu are also added to the figure to distinguish the eddies and rolls near the eyewall. We see that during the whole simulation period, strong negative inflow wind shear covers the main part of rolls (height of 30-200 m). The negative inflow wind shear is distributed all around the storm boundary layer and is consistent with the prevailing rolls shown in Figs. 1 and 2. The strong negative inflow wind shear is related to the strong dynamic instability. This inflow wind shear gradually enhances from an intensity of − 0.08 to over − 0.12 s −1 with hurricane intensification. During hurricane landfall, within a radius of 0-20 km in the hurricane eye region, there is also a strong negative inflow wind shear, with a minimum of less than − 0.12 s −1 . This strong  Fig. 5 in this strong negative inflow wind shear region, very stable or unstable stratification suppresses the generation of rolls there (Foster 2005). After landfall, the strong rolls almost disappear even though there is still negative inflow wind shear with a minimum of less than − 0.08 s −1 .
The inflow wind shear is the genesis of the hurricane roll (Foster 2005;Gao and Ginis 2014). However, the dynamic instability and thermodynamic conditions here cannot explain the large Iw associated with extremely intense rolls and eddies that exist near the hurricane eyewall but explain only the prevailing rolls in the entire hurricane boundary layer (Figs. 1  and 2). Other mechanisms of the hurricane, e.g., accumulation by horizontal convergence, may be critical for extremely intense TKE near the eyewall.

Convergence of Inflow
With strong inflow in the hurricane boundary layer, mass convergence and accumulation near the eyewall are observed. The impact of inflow convergence near the eyewall is analyzed in Fig. 7, which shows the horizontal divergence of the azimuthally averaged inflow (du r /dr; u r represents the inflow while the positive values for inflow) at the same time as in Fig. 4 before (Fig. 7a-c), during (Fig. 7d-f), and after hurricane landfall (Fig. 7g-i). The r-axis points to the center of the hurricane. Similarly, the azimuthally averaged Iw and Iu are also added to Fig. 7 as strong eddies and rolls with contour lines. From the figures, we see that there is always a strong negative du r /dr covering the main part of the strong rolls below 200 m with a minimum of less than − 2·× 10 -3 s −1 . This strong negative du r /dr indicates strong inflow convergence near the eyewall. During landfall, the strong inflow convergence region with a minimum of less than − 2 × 10 -3 ·s −1 still covers the main part of the strong rolls near the eyewall. After landfall, the convergence is weak, with a minimum of more than − Fig. 7 Similar to Fig. 4 but for the horizontal divergence of azimuthally averaged inflow (du r /dr). The azimuthally averaged Iw (solid line, contour interval 0.5 m 2 s −2 ) and Iu (dashed line, contour interval 4 m 2 s −2 ) are added as a contour line 1.2 × 10 -3 ·s −1 , related to the near disappearance of the strong eddies and rolls. Therefore, inflow convergence near the eyewall is critical for the strong eddies and rolls in the hurricane. Horizontal convergence is significantly related to the horizontal advection effect, and there is clearly roll horizontal advection through the horizontal wind from some distance away to near the eyewall, as shown in Figs. 1 and 2a-d. The inflow transports the rolls in the entire storm boundary layer and accumulates near the eyewall to support the intense rolls there. Once the horizontal transport, especially the inflow, gradually decays after hurricane landfall, the TKE accumulation mechanism is also gradually destroyed.
Based on this relationship of horizontal convergence and coherent eddies,  parameterized hurricane turbulence near the eyewall with the WRF model through convergence and added the vertical mixing contributions of rolls and eddies in a modified PBL scheme. The better hurricane forecast from the WRF model with the modified PBL scheme confirmed that horizontal convergence near the eyewall is critical for forming the intense rolls and eddies in a hurricane (see details in Li and Pu 2021).

Maintenance of the Vertical Motion of Rolls and Eddies
Although we have known that inflow convergence transports roll turbulence to near a storm's eyewall, it is unclear how this horizontal advection supports the vertical motion of turbulence near the eyewall, especially for the strong vertical component of eddies, which can strongly change vertical mixing near the eyewall (Zhang and Drennan 2012;Zhao et al. 2020;. To this end, Table 1 shows a summary of the magnitude of the wind shear term, buoyancy term, and pressure transport term, all of which can contribute to the TKE, during the simulation based on the following equations:  (3) On the right-hand side, terms 1 and 2 represent horizontal wind shear (in Eqs. 2 and 3) and vertical wind shear (in Eq. 4), term 3 and 4 in Eqs. 2-4 represent transport contributions that are very weak and neglected in this study. The last term represents the dissipation in Eqs. 2-4. In Eqs. 2-4, only a few terms contribute to the vertical motion of turbulence. The contribution details by the wind shear term, the buoyancy term, and the pressure transport term that are discussed in this section to verify the energy support for the vertical motion of rolls and eddies. Note that the equation should consider the contribution of both mesoscale and mean winds. In this study, the result from the Gaussian filter shows a very small change (less than 1%) to the detrend result, indicating the remove of the mesoscale impact. Therefore, the mesoscale contribution is neglected in Eq. 2-4 and in the following discussion. Figure 8 shows the azimuthally averaged horizontal wind shear term at the same time as in Fig. 4 before (Fig. 8a-c), during (Fig. 8d-f), and after hurricane landfall (Fig. 8g-i). Before landfall, in the region of strong rolls, the wind shear term is positive and gradually enhances, from a maximum of over 0.2-0.4 m 2 s −3 . This positive wind shear term supports the strong rolls in the hurricane boundary layer (Foster 2005;Gao and Ginis 2014). In the region of strong eddies, the wind shear term can only reach a maximum of over 0.2 m 2 s −3 . The maximum of 0.4 m 2 s −3 (0.2 m 2 s −3 ) is considered to be the magnitude of the contribution of the horizontal wind shear term to rolls (eddies), as summarized in Table 1. During landfall, the wind shear term is still positive around the region of strong rolls and eddies and gradually weakens. After landfall, similar to the strong eddies, wind shear is weak, with the maximum gradually decreasing to less than 0.1 m 2 s −3 .

Energy Support
Although the horizontal wind shear is strong, it does not directly generate vertical velocity associated with rolls and eddies (Fig. 4) (Zhou et al. 2019). The vertical wind shear term (not shown), which supports the vertical component of TKE associated with rolls and eddies (Fig. 4), is extremely weak, usually negative, and less than 0.001 m 2 s −3 (magnitude of the vertical wind shear term contribution shown in Table 1). Therefore, the wind shear does not support the vertical component of TKE associated with strong rolls and eddies.
To clarify the energy support for the vertical component of TKE, Fig. 9 shows the azimuthally averaged buoyancy term g < w θ L > /θ L at the same time as in Fig. 4 before (Fig. 9a-c), during (Fig. 9d-f), and after hurricane landfall (Fig. 9g-i). Effects of the water vapor, cloud water, and rain water are considered in the buoyancy term. Before landfall, around the region of strong eddies, the buoyancy term g < w θ L > /θ L is negative and gradually enhances, from a minimum of less than − 0.003 m 2 s −3 to less than − 0.008 m 2 s −3 (the  Table 1). This negative buoyancy term reflects a suppression effect on the vertical motion of strong eddies in the hurricane boundary layer. In the region of rolls, the buoyancy term g < w θ L > /θ L is very weak, with a minimum of − 0.001 m 2 s −3 (the magnitude of the buoyancy term contribution to rolls is shown in Table 1). During landfall, g < w θ L > /θ L is still negative around the region of strong eddies and gradually weakens. After landfall, similar to the strong eddies, g < w θ L > /θ L is weak, with a minimum of less than − 0.003 m 2 s −3 . Similar to rolls (Foster 2005;Lilly 1966;Brown 1980), the buoyancy term g < w θ L > /θ L is not a significant factor for eddies, as its value is significantly less than that of the wind shear term. The term g < w θ L > /θ L only slightly suppresses the generation of eddies.

Energy Transfer
Since the buoyancy term does not support the strong vertical motion of rolls and eddies, other transport terms in the TKE equation should be concerned with supporting the vertical  Figure 10 shows the azimuthally averaged pressure transport term d < w P > /ρdz at the same time as in Fig. 4 before (Fig. 10a-c), during (Fig. 10d-f), and after hurricane landfall (Fig. 10g-i). Before landfall, there is a strong pressure transport term in the same region as large Iw. With hurricane intensification, the pressure transport term is enhanced, with a maximum of over 0.4 m 2 s −3 or a minimum of less than − 0.4 m 2 s −3 , which enhances upward and downward turbulence, respectively. During landfall, the transport term is still strong, with a maximum of over 0.4 m 2 s −3 or a minimum of less than − 0.4 m 2 s −3 in the same region as large Iw, indicating a contribution of the pressure transport term of ± 0.4 m 2 s −3 (shown in Table 1). After landfall, the pressure transport term is weak but still exists, with a maximum of over 0.2 m 2 s −3 or a minimum of less than − 0.2 m 2 s −3 . The strong eddies almost disappear.
Compared to the weak buoyancy term (one-tenth of a percent of the pressure transport or less), the pressure transport term dominates the vertical component of rolls and eddies. Since the pressure transport term can only transfer the horizontal component of turbulence to the vertical component, the quick decay of strong eddies after landfall cannot be attributed to  Fig. 4 but for the azimuthally averaged pressure transport. The azimuthally averaged Iw is added as a contour line the energy transport of the pressure transport term but to the decay of the horizontal wind shear term (Fig. 8g-i). Before and during landfall, the pressure transport term transfers the horizontal component of rolls and eddies to support the strong Iw at a height of 75-3000 m. The energy transfer through the pressure transport term, which covers the region of rolls and eddies, connects the vertical component of rolls and eddies, resulting in a similar wavelength structure shown in Figs. 3a, b.
From the definition, pressure transport in the TKE budget is associated with the vertical perturbation pressure gradient force and then controlled by pressure perturbation. The pressure is related mainly to air density and temperature through ideal gas law. To verify the contributions of these two factors, the indexes P ρ ρ and P T T represent the contributions of air density and air temperature perturbations to the pressure perturbation, are analyzed in Fig. 11. The larger index indicates a larger impact on pressure perturbation from air density and/or temperature. Figure 11 shows an example of the azimuthally averaged P ρ ρ and P T T during hurricane landfall at 0100 UTC on 26 August 2017. P ρ ρ is always positive, with region are both smaller than that of P ρ ρ . This reflects the positive and strong relationship between pressure perturbation and air density, and the weak impact of air temperature perturbation on pressure perturbation. As shown in Fig. 7, there is usually strong convergence near the eyewall. The density perturbation may be associated with this strong convergence there. Therefore, due to the strong horizontal convergence in the storm, there is some density perturbation that leads to pressure perturbation and finally mixes the vertical perturbation to enhance the Iw associated with rolls and eddies. The horizontal advection of rolls throughout the entire hurricane boundary layer is usually combined with the strong wind shear near the eyewall to create intense rolls there. Then, the pressure transport term transfers the TKE horizontal component to organize the vertical component of rolls and eddies near the hurricane eyewall. After hurricane landfall, the weak horizontal advection and weak wind shear reduce the energy support and decrease the strong coherent eddies, including rolls and eddies.

Wavelength
In this section, the wavelength of strong eddies and that of pressure perturbation are compared to clarify the connection between eddies and pressure perturbation. Four cross sections from the hurricane center to the north, east, south, and west at a height of 1000 m are examined at each time to obtain pressure perturbation data. Only data from within a radius of 20-40 km are considered. With this cross-section data, pressure perturbation wavelength is detected by the distance of local maxima or minima. Figure 3c shows an example of the probability distribution of the wavelength of the pressure perturbation during hurricane landfall at 0100 UTC on 26 August 2017. Similar to the probability distribution of the wavelength of eddies, the probability distributions of the wavelengths of pressure perturbation from local maxima and minima are close to each other, with global maxima at around 0.5-1 km. The mean wavelength is 0.9 km. The similar mean wavelengths from eddies, rolls (Fig. 3a, b), and pressure perturbation (Fig. 3c) again confirm the relationship between pressure perturbation, rolls, and eddies. Once inflow convergence transports the rolls and creates strong rolls and eddies with strong wind shear near the eyewall, the strong pressure perturbation can then transfer the eddy horizontal component (Iu) into the vertical component (Iw). Even after landfall, when strong rolls and eddies almost disappear as the wind shear weakens, the pressure perturbation can still transfer weak horizontal turbulence (Iu) into weak vertical turbulence (Iw), with a mean wavelength of about 0.9 km.

Summary
In this study, the one-way nested WRF-LES model simulation results from  are used to examine the structure and generation mechanism of rolls and associated coherent eddies in the hurricane boundary layer around Hurricane Harvey's landfall from 00 UTC 25 to 18 UTC 27 August 2017. The simulation results indicate that the aligned downstream rolls prevail throughout the entire hurricane boundary layer from near the surface to a height of 200-400 m. Extremely intense rolls with associated intense eddies occur at a height of 200-3000 m, within a radius of 20 to 40 km from the storm center.
During the simulation period, the vertical and horizontal components of the turbulence intensity of the intense coherent eddies, including rolls and other eddies, are relates to the hurricane intensity. In other words, turbulence intensity increases with hurricane intensification. Before and during hurricane landfall, strong inflow convergence leads to horizontal advection of rolls throughout the entire hurricane boundary layer and combines with the strong wind shear there to generate the strong rolls and eddies near the eyewall with suitable thermodynamic (Rf at around − 0.2 to 0.2) and dynamic conditions (strong negative inflow wind shear). After landfall, the decayed inflow weakens the inflow convergence and wind shear, leading to the quick reduction of strong rolls.
Once strong rolls and eddies occur near the eyewall, atmospheric pressure perturbation, which is caused by horizontal convergence, forces some TKE horizontal components (Iu) to turn into vertical components (Iw). The pressure perturbation force connects the rolls and eddies vertically with a similar wavelength of about 1 km. Meanwhile, the buoyancy term in the TKE equation is weak and negative; thus, it suppresses the vertical component of rolls and eddies. Therefore, under the necessary thermodynamic and dynamic conditions, the hurricane flow will generate strong rolls near the eyewall and then combine with air pressure perturbation to generate strong vertical turbulence motion near the eyewall of the hurricane. The simulated wavelength structure is similar to previous observations. More case studies will be conducted in future work. and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Publisher's Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.