• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Modulation of the Wind Field Structure of Initial Vortex on the Relationship between Tropical Cyclone Size and Intensity

    2023-09-07 07:49:30KexinCHENGuanghuaCHENandDongleiSHI
    Advances in Atmospheric Sciences 2023年9期

    Kexin CHEN, Guanghua CHEN, and Donglei SHI

    1Key Laboratory of Cloud-Precipitation Physics and Severe Storms, Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China

    2Shanghai Typhoon Institute, China Meteorological Administration, Shanghai 200030, China

    3University of Chinese Academy of Sciences, Beijing 100049, China

    4Department of Atmospheric Science, School of Environmental Studies,China University of Geosciences, Wuhan 430074, China

    ABSTRACT This study investigates the modulation of initial wind field structure on the relationship between the size and intensity of a simulated vortex.A series of idealized experiments are conducted by varying the radius of maximum wind (RMW) and shape parameter b of the initial vortices.The size—intensity relationship is quantified by the linear regression coefficient of the azimuthally-averaged gale-force wind radius against the maximum wind during the development stage, reflecting the degree of size expansion at the same intensity increment.The regression coefficient increases with increased RMW and decreased b, with the RMW being the primary constraint.The effect of lowering b on the elevation of the regression coefficient gradually stands out when the initial RMW is large.Enlarging the RMW leads to a secondary circulation with a horizontally elongated structure, which retards the intensification while promoting size expansion, thus substantially enhancing size expansion as the vortex intensifies.Broadening the wind field outside the RMW by reducing b results in abounding convection in the outer region, which promotes size expansion.Based on the axisymmetric tangential wind tendency and Sawyer—Eliassen equations, when the RMW is large, the active convection in the outer region can weaken the radial inflow induced by the eyewall heating in the inner region, thus retarding the intensification by reducing the radial imports of vorticity near the RMW.

    Key words: tropical cyclone, size and intensity, initial wind field structure, tangential wind tendency equation Sawyer-Eliassen equation

    1.Introduction

    Observational studies of tropical cyclones (TCs), from the large-sample perspective that mixes up all TCs’records,have shown a generally weak correlation between intensity and size (Merrill, 1984; Weatherford and Gray, 1988;Chavas and Emanuel, 2010; Chan and Chan, 2012; Musgrave et al., 2012; Knaff et al., 2014; Wu et al., 2015; Chavas et al., 2016; Guo and Tan, 2017; Song et al., 2020; Sun et al.,2022).However, the large-sample framework has inherent limitations, which are neglected in previous literature.On one hand, the overall weak correlation obscures the intrinsic size—intensity correlation of a single TC.For example, size is observed in reality as well as in numerical simulation to correlate well with intensity during the early development stage in terms of individual TCs (Xu and Wang, 2010; Chan and Chan, 2014; Chavas et al., 2016; Schenkel et al., 2018;Chen et al., 2019; Martinez et al., 2020; Song et al., 2020).On the other hand, the large-sample framework fails to identify individual differences among TCs, thus neglecting the fundamental cause of the overall weak correlation.Recently,we pointed out that although size is correlated well with intensity during the early development stage of individual TCs,the size change at the same intensity increment is largely discrepant among TCs, therefore leading to a weak correlation if mixing up all records (Chen et al., 2022).

    To provide a comprehensive picture of TC size—intensity relevance, besides the overall correlation, the variation of size—intensity relationship from storm to storm, i.e., the degree of size expansion relative to the same intensity increment, should be investigated.A variety of factors in the environment and internal processes can affect the evolution of TC intensity and size.The environmental factors include latitude (Chan and Chan, 2014; Smith et al., 2015), deep-layer vertical wind shear (Frank and Ritchie, 2001; Wang and Toumi, 2019), underlying ocean temperature (Xu and Wang,2018a; Chen et al., 2021), and synoptic systems at higher or lower levels (Hanley et al., 2001; Liu and Chan, 2002).The internal dynamics consist of the potential vorticity mixing between eye and eyewall (Tsujino and Kuo, 2020; Shi and Chen, 2021), eyewall replacement cycle (Zhou and Wang,2011; Huang et al., 2012), spiral rainband activity (Wang,2009), and frictionally induced convergence and dissipation within the boundary layer (Li and Wang, 2021a).By conducting idealized numerical experiments, Wang and Lin (2020)found that a dry TC tends to be of a larger size than a moist one and further pointed out that the hydrometeors’falling within a TC circulation could limit the size and structure of TCs in nature.Additionally, based on the balanced vortex dynamics (e.g., Shapiro and Willoughby, 1982), the wind field structure is an important internal factor, since it determines how a TC responds to a given environmental or internal forcing.

    Numerous analyses have shown that the evolution of intensity and size depends on TC wind field structure.The radius of maximum wind (RMW) and the radial decay rate of tangential wind outside the RMW are crucial structural metrics.The RMW reflects the inner-core extent of a TC, and the radial decay of tangential wind determines the outer size of the wind field.The subsequent intensification rate of a TC possessing a smaller RMW or outer size (measured by the gale-force wind radius, R34) is usually higher, such that a TC is more likely to experience rapid intensification(Chen et al., 2011; Carrasco et al., 2014; Xu and Wang,2015, 2018a).In addition, TC intensity significantly correlates with the structure fullness measured by the ratio of the width of the outer annular wind ring (R34-RMW) to the R34 (Guo and Tan, 2017; Chen and Li, 2021).Guo and Tan(2022) further uncovered the existence of a preferable fullness for a specific intensity that raises the great intensification rate in nature.In terms of the impacts of the initial wind field on the size evolution, the initially large TC generally maintains its size superiority over the smaller counterpart throughout its lifetime (Lee et al., 2010), as it possesses higher absolute angular momentum import in the outer region (Chan and Chan, 2014).The higher momentum import arises mainly from a larger amount and larger-scale convection in the broader wind field (Martinez et al., 2020).

    The fundamental pathway leading to the spin-up of tangential wind in the inner and outer regions is driven by the convectively induced inflow at lower levels to advect angular momentum inwards (Schubert and Hack, 1982; Shapiro and Willoughby, 1982).Smith et al.(2009) pointed out that TC intensification is associated with the contraction of the absolute angular momentum surface within the boundary layer near the RMW, while the size expansion is related to the momentum convergence above the boundary layer in the outer region.Observations and numerical experiments show that the kinematic feature of TC boundary radial inflow varies with the wind field structure.Specifically, the TC with a narrow outer size or a small RMW is more likely to possess a deeper inflow layer with stronger updrafts close to the RMW (Kilroy and Smith, 2017; Zhang et al., 2023).

    The above review demonstrates that more attention has been paid to the influence of wind field structure on intensity and size.However, a quantitative analysis on the impact of vortex initial structure on the size—intensity relationship, i.e.,the degree of size expansion relative to the same intensity increment, is still lacking.Specifically, the modulation of the two structural parameters (RMW and decay rate of tangential wind outside the RMW) and their relative importance on the size—intensity relationship remains unknown, especially during early development.Following our previous work (Chen et al., 2022), the size—intensity relationship is quantified by the linear regression coefficient of size against intensity given that the size generally correlates well with the intensity during early development.This study aims to investigate how the initial wind field structure modulates the size—intensity relationship during early development and further evaluate the relative importance of initial RMW and radial decay rate of tangential wind outside the RMW.

    The remainder of this paper is organized as follows.Section 2 describes the numerical model settings and experimental design.Section 3 explains how the initial wind field structure modulates the size—intensity relationship.The thermodynamical process underlying the modulation of the RMW and radial decay of wind outside the RMW on the size—intensity relationship is investigated in section 4.Conclusions and discussions are given in section 5.

    2.Model configurations and experimental design

    2.1.Model configurations

    Idealized simulations are carried out with the Advanced Research Weather Research and Forecasting Model Version 3.6.1.The model domain is triple nested with grid dimensions of 301×301, 211×211, and 301×301, spacing by 27 km,9 km, and 3 km, respectively.Except for the outermost domain, two inner meshes are set to automatically follow the vortex.There are 35 levels in the vertical direction with the model top at 50 hPa.Experiments are initialized on anfplane at 20°N with a uniform SST of 28°C.

    For all experiments, the physics parameterizations used are as follows.The Thompson scheme (Thompson et al.,2008) and Yonsei University (YSU) scheme (Hong et al.,2004) are chosen for model microphysics and boundary layer physics, respectively.The Kain—Fritsch cumulus parameterization scheme (Fritsch and Kain, 1993) is used only for the outermost domain.The Rapid Radiative Transfer Model(RRTM) longwave radiation scheme (Mlawer et al., 1997)and shortwave radiation scheme (Dudhia, 1989) are used.Note that the radiation is fixed at noon in all experiments to remove the influence of diurnal solar forcing.

    The tangential wind (vt) profile of the axisymmetric initial vortex following Xu and Wang (2018b) is given by:

    2.2.Experimental design

    A series of experiments are conducted by varying the RMW and the shape parameterbto investigate how the initial wind field structure affects the size—intensity relationship of a vortex, that is, the degree of size expansion relative to the same intensity increment.Three settings of RMW of 80 km,120 km, and 160 km, andbof 0.75, 0.50, and 0.35, yield nine experimental combinations.Specifically, the initial vortex in experiment R80b75 possesses the smallest RMW of 80 km and the fastest decay of tangential wind outside the RMW with abof 0.75, while experiment R160b35 has the largest RMW of 160 km and the slowest wind decay outside the RMW with abof 0.35.To further validate the conclusions, experiments with the RMW reduced to 40 km are also carried out.The modulation of the initial wind field on the size—intensity relationship is qualitatively consistent with other experiments, and the experimental results with the RMW of 40 km are not shown for brevity.Figure 1 shows the radial wind profiles and corresponding inertial stability of the initial vortex in the experiment groups R80 and R160.All experiments were integrated over five days, long enough for a simulated vortex to complete early development.

    3.Experiment results

    3.1.Evolution of the vortex

    Figure 2 presents the evolutions of intensity and size, as well as the RMW of the simulated vortex in all experiments.Here, the intensity is defined as the maximum azimuthally averaged tangential wind speed (Vmax), and the size is measured by the azimuthally averaged R34; both are obtained at 10-m height.The hourly intensity change over a 6-hour period is converted to an equivalent value over 24 hours to determine the development stage.The development stage begins at the time when the intensification rate first exceeds 15 m s—1d—1and terminates at the first time when the intensification rate is less than 5 m s—1d—1for six consecutive hours,which is similar to Xu and Wang (2018b) and Martinez et al.(2020).The beginning and end times of the development stage in each simulation are marked by the red and black dots, respectively, shown in Fig.2.

    Fig.1.Radial profiles of (a) tangential wind and (b) inertial stability of the incipient vortex.For clarity, those only in the experiment groups R80 and R160 are displayed.

    Smaller RMW vortices have a shorter spin-up period followed by a faster intensification (Fig.2a), which is consistent with previous research (Kilroy and Smith, 2017; Xu and Wang, 2018b; Li and Wang, 2021a).Note that the smaller RMW contracts to smaller radii faster than the large RMW(Fig.2c); compared to the larger RMW vortices, the structure of the smaller RMW vortex changes into the preferable configuration for subsequent intensification faster.Explicitly,the fullness ratio of the vortex [measured by the ratio of TC fullness (T CF=1-RMW/R34) to the critical fullness(TCF0=1-17.0/Vmax)] reaches 1.4 (not shown, Guo and Tan, 2022).The RMW also serves as the primary constraint on size expansion.Although R34 appears later as RMW increases due to slow intensification, the size expansion rate increased with RMW, which is consistent with previous research (Kilroy and Smith, 2017).

    Fig.2.Temporal evolutions of (a) the maximum azimuthally averagedtangentialwindspeed(Vmax,ms—1),(b)the azimuthallyaveragedgale-forcewindradius (R34,km),and(c) the radius of maximum wind (RMW, km) for all experiments.The red and black dots denote the beginning and end times of the development stage of each experiment.

    The modulation of the shape parameterbon intensification is less pronounced as compared to the RMW and is furthermore dependent on it.Specifically, when the initial RMW is small (i.e., 80 km and 120 km), modifying thebbarely influences the intensification.However, when the initial RMW is large (i.e., 160 km), broadening the wind field in the outer region by decreasing thebslows down the intensification, especially during the second half of the development stage.The bifurcation point of intensity evolution in the R160 group in Fig.2a is approximately the time the discrepancy in size expansion occurs in Fig.2b.It is conceivable that when the inner inertial stability weakens with increasing RMW, the radial vorticity may be blocked from moving inward by the outer convection thus slowing intensification,which will be quantitively examined in section 4.Consistent with previous studies (Chan and Chan, 2014; Martinez et al.,2020), given the same RMW, the size expansion increases with the broadening outer wind field specified by decreasingb.

    3.2.Relationship between intensity and size

    In idealized simulations without external forcing, there is a positive linear correlation between size and intensity during the development stage with a correlation coefficient close to 1 for each vortex, as shown in Table 1.The linear regression coefficient of R34 againstVmax, which shows the degree of size expansion relative to the same intensity increment, can be another key metric to quantify the size—intensity relationship.The higher the regression coefficient is, the more pronounced the size expansion relative to the same intensity increment.Figure 3 displays the regression results for all experiments with the regression coefficients shown in Table 1.The varying size—intensity relationship among the vortices results in a wide range of size for a given inten-sity, which can account for the weakening of the overall size—intensity correlation if mixing up all experiment records.In particular, the correlation coefficient for all the experiment data in Fig.3 drops to 0.57 (Table 1), indicating that the large-sample framework can obscure the intrinsic size—intensity correlation in individual vortices.

    Table 1.Linear correlation and regression coefficients between size and intensity in each experiment and all experiment data during the development stage.

    Fig.3.Linear regressions of size (R34) against intensity (V max)during the development stage of all experiments.Specific linear correlation and regression coefficients between size and intensity for each experiment are given in Table 1.

    Different slopes of the regression lines demonstrate that the initial wind field structure can modulate the size—intensity relationship.Overall, either increasing the RMW or decreasing the shape parameterbelevates the regression coefficient(Table 1).In comparison, increasing the initial RMW elevates the regression coefficient more evidently than reducingb.The more pronounced effect of enlarging the RMW is due to its slowing intensification while accelerating expansion(Fig.2).The elevation by reducingbis of small magnitude and depends on the RMW.The larger the RMW is, the steeper the regression lines can become with decreasingb,which is associated with the more noticeable inhibition to intensification as RMW increases.In short, the initial RMW plays the leading role in modulating the evolution of intensity and size, and it is thus the primary factor regulating the size—intensity relationship.Therefore, compared to the wind structure in the outer region, the early RMW may warrant more attention when forecasting R34 changes.

    4.Quantitative evaluation of the sizeintensity relationship

    4.1.Tangential wind budget analysis

    Based on the above discussions, two questions are naturally raised: How does the initial RMW and shape parameterbaffect the evolution of intensity and size, and thus the size—intensity relationship? What is the quantitative contribution of RMW andbto the modulation of the size—intensity relationship? The evolution of intensity and size can be investigated by the tangential wind tendency equation.From an axisymmetric perspective, the tangential wind tendency equation can be written as

    wherevtandvrare the azimuthal-mean tangential and radial wind, respectively, η=ζ+f0is the azimuthal-mean absolute vertical vorticity,ζis the azimuthal-mean relative vertical vorticity,f0is the Coriolis parameter at the vortex center, andFis the surface friction from the model output.The three terms on the right-hand side of Eq.(2) are the radial advection of the azimuthal-mean absolute vertical vorticity, vertical advection of the azimuthal-mean tangential wind, and the surface friction, respectively.Note that contributions by eddy terms are quite small and negligible and thus were not included in Eq.(1).The modulation of the initial wind field structure on the size—intensity relationship can be understood and evaluated by comparing the three terms among different experiments.The typical experiments R80b75, R80b35, R160b75, and R160b35 are selected for further comparison.Moreover, the investigation focuses on the primary development stage of the vortex intensifying from 15 m s—1to 35 m s—1for a fair comparison.

    Figure 4 shows the radial profile ofvtaveraged vertically from 0 km to 1 km for every 5-m s—1increment in intensity to present how the size—intensity relationship results in wind profile changes.In the R80 group,vtamplifies mainly near the RMW, while outside the RMW amplification is not evident; thus, the slope ofvtoutside the RMW is steep with time (Figs.4a, b).The averaged tangential wind at lower levels of 17.5 m s—1(34 kt) expands from 40 km to 80 km as the vortices gradually intensify.In comparison, the expansion of the wind profile outside the RMW is more pronounced in the R160 group.Moreover, decreasing thebfrom R160b75 to R160b35 is more conducive to size expansion compared to going from R80b75 to R80b35 (Figs.4c, d).In the R160 group, the radius range of 120 km to 180 km is the main extension rage of the lower wind of 17.5 m s—1.For brevity,the radial ranges from 40 km to 80 km and 120 km to 180 km will henceforth be referred to as the outer region.

    The temporally averagedvtbudgets vertically from 0 km to 3 km based on the model outputs with 10-min time intervals are shown in Fig.5.Consistent with previous studies(i.e., Xu and Wang, 2018b), the radial advection of absolute vertical vorticity imposes the primary positive contribution below 1 km in height (second column in Fig.5) while the friction is the opposite (fourth column in Fig.5), both with the maximum residing near the RMW in all experiments.The vertical advection of the azimuthal-mean tangential wind mainly spins down thevtaround the RMW below 1 km in height (third column in Fig.5).

    Fig.4.The radial profile of azimuthal-mean tangential wind averaged vertically from 0 km to 1 km for every 5-m s—1 increment in intensity when the simulated vortex intensifies from 15 m s—1 to 35 m s—1 in (a) R80b75, (b) R80b35, (c) R160b75, and (d) R160b35.The blue, green, yellow, orange, and red lines denote the intensity of the vortex of 15 m s—1, 20 m s—1, 25 m s—1, 30 m s—1, and 35 m s—1,respectively.The horizontal black line in each panel indicates the 17.5-m s—1 (34-kt) wind speed.

    Compared to the R160 group, the radial vorticity advection near the RMW is higher in the R80 group (Figs.5b and 5f), which results in a greater positivevttendency in the inner region (Figs.5a and 5e) and a corresponding higher intensification rate.When the RMW is small, reducing thebincreases the positivevttendency slightly in the radial range of 40 km to 80 km going from R80b75 to R80b35, justifying the relatively higher expansion rate of the latter.In contrast,the positivevttendency of R160b35 radially within 60 km and vertically below 1 km noticeably weakened compared to R160b75 (Figs.6i and 6m), matching its slower intensification during the second half of the development stage.Note that the maximum of the positivevttendency and radial vorticity advection of R160b35 are relatively higher than those of R160b75, which is associated with the evaporation-cooling effects of convection outside the eyewall and will be illustrated in more detail later.Akin to the R80 group, contributed by the higher radial advection of vorticity in the outer region (Figs.6j and 6n), reducing thebfrom R160b75 to R160b35 increases the positivevttendency outside 120 km.

    The radial advection of absolute vertical vorticity dominates thevtchanges at lower levels.To examine the key process leading to the varying size—intensity relationship between different experiments, Fig.6 plots the time evolution of the radial vorticity advection, as well as the radial inflow averaged vertically from 0 km to 1 km during the period when the vortex intensifies from 15 m s—1to 35 m s—1in each experiment.As expected, the radial advection of vorticity evolves as a function of the radial inflow.The largest radial vorticity advection resides along the RMW, where the radial gradient of the inflow is the strongest.The radial gradient of the inflow is greater near the RMW with a much narrower extension in the R80 group than in the R160 group,which demonstrates the horizontally restricted secondary circulation driven by the eyewall heating with higher inertial stability of the former (Fig.1b) (Schubert and Hack, 1982;Tsuji et al., 2016).After the 47th and 67th hour, the radial inflow of R80b35 and R160b35 weakens in the inner region and extends to larger radii compared to R80b75 and R160b75, respectively.Moreover, the weakening is more noticeable in R160b35.

    Fig.5.Radius—height cross sections of the tangential wind tendency budget.The sum of all terms contributing to the tangential wind tendency (first column), radial advection of azimuthal-mean absolute vertical vorticity (second column),vertical advection of azimuthal-mean tangential wind (third column), and the friction term (fourth column) are averaged during the period when the vortex intensifies from 15 m s—1 to 35 m s—1 in R80b75 (first row), R80b35 (second row),R160b75 (third row), and R160b35 (fourth row).The grey line denotes the radius of 40/80 km in the R80 group and 120 km in the R160 group, respectively.

    Changes in radial inflow are closely associated with convection.Figure 7 compares the snapshots of convection represented by the simulated radar reflectivity after the 47th and 67th hour of the R80 and R160 groups.Similar to previous studies (Kilroy and Smith, 2017; Martinez et al.,2020), the convection gradually propagates towards the vortex center with a narrow wind field outside the RMW in R80b75 and R160b75, while the outer convection is more active outside the radius of 40 km and 120 km in R80b35 and R160b35, respectively.Previous studies (Kilroy and Smith, 2017; Zhang et al., 2023) have claimed that the broad wind field outside the RMW tends to generate upflow above the boundary layer to burst deep convection in the outer region.Additionally, the corresponding relatively higher inertial stability in the outer region can resist the radial inflow induced by eyewall heating such that the outer convection less easily gathers in the eyewall.The positive feedback between deep convection and radial inflow in the outer region favors the size expansion, while this process could restrict intensification by blocking mass and moisture supply toward the inner region (Powell, 1990; Wang, 2009).This justifies the weakening of radial inflow in the inner region when the radial extension broadens in R80b35 and R160b35.

    Fig.6.Temporal evolutions of the radial advection of the absolute vertical vorticity (m s—2,shading) and radial inflow (m s—1, contoured in black line every 3 m s—1) averaged vertically from 0 km to 1 km during the period when the vortex intensifies from 15 m s—1 to 35 m s—1 in(a) R80b75, (b) R80b35, (c) R160b75, and (d) R160b35.The radius of maximum wind is plotted as the green line.The white line denotes the radius of 40/80 km in (a, b) and 120 km in (c, d), respectively.

    In summary, the initial wind field specified by the RMW andbcan regulate the distribution of inertial stability and convection in the eyewall and outer region, thus modulating the extension and strength of the convectively induced radial inflow at lower levels.Therefore, the radial advection of absolute vertical vorticity and the changes in tangential wind differ with different wind field structures, leading to a varying size—intensity relationship.By solving the SE equation, we calculate in the following subsection the balanced response of radial inflow to heating released in the eyewall and outer convection to quantitatively evaluate the relative modulation of the RMW andbon the size—intensity relationship.

    4.2.Diagnostics of the SE equation

    Fig.7.Snapshots of simulated radar reflectivity (dBZ, shading) at 1-km height at different times for (a—d) R80b75, (e-h)R80b35, (i-l) R160b75, and (m-p) R160b35.The circles are placed at every 40-km radius from the vortex center.

    To quantitatively assess the contribution of the heating released in the eyewall and outer convection to thevttendency at lower levels, the SE equation following Bui et al.(2009) is adopted as follows:

    Since the SE equation is formulated based on the assumption of gradient wind balance (Bui et al., 2009), it is necessary to examine to what degree this assumption is satisfied.The agradient winds are calculated asvtminusvg, wherevgis the wind satisfying the gradient wind relation in the azimuthal average.Therefore, positive and negative agradient winds indicate the supergradient and subgradient winds,respectively.In general, a large portion of tangential winds at lower levels are in gradient wind balance, except for the subgradient wind in the boundary layer and the supergradient wind near the maximum of tangential winds under the eyewall(not shown), which is consistent with previous studies(Fudeyasu and Wang, 2011; Ge et al., 2016).Although the SE equation cannot fully explain the imbalance process in the boundary layer due to the assumption of gradient wind balance, it has been proven in previous studies to be a valuable tool to evaluate the individual balanced contribution of different forcings (Fudeyasu and Wang, 2011; Chen, 2018).Li and Wang (2021b) reported that the momentum forcing within the boundary layer is sensitive to the initial wind field structure and thus could influence the evolution of intensity and size.This study will focus on the contribution of heating forcing.

    Fig.8.Radius—height cross sections of the azimuthally averaged radial (m s—1, contour) and vertical velocities (m s—1,shading) based on the model output (first column), Sawyer—Eliassen equation (second column), and the balanced component induced by eyewall heating (third column) and outer convective heating (fourth column).All the results are averaged during the period when the vortex intensifies from 15 m s—1 to 35 m s—1 in (a—d) R80b75, (e—h) R80b35, (i—l) R160b75, and (m—p)R160b35.

    Figure 8 presents the azimuthal-mean vertical and radial velocities averaged during the 15—35-m s—1period from model output and diagnosed by the SE equation.In general, the SE equation can capture a large proportion of the circulation, except that it underestimates the inflow in the boundary layer and outflow in the upper level.Additionally, the vertical extension of the diagnosed vertical velocity is lower,and the strength is weaker than the model output.These differences are due to the unbalanced process and upper-level eddy effect (Bui et al., 2009; Fudeyasu and Wang, 2011).It is worth noting that the structure of vertical velocity greater than 0.5 m s—1is thinner in R160b35 than in R160b75,which might be due to the evaporation cooling effect of outer convection on the eyewall.As shown in Fig.7, the active convection in the outer region of R160b35 tends to restrict the inner rainbands compared to that in R160b35,reflected by the thinner reflectively high value in the eyewall convection.As purported by Li et al.(2015), the enhanced evaporation cooling from convection outside the eyewall can reduce the inner-core extension, such that thinner vertical velocity in the eyewall convection can be expected.

    We further divide heating into two regions by taking the azimuthal-mean diabatic heating rate within two times the RMW as the eyewall convection, and that which covers the area from two to four times the RMW is taken as the outer convection.Such a partition is determined by investigating the diabatic heating rate at each moment.Note that the RMW is calculated at each level since the eyewall convection tilts outward with height.The balanced response of the wind field to the eyewall heating and outer convective heating is also represented in Fig.8.The responses of the wind field to eyewall heating are the primary vertical motion within the eyewall and the radial inflow maximizing outside the eyewall at lower levels (third column in Fig.6).The heating released in outer convection drives a weak overturning circulation in the outer region (fourth column in Fig.6).Compared to the R80 group, the secondary and the overturning circulations in the R160 group have a weaker vertical motion at larger radii with a horizontally broad extension of the radial inflow at lower levels.The weak vertical velocity of the overturning circulation in the R160 group is mainly due to the loss of heating to gravity wave radiation with lower inertial stability at large radii.Additionally, the lower inertial stability implies weaker resistance to the radial movement, which facilitates the radial inflow of the overturning circulation to extend outside R120 in the R160 group (Figs.6l and p).

    Figure 9 compares the contribution of radial inflow induced by heating in the eyewall and outer convection to thevttendency.The temporal evolution of the primary contribution term, the radial advection of vorticity, is shown in Fig.10.The eyewall heating is the primary source for the positivevttendency at lower levels, except within the small radii from the vortex center (first column in Fig.9).Contributions of outer convective heating, albeit small in magnitude,acceleratevtoutside the heating forcing while deceleratingvtradially inward from the heating position (second column in Fig.9), which supports the claim that the outer convection favors size expansion while retarding intensification.Comparison within the R80 group implies that decreasing thebh ardly affects thevttendency supported by eyewall heating in the inner region (Figs.9a and 9c), which is also reflected in the quite similar evolution of radial inflow near the RMW (Figs.10a and 10b).However, there is a noticeable decrease in thevttendency and radial inflow induced by eyewall heating from R160b75 to R160b35 (Figs.9e and 9g,10c and 10d).When the RMW is large, decreasing thebpromotes the convection in the outer region, leading to the extension of the radial vorticity import to larger radii (Figs.10g and 10h).In this case, the radial inflow induced by the eyewall heating with lower inertial stability seems to be more vulnerable to the influences from active outer convection in a broader wind field.

    In summary, the secondary circulation induced by eyewall heating is the primary source for spinning up thevtat lower levels, which is modulated by the RMW since it influences the eyewall—convection organization and energy conversion efficiency.The influences of outer convective heating,although small in magnitude, spin up the tangential wind in the outer region to facilitate size expansion while spinning down the wind inside the heating position.When the RMW is large, the eyewall-induced radial inflow, thus the positivevttendency, in the inner region can be weakened by the outer convection.The stand-out negative influences of the outer convection to intensification as RMW increases explain why the modulation of varyingbstands out in the vortices with an initially large RMW.

    5.Conclusion and discussion

    This study investigates the modulation of initial wind field structure on the relationship between the size and intensity of a vortex during the development stage with a series of idealized simulations.The intensity and size are defined as the azimuthally averaged maximum wind (Vmax) and galeforce wind radius (R34), respectively.The initial structure is modified by the radius of maximum wind (RMW; 80 km,120 km, and 160 km) and the shape parameterbdetermining the radial decay of tangential wind outside the RMW (0.35,0.50, and 0.75).

    The vortex size is found to be highly correlated with the intensity during the development stage in the individual experiments without external forcing.The size—intensity relationship is quantified by the linear regression coefficient of R34 againstVmax, measuring the degree of size expansion at the same intensity increment.The RMW is the primary factor modulating the size—intensity relationship.The degree of size expansion relative to the same intensity increment increases profoundly with RMW.Comparatively, the effect by reducingbis of smaller magnitude and gradually stands out as the RMW becomes large.Various size—intensity relationships among the vortices with different initial structures divert the size distribution relative to the same intensity,reducing the overall size—intensity correlation if mixing up all experiments’data.

    Fig.9.Radius—height cross sections of the tangential wind tendency contributed by the eyewall heating (first line), outer convective heating (second line), and the momentum forcing (third line) based on the Sawyer—Eliassen equation.Results are averaged during the period when the vortex intensifies from 15 m s—1 to 35 m s—1 in R80b75 (first line), R80b35 (second line),R160b75 (third line), and R160b35 (fourth line).The grey line denotes the radius of 40 km and 80 km in the R80 group and 120 km in the R160 group, respectively.

    The size—intensity relationship is determined by changes in the tangential wind near the RMW and R34,which mainly arise from the advection of absolute vertical vorticity by the radial inflow at lower levels.The radial inflow induced by convection in the eyewall and outer region is modulated by the RMW andb, respectively.Based on the axisymmetric tangential wind tendency and Sawyer—Eliassen equations, the modulation of RMW andbon the size—intensity relationship can be quantitatively evaluated by calculating the balanced response of the wind field to heating released in the eyewall and outer convection.

    Results show that heating released in the eyewall is the primary source of radial inflow.Enlarging the RMW reduces the inertial stability within the RMW substantially,which retards the eyewall-convection enhancement and leads to a secondary circulation with a horizontally elongated structure, thus hindering the intensification while promoting size expansion.Broadening the wind field outside the RMW by lowering thebfacilitates the convection in the outer region.The influence arising frombdecrease, although minor relative to RMW, spins up the tangential wind in the outer region to facilitate size expansion while spinning down the tangential wind radially inwards of the heating location.When the RMW is large, the broad convection in the outer region weakens the radial inflow induced by the eyewall heating, thus playing an inhibiting role in TC intensification.This explains why the modulation of varyingbstands out in the vortices with an initially large RMW.

    As an extension of previous studies that solely focus on the impact of wind field structure on the intensity or size,we discuss the modulation of wind field structure on the size—intensity relationship, which might provide, to a certain degree, a more comprehensive picture of the evolution of storm structure.Results in this study suggest that considering the wind radius warning, caution should be paid to the incipient vortex with a large RMW since the wind field can considerably expand to large areas as the vortex intensifies.This study is a preliminary attempt to understand the physical mechanisms behind the varying size—intensity relationship,and more in-depth work is ongoing to incorporate more complicated internal and external forcings.For example, the unbalanced process at the convective scale and the interactions between the vortex and external forcings (i.e., the vertical wind shear) deserve more attention.To what degree and how the unbalanced process and external forcings can influence the size—intensity relationship during the development should be addressed in a future study.

    Acknowledgements.This study is supported by the National Natural Science Foundation of China (Grant Nos.42175073 and 41975071).

    亚洲精品乱久久久久久| 日韩一区二区三区影片| 国产精品麻豆人妻色哟哟久久| 免费看光身美女| 日本黄大片高清| 日日撸夜夜添| 午夜91福利影院| 飞空精品影院首页| 97超碰精品成人国产| 99视频精品全部免费 在线| 看非洲黑人一级黄片| 一级二级三级毛片免费看| 欧美精品亚洲一区二区| 精品视频人人做人人爽| 久久久亚洲精品成人影院| 七月丁香在线播放| 一本大道久久a久久精品| 高清毛片免费看| 国产一区二区三区综合在线观看 | xxx大片免费视频| 午夜老司机福利剧场| 久久久精品区二区三区| 亚洲在久久综合| 国产永久视频网站| 亚洲内射少妇av| 大香蕉97超碰在线| 精品一区二区三卡| 日本黄大片高清| 精品久久久久久电影网| 国产av国产精品国产| 成年人午夜在线观看视频| 日日啪夜夜爽| 国产精品.久久久| 汤姆久久久久久久影院中文字幕| 久久久精品94久久精品| 人妻少妇偷人精品九色| 免费人成在线观看视频色| 你懂的网址亚洲精品在线观看| 欧美精品一区二区免费开放| 国产精品一区二区在线不卡| 久久久久精品性色| 国产精品麻豆人妻色哟哟久久| 久久鲁丝午夜福利片| 蜜臀久久99精品久久宅男| 国产成人精品婷婷| 少妇猛男粗大的猛烈进出视频| 天堂俺去俺来也www色官网| 久久久国产欧美日韩av| 伦理电影大哥的女人| av播播在线观看一区| 日韩av免费高清视频| 久久国产亚洲av麻豆专区| 国产日韩欧美视频二区| 日韩在线高清观看一区二区三区| 免费高清在线观看视频在线观看| 国产亚洲午夜精品一区二区久久| 黄色视频在线播放观看不卡| 欧美精品高潮呻吟av久久| 最近中文字幕2019免费版| 国产一区二区在线观看日韩| 91久久精品国产一区二区三区| 能在线免费看毛片的网站| 欧美 日韩 精品 国产| 久久久欧美国产精品| 在线看a的网站| 色婷婷av一区二区三区视频| 久久国内精品自在自线图片| 欧美精品高潮呻吟av久久| 成人亚洲欧美一区二区av| 哪个播放器可以免费观看大片| 青春草亚洲视频在线观看| av不卡在线播放| 老女人水多毛片| 老司机影院毛片| 婷婷色综合大香蕉| 国产成人a∨麻豆精品| 91精品国产九色| 欧美老熟妇乱子伦牲交| 成人手机av| av一本久久久久| 国产精品一区www在线观看| 亚洲国产毛片av蜜桃av| 男女免费视频国产| 精品久久国产蜜桃| 亚洲欧美一区二区三区黑人 | 80岁老熟妇乱子伦牲交| a级毛片在线看网站| 九色成人免费人妻av| 能在线免费看毛片的网站| 自线自在国产av| 欧美+日韩+精品| 亚洲激情五月婷婷啪啪| 久久午夜综合久久蜜桃| 午夜激情av网站| 青春草视频在线免费观看| 免费观看av网站的网址| 国产欧美另类精品又又久久亚洲欧美| 男人爽女人下面视频在线观看| 精品视频人人做人人爽| 如日韩欧美国产精品一区二区三区 | 大香蕉97超碰在线| 免费少妇av软件| 寂寞人妻少妇视频99o| 国产亚洲欧美精品永久| 妹子高潮喷水视频| 最近中文字幕2019免费版| 免费不卡的大黄色大毛片视频在线观看| 99久久精品国产国产毛片| 国产日韩欧美亚洲二区| 一级,二级,三级黄色视频| 男女国产视频网站| 99热这里只有精品一区| 我的老师免费观看完整版| 国产精品免费大片| 特大巨黑吊av在线直播| 三级国产精品欧美在线观看| 男女边吃奶边做爰视频| 国产精品久久久久久久久免| 老司机影院毛片| 亚洲国产欧美在线一区| 丝袜脚勾引网站| 午夜免费男女啪啪视频观看| 亚洲欧美成人综合另类久久久| 99九九在线精品视频| 亚洲av成人精品一区久久| 成人毛片a级毛片在线播放| a级毛片免费高清观看在线播放| 久久99热6这里只有精品| 十八禁网站网址无遮挡| 中国三级夫妇交换| 麻豆精品久久久久久蜜桃| 国产成人freesex在线| 草草在线视频免费看| 亚洲av免费高清在线观看| 久久久精品免费免费高清| 18在线观看网站| 国产黄片视频在线免费观看| 亚洲欧美色中文字幕在线| 久久久久久久久久久久大奶| 啦啦啦中文免费视频观看日本| 亚洲不卡免费看| 亚洲欧美一区二区三区黑人 | 啦啦啦在线观看免费高清www| 精品人妻熟女av久视频| 精品酒店卫生间| 99久久中文字幕三级久久日本| 国产片特级美女逼逼视频| 熟女av电影| 69精品国产乱码久久久| 午夜老司机福利剧场| 久久婷婷青草| 亚洲色图 男人天堂 中文字幕 | 亚洲欧美一区二区三区国产| 伦精品一区二区三区| 搡女人真爽免费视频火全软件| 婷婷色综合www| 少妇被粗大的猛进出69影院 | 亚洲av男天堂| 99热6这里只有精品| 欧美激情极品国产一区二区三区 | 一区二区三区乱码不卡18| 精品视频人人做人人爽| 成人国产麻豆网| 热99国产精品久久久久久7| 国产精品人妻久久久久久| 久久精品国产自在天天线| 国产精品久久久久成人av| 18+在线观看网站| 看非洲黑人一级黄片| 国产精品久久久久成人av| 18+在线观看网站| 看非洲黑人一级黄片| 中文字幕精品免费在线观看视频 | 亚洲精品亚洲一区二区| 亚洲,一卡二卡三卡| 又黄又爽又刺激的免费视频.| 日韩人妻高清精品专区| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说 | 国产亚洲一区二区精品| 五月开心婷婷网| 各种免费的搞黄视频| 狂野欧美白嫩少妇大欣赏| 亚洲精品美女久久av网站| 国产欧美另类精品又又久久亚洲欧美| 美女视频免费永久观看网站| 中文天堂在线官网| 丰满少妇做爰视频| 亚洲国产精品999| 美女国产视频在线观看| 丰满饥渴人妻一区二区三| 国产白丝娇喘喷水9色精品| 啦啦啦在线观看免费高清www| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 美女xxoo啪啪120秒动态图| 丝瓜视频免费看黄片| 国产精品 国内视频| 91在线精品国自产拍蜜月| 国产一区有黄有色的免费视频| 亚洲精品日韩在线中文字幕| 18禁在线无遮挡免费观看视频| 秋霞伦理黄片| 日本vs欧美在线观看视频| 欧美日韩av久久| 黄色一级大片看看| 成人漫画全彩无遮挡| 丁香六月天网| 久久精品久久久久久噜噜老黄| 国产老妇伦熟女老妇高清| 精品少妇黑人巨大在线播放| 激情五月婷婷亚洲| 精品久久久噜噜| 国产 一区精品| 春色校园在线视频观看| 久久久久久久久久久免费av| 丰满迷人的少妇在线观看| 涩涩av久久男人的天堂| 观看av在线不卡| 国产老妇伦熟女老妇高清| 国产一区二区三区综合在线观看 | 亚洲无线观看免费| 精品久久蜜臀av无| 国产深夜福利视频在线观看| 午夜av观看不卡| 国产精品久久久久成人av| 亚洲在久久综合| 国产亚洲欧美精品永久| 涩涩av久久男人的天堂| 亚洲精品日韩av片在线观看| 免费黄色在线免费观看| 亚洲欧美精品自产自拍| 制服人妻中文乱码| 精品视频人人做人人爽| 亚洲av.av天堂| 国产精品国产三级国产专区5o| 九九在线视频观看精品| 国产av国产精品国产| 久久久精品免费免费高清| 一本久久精品| 国产黄色视频一区二区在线观看| 校园人妻丝袜中文字幕| 欧美丝袜亚洲另类| 国产探花极品一区二区| 草草在线视频免费看| 91久久精品国产一区二区成人| 日韩欧美精品免费久久| 狂野欧美激情性bbbbbb| 飞空精品影院首页| 亚洲丝袜综合中文字幕| 成人毛片60女人毛片免费| 国产精品国产三级国产av玫瑰| 国产精品久久久久久av不卡| 嘟嘟电影网在线观看| 熟女av电影| 日本与韩国留学比较| 在线观看三级黄色| 一区二区三区乱码不卡18| 999精品在线视频| 欧美老熟妇乱子伦牲交| 国产极品天堂在线| a级毛片免费高清观看在线播放| 大香蕉97超碰在线| 亚洲av成人精品一区久久| 欧美日韩一区二区视频在线观看视频在线| av一本久久久久| 丝袜在线中文字幕| 老司机影院毛片| 成人漫画全彩无遮挡| 在线观看人妻少妇| 欧美精品一区二区免费开放| 午夜激情福利司机影院| 国产成人精品无人区| 一二三四中文在线观看免费高清| 人妻一区二区av| 国产国拍精品亚洲av在线观看| 另类精品久久| 日韩av不卡免费在线播放| 好男人视频免费观看在线| 国产毛片在线视频| 成人18禁高潮啪啪吃奶动态图 | 两个人的视频大全免费| 中文字幕精品免费在线观看视频 | 精品久久久久久久久亚洲| 51国产日韩欧美| 老司机影院毛片| 欧美性感艳星| 国产不卡av网站在线观看| 成人国语在线视频| 内地一区二区视频在线| 国产一区二区三区av在线| 久久精品夜色国产| av网站免费在线观看视频| 毛片一级片免费看久久久久| 欧美+日韩+精品| 免费久久久久久久精品成人欧美视频 | 亚洲av免费高清在线观看| videosex国产| 18禁在线播放成人免费| 另类亚洲欧美激情| 高清av免费在线| 亚洲图色成人| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| 久久免费观看电影| 国产精品一区www在线观看| 欧美三级亚洲精品| 午夜老司机福利剧场| 少妇高潮的动态图| 久久青草综合色| 亚洲欧美成人精品一区二区| 国国产精品蜜臀av免费| 夜夜爽夜夜爽视频| 欧美精品国产亚洲| 国产高清三级在线| 亚洲一级一片aⅴ在线观看| 亚洲人成网站在线播| 午夜av观看不卡| 三级国产精品片| 国产永久视频网站| 国产精品久久久久久精品电影小说| 肉色欧美久久久久久久蜜桃| 岛国毛片在线播放| 久久久午夜欧美精品| 高清av免费在线| 亚洲av综合色区一区| 国产成人freesex在线| 国产精品秋霞免费鲁丝片| 99久国产av精品国产电影| 80岁老熟妇乱子伦牲交| 只有这里有精品99| 国产综合精华液| 嘟嘟电影网在线观看| 九草在线视频观看| 18禁裸乳无遮挡动漫免费视频| av不卡在线播放| 男女无遮挡免费网站观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产欧美在线一区| 亚洲欧美一区二区三区国产| 中国美白少妇内射xxxbb| 国产一区二区三区av在线| 美女福利国产在线| 国产日韩欧美在线精品| 2022亚洲国产成人精品| 国产视频内射| av在线观看视频网站免费| 久久久久精品性色| 两个人的视频大全免费| 亚洲国产色片| 亚洲激情五月婷婷啪啪| 97精品久久久久久久久久精品| 成人毛片60女人毛片免费| 高清欧美精品videossex| 久久国产精品大桥未久av| 久久国产亚洲av麻豆专区| 十八禁网站网址无遮挡| 一二三四中文在线观看免费高清| 欧美激情极品国产一区二区三区 | 婷婷色av中文字幕| 亚洲天堂av无毛| 成人二区视频| 观看av在线不卡| 18禁动态无遮挡网站| 91国产中文字幕| 99热全是精品| av播播在线观看一区| 在线观看www视频免费| 18在线观看网站| 黄色视频在线播放观看不卡| 亚洲av中文av极速乱| 午夜激情av网站| 日韩人妻高清精品专区| 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 三上悠亚av全集在线观看| 最新中文字幕久久久久| 自线自在国产av| 日韩制服骚丝袜av| av在线app专区| 亚洲欧美成人综合另类久久久| 成年av动漫网址| 国产欧美日韩综合在线一区二区| 亚洲精品一区蜜桃| 日韩大片免费观看网站| 十八禁高潮呻吟视频| 精品久久久久久久久亚洲| 婷婷成人精品国产| 国产欧美日韩综合在线一区二区| 免费黄网站久久成人精品| 国产精品秋霞免费鲁丝片| 蜜桃久久精品国产亚洲av| 九九久久精品国产亚洲av麻豆| 纯流量卡能插随身wifi吗| 久久精品久久久久久噜噜老黄| 国产欧美日韩一区二区三区在线 | 丝袜喷水一区| 中文精品一卡2卡3卡4更新| av在线老鸭窝| 夫妻午夜视频| 18+在线观看网站| 黑人欧美特级aaaaaa片| 亚洲国产日韩一区二区| 精品少妇内射三级| 人体艺术视频欧美日本| 久久久久网色| 久久久久久久久久久久大奶| 在线观看www视频免费| 91aial.com中文字幕在线观看| 飞空精品影院首页| 国产伦精品一区二区三区视频9| 亚洲美女搞黄在线观看| 国产成人免费观看mmmm| 免费黄频网站在线观看国产| 国产亚洲欧美精品永久| 久久亚洲国产成人精品v| 国产精品女同一区二区软件| 日韩熟女老妇一区二区性免费视频| 久久 成人 亚洲| videosex国产| 美女cb高潮喷水在线观看| 永久网站在线| 亚洲欧美中文字幕日韩二区| 91国产中文字幕| 999精品在线视频| 亚洲欧洲国产日韩| 一个人免费看片子| 黑人猛操日本美女一级片| 免费看不卡的av| 欧美日韩视频高清一区二区三区二| 十八禁网站网址无遮挡| 亚洲国产最新在线播放| 亚洲av日韩在线播放| 免费黄色在线免费观看| 午夜av观看不卡| 亚洲美女搞黄在线观看| 爱豆传媒免费全集在线观看| 热99久久久久精品小说推荐| 人人妻人人添人人爽欧美一区卜| 老司机影院毛片| 我要看黄色一级片免费的| 久久久久久久久久久久大奶| 韩国高清视频一区二区三区| 一区二区av电影网| 午夜影院在线不卡| 中文字幕最新亚洲高清| 男人爽女人下面视频在线观看| 国产成人精品无人区| 熟女人妻精品中文字幕| 免费人妻精品一区二区三区视频| 丝袜在线中文字幕| 日日撸夜夜添| 久久精品熟女亚洲av麻豆精品| 久久ye,这里只有精品| 熟女av电影| 精品少妇久久久久久888优播| 国产精品不卡视频一区二区| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区三区四区免费观看| 18禁在线无遮挡免费观看视频| 免费高清在线观看日韩| 国产精品嫩草影院av在线观看| 伊人亚洲综合成人网| 久久精品久久久久久久性| 日本午夜av视频| 国产在线一区二区三区精| 欧美精品国产亚洲| kizo精华| 久久97久久精品| 熟女电影av网| 精品人妻熟女av久视频| 精品久久久久久久久亚洲| 十八禁高潮呻吟视频| 亚洲av福利一区| 久久av网站| 久久国产亚洲av麻豆专区| 天堂中文最新版在线下载| 最近的中文字幕免费完整| 汤姆久久久久久久影院中文字幕| 国产爽快片一区二区三区| 最近最新中文字幕免费大全7| 丁香六月天网| 美女国产视频在线观看| 久久精品国产亚洲av天美| 亚洲精品日韩av片在线观看| 少妇 在线观看| 22中文网久久字幕| 精品一区二区三区视频在线| 熟女电影av网| 国产精品久久久久久精品电影小说| 超色免费av| 丝瓜视频免费看黄片| 亚洲精品乱久久久久久| 蜜桃在线观看..| 卡戴珊不雅视频在线播放| 99精国产麻豆久久婷婷| 最近2019中文字幕mv第一页| 制服诱惑二区| 极品少妇高潮喷水抽搐| 精品一区在线观看国产| 久久鲁丝午夜福利片| 午夜老司机福利剧场| 亚洲精品视频女| 国产一级毛片在线| 99久久中文字幕三级久久日本| 久久毛片免费看一区二区三区| 精品少妇久久久久久888优播| 99国产综合亚洲精品| 水蜜桃什么品种好| 如日韩欧美国产精品一区二区三区 | 免费看光身美女| 国语对白做爰xxxⅹ性视频网站| 人成视频在线观看免费观看| 999精品在线视频| 亚洲精品色激情综合| 欧美变态另类bdsm刘玥| 国产成人a∨麻豆精品| 亚洲av日韩在线播放| 韩国高清视频一区二区三区| 99久久中文字幕三级久久日本| 99久久精品一区二区三区| 午夜免费鲁丝| 成人免费观看视频高清| 老司机亚洲免费影院| 亚洲精品日韩av片在线观看| 精品亚洲成国产av| 国产片内射在线| 免费av不卡在线播放| 精品人妻熟女av久视频| 亚洲av成人精品一二三区| 日韩视频在线欧美| 青青草视频在线视频观看| 老司机影院毛片| 欧美日韩av久久| 国产国语露脸激情在线看| 国产日韩欧美在线精品| 日韩成人av中文字幕在线观看| 精品国产一区二区三区久久久樱花| 欧美日韩av久久| 成人国产麻豆网| 青春草视频在线免费观看| 美女xxoo啪啪120秒动态图| 国产有黄有色有爽视频| 一边摸一边做爽爽视频免费| 我的女老师完整版在线观看| 少妇被粗大的猛进出69影院 | 欧美 亚洲 国产 日韩一| 成人毛片a级毛片在线播放| 狠狠精品人妻久久久久久综合| 一本—道久久a久久精品蜜桃钙片| 久久毛片免费看一区二区三区| 嫩草影院入口| 久久久亚洲精品成人影院| 亚洲图色成人| 伦理电影免费视频| 久久女婷五月综合色啪小说| 三级国产精品片| 亚洲av男天堂| 在线精品无人区一区二区三| 夜夜骑夜夜射夜夜干| 国产 一区精品| 在线观看一区二区三区激情| 亚洲国产精品一区三区| 永久网站在线| 婷婷成人精品国产| 精品久久久久久久久av| 亚洲国产日韩一区二区| 九色亚洲精品在线播放| 精品国产一区二区久久| 少妇被粗大的猛进出69影院 | 下体分泌物呈黄色| 亚洲精品日韩av片在线观看| 男人添女人高潮全过程视频| 毛片一级片免费看久久久久| 久久久亚洲精品成人影院| 极品少妇高潮喷水抽搐| 中文字幕最新亚洲高清| 老司机影院毛片| 国产69精品久久久久777片| 色吧在线观看| 免费看光身美女| 成年人午夜在线观看视频| 亚洲av成人精品一二三区| 人妻 亚洲 视频| 国产免费一区二区三区四区乱码| 国产精品国产三级专区第一集| 精品亚洲成国产av| 日本爱情动作片www.在线观看| 国产成人精品久久久久久| 纯流量卡能插随身wifi吗| 免费久久久久久久精品成人欧美视频 | 久久精品熟女亚洲av麻豆精品| 国精品久久久久久国模美| 一级,二级,三级黄色视频| 亚洲图色成人| 亚洲欧洲精品一区二区精品久久久 | 日产精品乱码卡一卡2卡三| 国产乱人偷精品视频| 18禁在线播放成人免费| 夜夜看夜夜爽夜夜摸| 99热这里只有是精品在线观看| 亚洲综合色网址| 国产成人精品福利久久| 少妇熟女欧美另类| 一本色道久久久久久精品综合| 不卡视频在线观看欧美| 人妻少妇偷人精品九色| 日韩制服骚丝袜av| 日韩成人伦理影院| 久久国内精品自在自线图片| 视频区图区小说| 91精品一卡2卡3卡4卡| 成年人午夜在线观看视频| 欧美bdsm另类| 国产伦理片在线播放av一区|