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

    The Roles of Barotropic Instability and the Beta Effect in the Eyewall Evolution of Tropical Cyclones※

    2022-12-07 10:26:54JieJIANGandYuqingWANG
    Advances in Atmospheric Sciences 2022年11期

    Jie JIANG and Yuqing WANG

    1Department of Atmospheric Sciences, University of Hawaii at Manoa Honolulu, Hawaii 96822, USA

    2International Pacific Research Center and Department of Atmospheric Sciences,University of Hawaii at Manoa Honolulu, Hawaii 96822, USA

    (Received 12 June 2021; revised 2 November 2021; accepted 25 November 2021)

    ABSTRACT Diabatic heating by convection in the eyewall often produces an annular region of high potential vorticity (PV) around the relatively low PV eye in a strong tropical cyclone (TC). Such a PV ring is barotropically unstable and can encourage the exponential growth of PV waves. In this study, such instability and the subsequent nonlinear evolution of three TC-like vortices having PV rings with different degrees of hollowness on an f-plane are first examined using an unforced, inviscid shallow-water-equation model. Results show that the simulated eyewalls evolve similarly to those in the nondivergent barotropic model. It is also found that the polygonal eyewall structure can be decomposed into vortex Rossby waves(VRWs) of different wavenumbers with different amplitudes, allowing for wave-wave interactions to produce complicated behaviors of mesovortices in the TC eyewall. The same set of PV rings has been examined on a beta-plane. Although the beta effect has been rendered unimportant to the eyewall evolution due to the relatively small scale of the inner-core circulation, this study shows that the beta effect may erode the coherent structure of mesovortices in the eyewall of an initially hollow PV-ring vortex. Mesovortices modeled on the beta-plane with a greater beta parameter tend to experience an earlier breakdown and enhanced radial gradients of the basic-state (azimuthally mean) angular velocity, followed by wave-wave, wave-flow interactions, leading to earlier merger and axisymmetrization processes. This result implies that the beta effect could be one of the forcings that shorten the lifetime of quasi-steady mesovortices in the eyewall of real TCs.

    Key words:tropical cyclones, vortex rossby waves, barotropic instability, Beta effect

    1.Introduction

    In a mature tropical cyclone (TC), strong convection and associated diabatic heating rates in the mid-tropospheric eyewall produce high cyclonic potential vorticity (PV) in the lower troposphere, which can be transported upward by the eyewall updraft. However, subsidence in the eye region suppresses convection, inhibits latent heat release, and upsets PV generation. As a result, a strong TC is often characterized by a high PV tower in the eyewall, with relatively low PV in the eye region. This off-centered PV ring structure is barotropically unstable; a phenomenon previously studied in a nondivergent barotropic model by Schubert et al.(1999). They showed that in an initially axisymmetric,barotropically unstable vortex, a small perturbation (deviation from the azimuthal mean) could grow exponentially with its growth rate depending on both the basic-state (azimuthally mean) vortex structure and the wavelength of the perturbation itself. Such oscillations propagate azimuthally around the vortex directed to the right of the radial basic-state PV gradients. Since the underlying dynamics of PV waves in a circular vortex are analogous to the planetary Rossby waves, they are termed vortex Rossby waves (VRWs; MacDonald,1968; Guinn and Schubert, 1993; Montgomery and Kallenbach, 1997; Wang, 2002). After the VRWs reach a finite amplitude, the nonlinear wave-wave and wave-flow interactions may result in the wave phase-locking, wave breaking,mesovortex merger, and axisymmetrization processes. In this paper, we extend the work of Schubert et al. (1999) to the divergent condition to examine the significance of divergence on both the linear and nonlinear development of VRWs.

    Kossin and Schubert (2001) adopted a thinner high-PV annulus than that used in Schubert et al. (1999) to better mimic intense or intensifying TCs in the real world, such as Hurricanes Gilbert (1988) and Guillermo (1997). They found that a thinner and larger PV annulus can quickly experience the development of an unstable mode and breakdown of the PV annulus and later roll up into mesovortices that may merge and end up with a steady vortex crystal, which gives the related streamlines a polygonal appearance. On the contrary, a smaller or radially broader PV annulus tends to relax into a monopole with the maximum PV residing at the vortex center, accompanied by a significant central pressure drop and a decrease of the maximum wind speed. They also suggested that such a PV redistribution process can occur without moist dynamics. Polygonal eyewalls and mesovortices have been observed in real TCs (Lewis and Hawkins, 1982; Muramatsu, 1986; Reasor et al., 2000;Kossin and Schubert, 2004).

    Other than unforced two-dimensional simulations where initial PV annuli are presumed, Hendricks et al.(2014) further explored TC eyewall evolution forced by continuous diabatic heating. Using a shallow-water-equation model where deep convective heating is parameterized as a mass sink, they showed that the linear evolution of vortices in a forced, divergent condition is qualitatively similar to the unforced, nondivergent case studied by Schubert et al.(1999). However, the TC intensity variability is also significantly affected by the radial structure, temporal evolution,and the radial location of diabatic heating.

    Most studies on the barotropic instability of TC eyewalls have focused on the internal dynamics without considering any environmental effects. One of the interactions between TCs and the environment is the so-called beta effect that results from PV conservation and the differential advection of planetary vorticity by the TC flow (Holland, 1983; Chan and Williams, 1987; Fiorino and Elsberry, 1989; Wang and Holland, 1996a, b). On the beta-plane, a pair of mesoscale,northeast-southwest-oriented, counter-rotating gyres (beta gyres), centered hundreds of kilometers away from the vortex center, would present themselves with southeasterly ventilation flow in between. These activities would have the net effect of translating the entire vortex northwestward in the Northern Hemisphere (the beta drift). Previous theoretical and numerical studies on the TC inner-core dynamics are performed on anf-plane with the beta effect ignored. Ignoring the beta effect on the eyewall evolution is justified by the relatively small radial scale of motion in the compact TC innercore region. However, some real TCs have large circulations with eye sizes ~100 km (such as Typhoon Carmen in 1960,Typhoon Winnie in 1997, etc.) at a scale by which the beta effect on the inner-core structure may indeed be relevant.Besides, careful consideration of any asymmetric forcing is desirable in the regime of nonlinear dynamics, under which even small perturbations can cause irreversible changes.This study intends to evaluate whether the beta effect exerts a significant impact on the development of TC eyewalls. To this end, barotropically unstable vortices are numerically simulated on both thef-plane and beta-plane to examine the beta effect on the eyewall evolution.

    The remainder of this paper is organized as follows. In section 2, we derive the two necessary conditions of barotropic instability in the context of the two-dimensional asymmetric balance (AB) framework. Description of experiments and initialization of the shallow-water-equation model are illustrated in section 3. Section 4 presents the eyewall evolution of vortices with different basic-state PV structures on thef-plane. We will show that both the linear eyewall evolution and the subsequent nonlinear behaviors caused by wave-wave interactions in the shallow-water-equation model are similar to those in the nondivergent barotropic model. The beta effect on eyewall instability is examined by comparing the vortex evolution in thef-plane and betaplane experiments, the details of which are discussed in section 5. Finally, a brief summary is given.

    2.Wave instability

    Here we choose the asymmetric balance (AB) shallowwater-equation model that allows divergence (Shapiro and Montgomery, 1993) to be included to address the wave instability in a two-dimensional shallow-water system. The linearized perturbation geopotential tendency equation deduced from the AB shallow-water equations in the cylindrical coordinates can be written as

    wherer,λ,denote radius, azimuthal angle, and tangential wind of the basic vortex, respectively; ?′is the perturbation geopotential;is the inertial parameter;fis the Coriolis parameter;is the basic-state PV;is thebasic-stateabsolute verticalvorticity;His thebasic-stategeopotentialheight;isthesquare of the inverse of local Rossby number; andis the basicstate geopotential. Assume that the perturbation geopotential satisfies a wave-type solution ?′= Φ(r)exp[i(nλωt)] (Φ is the wave amplitude;iis the imaginary number;nis the azimuthal wavenumber;ωis the time-frequency;tis time).Substitute it into Eq. (1), and one can get

    Given that the perturbation would vanish when the radial distance approaches infinity, then the term on the right-hand side of Eq. (3) should disappear when we integrate Eq. (3)fromr= 0 tor= ∞,

    Subtracting the imaginary part of Eq. (4) yields

    In the unstable condition, ωi≠0, then we have

    Another necessary condition, the Fjortoft theorem, can be obtained by subtracting the real part of Eq. (4), which gives

    Utilizing Eq. (6), Eq. (7) can be reduced to

    Since the left-hand side of Eq. (8) is always positive,andshould be anticorrelated on average (whento make the term on the right-hand side positive. Under the premise of the necessary condition, as dictated by the Charney-Stern theorem, that allows for the generation of counterpropagating VRWs, the Fjortoft theorem can ensure phaselocking between the pair of waves and lead to an enhanced perturbation amplitude (Hoskins, 1985). Note that the linear analysis above is only applicable to the early stage of exponential growth, which can profoundly impact the subsequent nonlinear evolution of TC eyewalls. Although the criteria only disclose the nature of waves in terms of global stability and can hardly indicate the most unstable modes, we will show how excitations of the most unstable modes depend on the basic-state PV structures through numerical simulations.

    3.Model description and experimental design

    The shallow-water-equation model used in this study includes the following equations in the Cartesian coordinates:whereu,v,Vdenote the zonal, meridional wind components, and wind vector, respectively;f= 2 Ω sinφ with φ

    being the latitude. The model domain has 3203 × 2403 grid points with a grid spacing of 3 km. Results within an area of 601 × 601 grid points centered at the vortex center are saved for analysis. Rigid lateral boundary conditions are used at the north and south boundaries, and cyclic conditions are used at the east and west boundaries. The combined scheme of Euler forward- and Matsuno time differencing is adopted for the model time integration with a time step of 9 s, while the simulated data are saved at 15-minute intervals. Note that our attention is given to the internal TC dynamics and the beta effect. Thus, we assume a quiescent environment.

    Observations suggest that intensifying TCs generally have very thin, high-PV annuli accompanied by nearly irrotational flow at the center (Schubert et al., 1999; Kossin and Eastin, 2001; Kossin and Schubert, 2001). To meet the two necessary conditions for barotropic instability presented in section 2, the model is initialized with an axisymmetric annulus of elevated PV (i.e., a PV ring), which has a radial distribution given as

    whereqr0is the maximum PV of the PV ring;rris the radial distance to the TC center;rois the radius at which the annulus is centered;lis the width of the annulus; the Coriolis parameterfis set at the latitude of 15°N; andHis the basic-state geopotential height. This radial profile indicates a nearly zero PV in the eye region with an increase in PV with increasing radius until a threshold radius given byro-l/2, then remains nearly constant before decreasing outward from another threshold radius given byro+l/2. Since PV and the inertial parameter are always positive in this case, such a structure satisfies the Charney-Stern theorem.

    Table 1 lists the parameters used for the PV distributions in the different experiments in this study. Because the model describes inviscid, unforced motion, PV is rearranged without any new PV generation in all experiments. Hence, it is of interest to see the redistribution of the same amount of PV under varying basic-state conditions. The total integrated PV over the model domain is to be nearly equal, by design,in all experiments; therefore, the maximum amplitudes of PV rings located at larger radii are smaller. All the model vortices have elevated PV rings of 6 km in width but have different central radii of 18 km, 45 km, and 90 km to mimic real unchanged.Referring to the segmentation function concerning the basic-state relative vorticity in Schubert et al. (1999),ξ1+ξ2and ξ2represent the relative vorticity in the axisymmetric eye (0 ≤r≤r1) and eyewall (r1≤r≤r2) regions, respectively. Two parametersγ= (ξ1+ξ2)/ξave(ξaveis the average relative vorticity over the region 0 ≤r≤r2) andδ=r1/r2are TCs of different eyewall sizes. Although more experiments can be done by further varying the widths, amplitudes, and ring sharpnesses, we will focus on a small parameter space to address the key scientific questions mentioned in section 1. Several three-point smoothing operations are performed to reduce the computational modes related to the sharp PV gradients at both the inner and outer edges of PV rings. As a result, the maximum amplitudes of PV rings in all experiments are reduced slightly, and the uniform PV distributions of eyewall are deconstructed due to the limited spatial resolution of 3 km (Fig. 1). Nevertheless, the general characteristics of the radial PV distributions are conserved, and the constraint of Charney-Stern theorem is still met since the PV structure with respect to the off-centered maximum is kept used to measure the contrast of relative vorticity between the central region and the enhanced ring and the degree of hollowness, respectively. In the shallow-water case, we replace relative vorticity with PV. Since the initial PV in the eye region has been set to zero, all the vortices haveγ=0 but differ in terms ofδ. Here the expression ofδis given asri/re,whereriis the radius of maximum basic-state PV gradient(positive) andreis the radius of the minimum basic-state PV gradient (negative). The smoothed PV profile is used to calculate this parameter shown in Table 1.

    Fig. 1. Radial distributions of the initial axisymmetric PV annuli (units: 0.1 PVU) given in Eq. (13) before smoothing(solid) and after smoothing (dashed) for three cases with an ro of 18 km (a), 45 km (b), and 90 km (c), respectively.

    Table 1. Description of the eight experiments in this study. The second and third columns give information on the maximum potential vorticity (PV) before and after smoothing; the fourth, fifth, and sixth columns present the widths of the PV rings, central radii of the location, and degree of hollowness, respectively. The last column shows the Coriolis parameters used in simulations. The first three experiments are conducted on the f-plane, while the last five are on the beta-plane.

    It is difficult to determine whether the basic-state flow fulfills the constraint of the Fjortoft theorem because the term on the right-hand side of Eq. (8) includes unknown parameters related to the perturbation. However, the integral of the part that only includes the basic-state parameters,given by can help determine the sign of the entire integral on the right-hand side of Eq. (8). Negative values (-3400 for 18 km; -1073 for 45 km; -308 for 90 km) for Eq. (14) suggest that all the experiments are likely to satisfy the requirement of the Fjortoft theorem.

    Given the basic-state PV distribution, the basic-state wind field is obtained by solving Eq. (12) iteratively,together with the geopotential height obtained from the nonlinear balance equation, which yields the relationship between the geopotential and streamfunction fields as given below

    whereψis the streamfunction. Equation (15) is numerically solved using a simple over-relaxation method based on the basic-state wind field (streamfunction) to obtain the initial geopotential height in Eq. (12). Ten iterations are deemed sufficient to obtain convergent numerical solutions.

    Threef-plane experiments are performed to serve as a basis for comparison with results from other experiments that include the beta effect (Table 1). Note that we initialize the model with a vortex structure that can enable the barotropic instability, which may develop freely with even small model numerical truncation errors. The opposite signed sharp radial PV gradients at radii ofro±l/2 permit the emergence of a pair of azimuthally counter-propagating Rossby edge-waves (Schubert et al., 1999). In contrast, the uniform distribution of PV within the interior and beyond the exterior boundary of the annulus is unfavorable for wave generation in our simulations.

    4.Development of instability and evolution on the f-plane

    The evolution of three kinds of vortices with different basic-state PV rings and associated wind structures on thefplane are examined first to explore the effects of barotropic instability under the weak-divergent condition. The azimuthal Hovm?ller diagrams of the perturbation PV at the radii where the most significant disturbances appear illustrate the overall development of asymmetries in each experiment(Fig. 2). The vortex in F18 first displays a short-period,nearly stationary, and weak wavenumber-4 (WN-4) asymmetry (Figs. 2a and 3b), which might be a consequence of the model’s intrinsic limitation of finite grid spacing since similar fluctuations have been detected in other experiments as well(Figs. 4b and 5b). After the initial slow development,barotropic instability becomes evident with the wave gradually intensifying to form a WN-2 structure (see the elliptical PV structure in Fig. 3e and the WN-2 geopotential height anomaly in Fig. 3f). However, due to the radial expansion of the nonzero PV area during the first few hours (Figs. 3a and 3c, fromt= 0 h tot= 6 h), the radial position of the growing VRW deviates from the initial location of maximum PV.In this case, when the VRW becomes significant, the maximum amplitude of the wave lies at a radius of about 30 km rather than at a radius of 18 km, justifying the azimuthal Hovm?ller diagram in Fig. 2a to be drawn at a radius of 30 km. Although an elliptically-shaped eyewall (corresponding to a WN-2 perturbation) survives till the end of the simulation, the perturbation gradually weakens and is anticipated to dissipate.

    Compared with F18, the predominant disturbance in F45 has greater growth and decay rates. In F45, a WN-4 asymmetry develops and experiences a rapid amplification during the first 2-4 h (Fig. 2b), then the asymmetry breaks down into four PV mesovortices (Fig. 4c). Later, a transition from the WN-4 perturbation to a lower-wavenumber perturbation occurs through the merger process (sincet= 10 h, see Fig. 2b).From then on, all disturbances remain insignificant, and ultimately, the initially hollow PV ring relaxes into a monopole,which is dynamically stable (Fig. 4g). Along with the axisymmetrization process, the central geopotential height drops from 1345 gpm at the beginning to 1263 gpm at the end of the simulation, implying an intensification of the vortex. It is consistent with the results of Kossin and Schubert (2001),who proposed that a large fall in central pressure is likely to occur after the merger of mesovortices in a two-dimensional nondivergent vortex.

    Fig. 2. Azimuthal Hovm?ller diagrams of perturbation PV (grey shading; units: 0.1 PVU) in F18 at the radius of 30 km (a), F45 at the radius of 45 km (b), and F90 at the radius of 90 km (c), respectively.

    Unlike previous experiments, two waves of different wavenumbers dominate in F90 simultaneously. The initial PV ring promises both WN-4 and WN-8 VRWs to grow exponentially. Their superposition enhances half of the anomalies while the other half is suppressed (Fig. 5d), leading to an asymmetric structure of inflow/outflow. The strengthened anticyclonic perturbations are located downstream (relative to the TC cyclonic circulation) relative to the strengthened cyclonic perturbations, while they are located upstream relative to the suppressed cyclonic perturbations. The outflow associated with the downstream weaker cyclonic perturbations, imposed on the enhanced anticyclonic perturbations,would be opposed and offset by the stronger inflow associated with the upstream enhanced cyclonic perturbations. The net effect is an inward advection of enhanced anticyclonic perturbations. A similar analysis can be applied to the other mesovortices.

    Consequently, those strengthened perturbations (both cyclonic and anticyclonic) would be advected inward,whereas those suppressed perturbations would be advected outward (Fig. 5f). Therefore, due to the vortex-vortex interaction, the eight pairs of cyclonic/anticyclonic anomalies split into two groups, with an inward-moving group of enhanced amplitudes and the other outward-moving group of reduced amplitudes. When the two groups of disturbances propagate into a proper position, the outer group would be translated inward by the inflow between the strengthened cyclonic and anticyclonic perturbations in the inner group, initiating the merger between the two groups as shown at the time of 15 h(Figs. 5g and 5h). The time-radius Hovm?ller diagrams of the WN-4 and WN-8 components (Fig. 6) demonstrate the evolution of the two unstable modes and the wave-wave interaction discussed above. Exponential growths of WN-4 and WN-8 waves occur simultaneously at 3 h (Figs. 6a and 6b),followed by the radial separation of the enhanced and suppressed WN-4 perturbations (Fig. 6a). The merger of the oppositely moving groups at 15 h leads to the re-intensification of the WN-4 perturbation, which can maintain itself for hours. After the radial movement and the merger process among mesovortices, the WN-4 perturbation keeps its maximum amplitude at a radius of 75 km instead of the initial radius of 90 km after 18 h (Fig. 6a). Though the WN-8 VRW also persists for a long time, its intensity is much weaker than the WN-4 VRW. Vortex development in F90 implies that a certain PV structure can enable VRWs with multiple wavenumbers emanating from barotropic instability, and the wave-wave interactions can cause complex nonlinear behaviors of mesovortices.

    Fig. 3. PV (left column; grey shading; units: 0.1 PVU) and perturbation geopotential height (right column;color shading; units: gpm) distributions over the domain of 45 × 45 (units: km), at 0 h, 6 h, 12 h, and 24 h in F18. Black contours, grey arrows in the left column are the total wind speeds and wind vectors, respectively.Grey arrows in the right column are the perturbation wind vectors. In each panel, the solid grey circle is drawn at a radius of 30 km, and the black velocity vector in lower right indicates the max wind. Note that contour levels are different in the right column for better graphic exhibition.

    Fig. 4. Same as Fig. 3 but over the domain of 60 × 60 (units: km) in F45. The grey circle in each panel is drawn at a radius of 45 km.

    Fig. 5. Same as Fig. 3 but over the domain of 120 × 120 (units: km) at 0 h, 3 h, 7 h, and 15 h in F90. The grey circle in each panel is drawn at a radius of 90 km.

    Results from numerical experiments based on the shallow-water-equation model demonstrate that the basic-state PV structure has vital influences on the excitation of barotropic instability and the subsequent nonlinear evolution. A less hollow vortex undergoes exponential growth at a smaller wavenumber with a relatively slower growth rate and gradually decays, ultimately resulting in a PV monopole. On the contrary, a hollower vortex exhibits a greater growth rate at a larger wavenumber and can maintain its coherent asymmetric structure for hours and days, which was also documented by Kossin and Schubert (2001) under a nondivergent framework. These results are in good agreement with those of Schubert et al. (1999), who employed a nondivergent barotropic model to study eyewall evolution driven by barotropic instability, suggesting that weak divergence has very limited effects on the linear instability of barotropic vortices. Note that although the shallow-waterequation system allows for the existence of inertia-gravity waves, disturbances in all experiments exhibit the typical features of VRWs that are dominated by the quasi-balanced dynamics, namely, the positive vorticity anomaly is in phase with the negative geopotential height perturbation(Figs. 3, 4, and 5).

    5.Influences of the beta effect

    5.1.Beta gyres and vortex movement

    In a barotropic TC-like vortex developing in a quiescent environment on a beta-plane, a nearly stationary (relative to the vortex) WN-1 perturbation (i.e., the beta gyres) over a large area will be present. Namely, a positive geopotential height anomaly with anticyclonic circulation appears to the east of the storm center, and a negative geopotential height anomaly with cyclonic circulation appears to the west, with southeasterly ventilation flow in between, which will advect the whole vortex northwestward. The translational speed and direction are sensitive to the wind structure at a radial range of 300 km to 1000 km from the vortex center (Fiorino and Elsberry, 1989). Overall, the northwestward displacements are quite similar among experiments B18, B45, and B90 except for the wobble of the TC vortex (a substantial oscillation of the vortex center) from the time around 60 h in B90 (Fig. 7b), which is, in fact, a manifestation of algebraic instability (Smith and Rosenbluth, 1990; Nolan and Montgomery, 2000; Nolan et al., 2001). The F90 vortex also experiences a TC wobble from 90 h of the simulation (Fig. 7a),whereas the vorticies in F18 and F45 hardly move during the simulation period (not shown).

    Fig. 6. Time-radius Homv?ller diagrams of the maximum PV amplitudes (grey shading; units: 0.1 PVU) of WN-4 (a) and WN-8 (b) waves in F90. Yellow lines denote the radius of 75 km and 90 km, respectively.

    5.2.Influences of the beta effect on eyewall evolution

    With the occurrence of the beta gyres and associated beta drift, there are little changes in the evolution of unstable VRWs (presented as elliptic or polygonal eyewalls) and basic flows on the beta-plane in B18 and B45 compared to their counterparts on thef-plane in F18 and F45 (Figs. 8a and 8b). In contrast, the development of the inner-core structure of the vortex in B90 is influenced by the beta effect.Although four enduring mesovortices in both F90 and B90 have undergone merger processes and eventually evolved into monopoles, there is a large discrepancy between the timing of the breakdown of asymmetric structures in the two experiments. The intensification of the basic circulation,which results from the axisymmetrization of asymmetric motion, takes place around 70 h in B90 but does not appear until 100 h in F90. Besides, a decrease of perturbation energy of the unstable modes (WN-4 and WN-8 waves)begins from the time of 60 h in B90 but about 30 hours later in F90 (Fig. 8c), which are consistent with the times when the vortices start to wobble substantially in the two experiments (Fig. 7).

    The wobble of the vortex center is part of the asymptotic solution for a WN-1 disturbance on a circular vortex with the existence of a local maximum in the basic-state angular velocity off the vortex center (Nolan and Montgomery,2000). Nevertheless, due to the active interaction between the eye and the eyewall in experiments of 18-km and 45-km PV-ring vortices, there is continuous transport of large angular momentum to the central area, resulting in the maximum angular velocity near the vortex center (Fig. 9a). This may explain the absence of algebraic instability in these experiments. In contrast, the asymmetries in F90 and B90 can maintain their intensity and structure for a long time without pronounced mixing of angular momentum within the vortices,thus a local maximum of angular velocity outside the eyewall may exist (Fig. 9a). According to Nolan and Montgomery(2000), who first applied the concept of algebraic instability in plasma physics (Smith and Rosenbluth, 1990) to twodimensional inviscid TC-like vortices, the asymptotic solution consists of three components: (1) a pseudo-mode,which accounts for the wobble displacement of the vortex center; (2) a growing mode (modal mode) with a secular growth rate that is proportional tot1/2, whose perturbation vorticity is exactly proportional to the basic-state radial vorticity gradient up to the radius of the angular velocity maximum;(3) the residual part, comprised of a collection of sheared VRWs, without the support of which the growing mode is unable to develop. Figure 8c presents the temporal evolution of maximum PV amplitudes of the WN-1 disturbances associated with algebraic instability until the axisymmetrization takes place in experiments of F90 and B90aThe distortion of the vortex caused by axisymmetrization process would present as a superposition of waves of different wavenumbers including the WN-1 wave in the Fourier space, which is not relevant to algebraic instability. Hence, we only show the time evolution of a WN-1 disturbance before the distortion begins to avoid any confusion.. The amplification of the WN-1 perturbations in F90 and B90 is accompanied by the emergence of significant trochoidal movements of the vortices and the weakening of the two unstable modes.Note that there are only slight changes in the strength of the basic circulation until a great surge occurs at 70 h in B90 and 100 h in F90. This suggests that the unstable modes, not the basic flow, are the energy source for the algebraic instability.

    The evolution of the algebraic instability is interfered with by the sheared VRWs. The energy transfer efficiency to the modal mode in a two-dimensional inviscid axisymmetric vortex can be computed by the formula given below(Nolan and Montgomery, 2000)

    whereE′is the perturbation energy; Ω is the basic-state angular velocity;are the radial and azimuthal velocities of the modal part, respectively;are the radial and azimuthal velocities of the residual part, respectively. Based on Eq. (16), the growth rate of the modal mode is affected by the radial gradient of the basic-state angular velocity in the area where algebraic instability emerges.Figure 9b shows the evolution of the averaged radial gradient of the basic-state angular velocity over the annular region between the radii of 75 km and 90 km, which coincides with the location of the largest amplitudes of the unstable VRWs in 90-km PV-ring vortices. The radial shear of Ω is nearly identical for the first 60 hours in F90 and B90, before it gradually sharpens and becomes negative in B90, but this gradual sharpening occurs 30 hours later in F90. Simultaneously, there is a tremendous drop in the energy of the unstable modes from 60 h in B90 and from 90 h in F90 (Fig. 8c),which justifies the role of the radial shear of rotation in wave-wave interactions.

    Fig. 7. Trajectories of vortex centers over the simulation time of 120 h. The vortex center is identified as the origin according to which the vortex has the largest azimuthally averaged winds. The left panel (a) displays the track of the vortex center in F90, in which black dots are the locations of the vortex center with a time interval of 30 min; grey curves indicate displacements of the vortex center during the time intervals; and the red dot [at the position of (0 km, 0 km)], yellow dot, and blue dot denote the vortex center positions at the times of 60 h, 90 h, and 91 h, respectively, indicating the appearance of significant wobble occurring at 90 h of simulation. The right panel (b) shows three trajectories for the vortices in B18 (blue dots), B45 (red dots),and B90 (green dots). Dots from faint to dark imply the time from the beginning to the end of simulations with an interval of 30 min. The yellow dot denotes the position of the vortex center at 60 h when the vortex begins to wobble substantially in B90.

    Fig. 8. Temporal evolution of the maximum PV amplitudes (units: 0.1 PVU) of the unstable modes and the basic circulation in each experiment. The maximum PV amplitude of the basic flow in each experiment corresponds to the left axis and is summed over the circular region within the radius of 120 km to cover the high-PV annulus. The maximum PV amplitudes of the unstable modes at the radius of the maximum intensity of perturbations (top panel: 30 km; middle panel: 45 km; bottom panel: 75 km) in each experiment correspond to the right axis. In the bottom panel, solid curves denote the maximum PV amplitudes of the basic flow, dashed curves denote the sum of maximum PV amplitudes of WN-4 and WN-8 waves at a radius of 75 km, and dotted curves denote the maximum PV amplitudes of WN-1 disturbances that are associated with the algebraic instability at a radius of 75 km. Since the modal mode in the asymptotic solution of algebraic instability is exactly proportional to the basic-state vorticity gradient, we draw the evolution of maximum PV amplitudes of the WN-1 disturbances at a radius of 75 km where the basic-state PV gradient is maximized in the 90-km PV-ring experiments. The magnitude at the upper-left corner in each panel corresponds to the left y-axis.

    The ability to sustain convection (high-PV structure)depends on the relative contributions of strain and vorticity.In a strain-dominated area, PV perturbations are short-lived and end up as elongated and deformed filaments of anomalous vorticity (Rozoff et al., 2006; Wang, 2008). In the 90-km PV-ring vortices, PV disturbances (VRWs) sheared by the radially decreasing angular velocity feed their energy to the modal mode of algebraic instability (starting from 60 h in B90, from 90 h in F90). This interaction could enhance the PV mixing between the eye and eyewall regions (Nolan and Montgomery, 2000), leading to the intensification of the azimuthal mean circulation (starting from 70 h in B90 and 100 h in F90) and eventually leading to the formation of PV monopoles. As shown above, although the destruction of the coherent structure of mesovortices also occurs in F90 wherefis constant across the domain, when the beta effect is involved, such a breakdown appears much earlier, as in B90.

    To further address the role of the beta effect in the nonlinear interaction between perturbations, two additional experiments on the beta-plane are performed. Both are centered at 15°N, but use different beta parameters, β at 60°N (B90-60)and β at 75°N (B90-75). The same central latitude is adopted to ensure that the vortex evolutions are influenced only by the meridional change of the Coriolis parameter rather than the Coriolis parameter itself. Vortices modeled with larger beta parameters (lower latitudes) are less likely to sustain coherent asymmetric structures. The azimuthally averaged radial shear of rotation in the region where VRWs exist becomes negative in each experiment but at different times (Fig. 9c). The radial profile of angular velocity within the annular region, between the radii of 75 km and 90 km in B90, experiences an earlier alteration from oscillating around a positive value to gradual decreasing. Such changes in experiments with smaller beta parameters occur in later hours. Subsequently, the mesovortices begin to merge at 70 h of simulation in B90, while the merger occurs approximately 10 and 20 hours later in B90-60 and B90-75, respectively. The quasi-steady structure in F90 maintains itself until about 100 h, when a final monopole state is nearly achieved in B90 (Fig. 10). The evolution of the azimuthal mean tangential circulation of a TC moving on a beta-plane has been previously studied by Smith et al. (1995). They proposed that the beta drift of an initially symmetric, cyclonic vortex in the Northern hemisphere would experience a weakening of its azimuthal mean tangential wind at most radii,with the tangential wind decreasing along the radius more rapidly than the decay rate of 1/(radius), which is mainly attributed to the radial flux of planetary vorticity and the meridional change of Coriolis parameter. Compared to those in thef-plane experiments, TC circulations in the betaplane experiments are subject to an extra, instantaneous flux of planetary vorticity due to the northwestward displacement of the vortex, which would spin down the azimuthal mean tangential winds, especially in the larger radii. The greater spindown rate of vortex circulation at larger radii will result in a negative gradient of basic-state angular velocity, and thus an earlier development of the algebraic instability, making the coherent structure of mesovortices less sustainable. The larger the beta parameter is, the faster the destruction will occur.

    Fig. 9. Radial profiles of the angular velocity in the experiments of 18-km PV-ring vortices, 45-km PV-ring vortices, F90 and B90 at 60 h when the vortex starts to wobble substantially in B90 (a), and the temporal evolution of the averaged gradient of basic-state angular velocity over the annular area between the radii of 75 km and 90 km (b). In the top panel, black and yellow lines correspond to the left axis, while blue lines correspond to the right axis. The magnitude at the upper left (right) corner corresponds to the left (right) yaxis in each panel, respectively.

    The above results demonstrate that the asymmetric forcing arising from the beta effect can perturb nonlinear TC evolutions. The beta effect may destabilize the coherent mesovortices resulting from the barotropic instability of a TC-like vortex with a large hollow PV structure. Subsequent vortex structural changes include sharpening the radial gradient of the basic-state angular velocity, which would deform VRWs and transfer energy from the sheared VRWs to the modal mode of algebraic instability. Differential rotation and innercore PV mixing by the algebraic instability can invoke mesovortex merger and vortex axisymmetrization. However,it must be clarified that the impact of the beta effect on the eyewall evolution is confined only to large-size TCs and is limited in less hollow TC-like vortices in which the resultant asymmetries are often axisymmetrized to form the monopole structure in a relatively short time scale, such as in 18-km and 45-km PV-ring vortices.

    6.Summary

    This study has examined the roles of barotropic instability and the beta effect on the evolution of TC-like vortices with different PV-ring structures using a shallow-water-equation model. Two necessary conditions of barotropic instability are deduced based on the two-dimensional AB shallowwater equations following Montgomery and Shapiro (1995).Three kinds of barotropically unstable PV rings with different degrees of hollowness are simulated, first on thef-plane, all of which exhibit polygonal shapes at different stages. More intense annuli with lower degrees of hollowness develop lower-wavenumber asymmetries with smaller growth rates.At the same time, hollower PV rings are more likely to feature higher-wavenumber disturbances experiencing transient growth and decay, consistent with those found by Schubert et al. (1999) in a nondivergent barotropic model. This implies that the weak divergence in a shallow-water-equation model would generally not alter the unstable nature of nondivergent motion although some differences may exist in detail. The total asymmetry excited by the barotropic instability can be decomposed into VRWs of different wavenumbers with various amplitudes, and nonlinear wave-wave interactions could substantially complicate the eyewall evolution.

    Although the beta effect is often neglected in studies that concern the internal dynamics of the inner-core region of a TC because of the comparatively small scale of the inner-core circulation, we have successfully demonstrated that the beta effect may undermine the stability of the coherent mesovortices in the eyewall of an extremely hollow, largesize PV ring. In this case, the asymmetric mesovortices can maintain themselves for an extended period in a nondivergent barotropic frame with a constant Coriolis parameter (Kossin and Schubert, 2001). However, observed mesovortices in the eyewalls of real TCs often have a much shorter lifetime,i.e., several to ten hours. In this study, we have shown that the beta effect could lead to an earlier breakdown of coherent mesovortices by modifying the radial shear of rotation and initiating energy transfer among the waves and also between the waves and the basic vortex flow. Therefore, the quasisteady mesovortices may undergo much earlier merger and axisymmetrization processes on a beta-plane, which suggests that the beta effect may be one of the forcings that shorten the lifetime of mesovortices in real strong TCs.

    Acknowledgements. This study was supported in part by the National Natural Science Foundation of China under Grant No.41730960 and in part by NSF grant AGS-1834300.

    淫妇啪啪啪对白视频| 搡老熟女国产l中国老女人| 国产日本99.免费观看| 男女床上黄色一级片免费看| 国产精品1区2区在线观看.| 一区二区三区国产精品乱码| 91大片在线观看| 成年免费大片在线观看| 久久久精品国产亚洲av高清涩受| 欧美日韩中文字幕国产精品一区二区三区| 午夜免费观看网址| 啦啦啦免费观看视频1| 啦啦啦观看免费观看视频高清| 制服诱惑二区| 一级片免费观看大全| 亚洲中文av在线| 亚洲第一青青草原| 亚洲久久久国产精品| 亚洲第一av免费看| 亚洲专区国产一区二区| 伦理电影免费视频| 中亚洲国语对白在线视频| 在线视频色国产色| 亚洲中文字幕一区二区三区有码在线看 | 熟女少妇亚洲综合色aaa.| 99久久国产精品久久久| 亚洲七黄色美女视频| 男女下面进入的视频免费午夜 | 亚洲av电影不卡..在线观看| 国产成人影院久久av| 国内揄拍国产精品人妻在线 | 成人特级黄色片久久久久久久| av在线天堂中文字幕| 国产欧美日韩精品亚洲av| 麻豆成人av在线观看| 国产1区2区3区精品| 特大巨黑吊av在线直播 | 啦啦啦韩国在线观看视频| 两人在一起打扑克的视频| 久久人妻福利社区极品人妻图片| 久久亚洲精品不卡| 精品久久久久久,| 日韩成人在线观看一区二区三区| 老司机在亚洲福利影院| 热re99久久国产66热| 人妻丰满熟妇av一区二区三区| 日韩欧美在线二视频| 精品久久蜜臀av无| 伦理电影免费视频| 欧美色欧美亚洲另类二区| 一本久久中文字幕| 18禁美女被吸乳视频| 波多野结衣高清作品| 国产一区在线观看成人免费| 久久久久久九九精品二区国产 | 脱女人内裤的视频| 午夜成年电影在线免费观看| 无人区码免费观看不卡| 视频在线观看一区二区三区| 天天添夜夜摸| 国产精品,欧美在线| 亚洲美女黄片视频| 身体一侧抽搐| 欧美 亚洲 国产 日韩一| 成人永久免费在线观看视频| 看免费av毛片| 男男h啪啪无遮挡| 看片在线看免费视频| 久久午夜综合久久蜜桃| svipshipincom国产片| 色av中文字幕| www.精华液| 色av中文字幕| 色av中文字幕| 精品久久蜜臀av无| 欧美成人免费av一区二区三区| 午夜老司机福利片| 精品熟女少妇八av免费久了| 久久中文字幕一级| 男女之事视频高清在线观看| 久久久国产精品麻豆| 天堂动漫精品| 国产黄a三级三级三级人| 男女做爰动态图高潮gif福利片| 美女大奶头视频| 成人三级黄色视频| 18禁观看日本| 欧美日本视频| 亚洲在线自拍视频| 亚洲av电影不卡..在线观看| 99久久精品国产亚洲精品| 成年免费大片在线观看| 两人在一起打扑克的视频| 很黄的视频免费| 成人三级做爰电影| xxxwww97欧美| 人人妻人人澡人人看| 欧美乱妇无乱码| 亚洲av熟女| 欧美日韩乱码在线| 精品久久久久久久毛片微露脸| 亚洲av日韩精品久久久久久密| 欧美激情高清一区二区三区| 满18在线观看网站| 99久久无色码亚洲精品果冻| 久久精品91无色码中文字幕| 午夜福利视频1000在线观看| 亚洲成av人片免费观看| 成人手机av| 男女之事视频高清在线观看| 久热这里只有精品99| 一进一出抽搐gif免费好疼| 侵犯人妻中文字幕一二三四区| 久久人妻av系列| 最近最新中文字幕大全电影3 | 正在播放国产对白刺激| 国产精品野战在线观看| 成人永久免费在线观看视频| 中文字幕av电影在线播放| 在线观看舔阴道视频| 在线国产一区二区在线| 9191精品国产免费久久| 中文字幕精品亚洲无线码一区 | 国产高清videossex| 在线播放国产精品三级| 中文在线观看免费www的网站 | www.精华液| 久久青草综合色| 两性午夜刺激爽爽歪歪视频在线观看 | 91成年电影在线观看| 欧美zozozo另类| 99久久国产精品久久久| 国产麻豆成人av免费视频| aaaaa片日本免费| 夜夜夜夜夜久久久久| 美女扒开内裤让男人捅视频| 在线播放国产精品三级| 国产激情欧美一区二区| 国产精品国产高清国产av| 欧美精品亚洲一区二区| 变态另类成人亚洲欧美熟女| 国产一区二区三区视频了| 2021天堂中文幕一二区在线观 | 久久香蕉激情| aaaaa片日本免费| 看免费av毛片| 亚洲av五月六月丁香网| 久久精品91蜜桃| 日韩成人在线观看一区二区三区| 亚洲国产精品sss在线观看| 草草在线视频免费看| 国产av不卡久久| 性色av乱码一区二区三区2| 国产精品 国内视频| 青草久久国产| 婷婷丁香在线五月| 黄色丝袜av网址大全| svipshipincom国产片| 日本三级黄在线观看| 亚洲欧洲精品一区二区精品久久久| 观看免费一级毛片| 高潮久久久久久久久久久不卡| 精品国产乱子伦一区二区三区| 啦啦啦韩国在线观看视频| 免费观看精品视频网站| 窝窝影院91人妻| 色播在线永久视频| 在线观看舔阴道视频| 免费在线观看视频国产中文字幕亚洲| 老司机深夜福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 女人高潮潮喷娇喘18禁视频| 一区福利在线观看| 岛国视频午夜一区免费看| 欧美性猛交黑人性爽| 一边摸一边抽搐一进一小说| 久久午夜亚洲精品久久| 大型av网站在线播放| 日韩有码中文字幕| 精品无人区乱码1区二区| 日本在线视频免费播放| 精品第一国产精品| 香蕉久久夜色| 亚洲一区二区三区色噜噜| 女人高潮潮喷娇喘18禁视频| 久热这里只有精品99| 91成年电影在线观看| av在线天堂中文字幕| 一本久久中文字幕| 别揉我奶头~嗯~啊~动态视频| 校园春色视频在线观看| 国产又黄又爽又无遮挡在线| 精品久久久久久久人妻蜜臀av| 亚洲专区国产一区二区| 亚洲 国产 在线| 一区福利在线观看| 波多野结衣高清作品| 一夜夜www| 亚洲五月婷婷丁香| xxxwww97欧美| 50天的宝宝边吃奶边哭怎么回事| 麻豆一二三区av精品| 国产精品一区二区精品视频观看| 午夜福利高清视频| 婷婷六月久久综合丁香| 精品久久久久久久毛片微露脸| 国产精品 国内视频| 久久久精品国产亚洲av高清涩受| 色尼玛亚洲综合影院| 成人精品一区二区免费| 侵犯人妻中文字幕一二三四区| 中文字幕人成人乱码亚洲影| 日韩欧美一区视频在线观看| 亚洲国产欧洲综合997久久, | 女人被狂操c到高潮| 九色国产91popny在线| 无限看片的www在线观看| 亚洲五月色婷婷综合| 日韩精品青青久久久久久| 国产免费av片在线观看野外av| videosex国产| 午夜福利在线观看吧| 欧美黄色淫秽网站| 亚洲天堂国产精品一区在线| 免费一级毛片在线播放高清视频| 色老头精品视频在线观看| 一级毛片精品| 成人一区二区视频在线观看| 亚洲自偷自拍图片 自拍| 亚洲成人久久性| 精品乱码久久久久久99久播| 99国产综合亚洲精品| 国产极品粉嫩免费观看在线| 制服人妻中文乱码| 亚洲午夜理论影院| 亚洲国产欧美日韩在线播放| 色综合亚洲欧美另类图片| 久久久精品欧美日韩精品| 国产国语露脸激情在线看| 看片在线看免费视频| 少妇粗大呻吟视频| 香蕉av资源在线| 听说在线观看完整版免费高清| 99riav亚洲国产免费| 老司机福利观看| 亚洲一区二区三区不卡视频| 日本免费一区二区三区高清不卡| 香蕉av资源在线| 国产精品自产拍在线观看55亚洲| 中国美女看黄片| 99热这里只有精品一区 | 亚洲电影在线观看av| 人人妻,人人澡人人爽秒播| 久久久水蜜桃国产精品网| 视频在线观看一区二区三区| √禁漫天堂资源中文www| 婷婷亚洲欧美| 搡老岳熟女国产| 国产精品精品国产色婷婷| 一本一本综合久久| 亚洲av电影不卡..在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 很黄的视频免费| 国产精品电影一区二区三区| 变态另类成人亚洲欧美熟女| 在线永久观看黄色视频| 久久中文看片网| 69av精品久久久久久| 女人被狂操c到高潮| 国产99久久九九免费精品| 精品福利观看| 午夜免费成人在线视频| 国产高清激情床上av| 国产精品精品国产色婷婷| 国产熟女午夜一区二区三区| 午夜福利成人在线免费观看| 久久久久免费精品人妻一区二区 | 日韩 欧美 亚洲 中文字幕| 国产高清有码在线观看视频 | 天堂动漫精品| 久久久久国内视频| 午夜久久久在线观看| 无限看片的www在线观看| 亚洲精品在线观看二区| 亚洲精品国产精品久久久不卡| 国产视频内射| 脱女人内裤的视频| 成人特级黄色片久久久久久久| 精品国产一区二区三区四区第35| 亚洲无线在线观看| 国产国语露脸激情在线看| 黄片播放在线免费| 黑人巨大精品欧美一区二区mp4| 国产真人三级小视频在线观看| 亚洲欧美一区二区三区黑人| 久久久久久亚洲精品国产蜜桃av| 真人做人爱边吃奶动态| 一进一出抽搐gif免费好疼| 他把我摸到了高潮在线观看| 亚洲狠狠婷婷综合久久图片| 国产国语露脸激情在线看| 两个人看的免费小视频| 亚洲色图av天堂| 国产高清激情床上av| 色精品久久人妻99蜜桃| 色尼玛亚洲综合影院| 久久香蕉精品热| 久热这里只有精品99| 欧美zozozo另类| 老司机福利观看| 国产亚洲精品一区二区www| 12—13女人毛片做爰片一| 午夜福利欧美成人| 久久久久久久精品吃奶| 国产v大片淫在线免费观看| 国产免费男女视频| 亚洲国产精品sss在线观看| 日韩欧美免费精品| 亚洲人成伊人成综合网2020| 久久精品国产综合久久久| 最好的美女福利视频网| 国产三级黄色录像| 香蕉国产在线看| 国产aⅴ精品一区二区三区波| 老司机在亚洲福利影院| 欧美在线黄色| 啪啪无遮挡十八禁网站| 99热这里只有精品一区 | 女同久久另类99精品国产91| 亚洲国产精品999在线| 免费观看人在逋| 91国产中文字幕| 久久人妻av系列| 午夜两性在线视频| 国产蜜桃级精品一区二区三区| 夜夜爽天天搞| 黄色丝袜av网址大全| 一边摸一边做爽爽视频免费| 91麻豆精品激情在线观看国产| 18禁裸乳无遮挡免费网站照片 | 色播在线永久视频| 久久这里只有精品19| 看黄色毛片网站| 成人国产一区最新在线观看| 色综合站精品国产| 午夜精品久久久久久毛片777| 淫妇啪啪啪对白视频| 夜夜看夜夜爽夜夜摸| 欧美黑人精品巨大| 免费搜索国产男女视频| 熟女少妇亚洲综合色aaa.| 欧美日韩精品网址| 国产一卡二卡三卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品av麻豆狂野| 搡老岳熟女国产| 亚洲av熟女| 制服人妻中文乱码| 国产一区二区三区视频了| 久久香蕉国产精品| 久久久精品国产亚洲av高清涩受| 亚洲国产欧美网| 一a级毛片在线观看| a在线观看视频网站| 中文在线观看免费www的网站 | 757午夜福利合集在线观看| 美女大奶头视频| 国产高清有码在线观看视频 | 视频区欧美日本亚洲| 免费在线观看亚洲国产| 欧美乱色亚洲激情| 天天躁狠狠躁夜夜躁狠狠躁| 成人免费观看视频高清| e午夜精品久久久久久久| 99国产精品一区二区三区| 久久精品aⅴ一区二区三区四区| 国产成人系列免费观看| 欧美成人一区二区免费高清观看 | 曰老女人黄片| 波多野结衣高清作品| 麻豆国产av国片精品| 婷婷亚洲欧美| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全电影3 | 免费电影在线观看免费观看| 一本久久中文字幕| 窝窝影院91人妻| 午夜日韩欧美国产| 曰老女人黄片| 亚洲熟妇熟女久久| 可以免费在线观看a视频的电影网站| 黄网站色视频无遮挡免费观看| 中文字幕另类日韩欧美亚洲嫩草| 国产aⅴ精品一区二区三区波| 精品国产超薄肉色丝袜足j| 亚洲av日韩精品久久久久久密| 夜夜夜夜夜久久久久| 亚洲av中文字字幕乱码综合 | 91国产中文字幕| 久久久久久久久免费视频了| 亚洲国产精品sss在线观看| 丝袜人妻中文字幕| 丁香六月欧美| 国产精品久久久久久亚洲av鲁大| 日日干狠狠操夜夜爽| 国产激情久久老熟女| 999久久久精品免费观看国产| 美女高潮到喷水免费观看| 老司机福利观看| 18禁观看日本| 免费人成视频x8x8入口观看| 观看免费一级毛片| 久久婷婷成人综合色麻豆| 露出奶头的视频| 国产久久久一区二区三区| ponron亚洲| а√天堂www在线а√下载| 日韩欧美国产在线观看| ponron亚洲| av超薄肉色丝袜交足视频| 在线国产一区二区在线| 久久中文看片网| 午夜激情福利司机影院| 国产精品1区2区在线观看.| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 成人特级黄色片久久久久久久| 成人免费观看视频高清| 色av中文字幕| 久久精品91蜜桃| 黄色视频,在线免费观看| 午夜久久久在线观看| 18禁观看日本| 欧美成人性av电影在线观看| 午夜福利成人在线免费观看| 熟女电影av网| 精品久久久久久久毛片微露脸| 精品日产1卡2卡| 亚洲欧美精品综合一区二区三区| av有码第一页| 99久久综合精品五月天人人| 亚洲精品中文字幕在线视频| 久久中文字幕一级| 欧美大码av| 亚洲真实伦在线观看| 精品久久久久久久末码| 波多野结衣巨乳人妻| 在线永久观看黄色视频| 国产麻豆成人av免费视频| 久久久久国产精品人妻aⅴ院| 国产私拍福利视频在线观看| 色综合亚洲欧美另类图片| 满18在线观看网站| 久久热在线av| 最好的美女福利视频网| 青草久久国产| 天天一区二区日本电影三级| 天堂√8在线中文| svipshipincom国产片| 亚洲精品一区av在线观看| 国产精品久久视频播放| 亚洲中文日韩欧美视频| 国产av一区二区精品久久| 非洲黑人性xxxx精品又粗又长| 老熟妇乱子伦视频在线观看| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 日韩中文字幕欧美一区二区| 午夜福利免费观看在线| 国产人伦9x9x在线观看| 91成人精品电影| 后天国语完整版免费观看| 亚洲最大成人中文| 波多野结衣高清作品| 国产亚洲欧美精品永久| 欧美+亚洲+日韩+国产| 少妇裸体淫交视频免费看高清 | 欧美黑人欧美精品刺激| 国产97色在线日韩免费| 黑人巨大精品欧美一区二区mp4| 成人亚洲精品一区在线观看| 国内久久婷婷六月综合欲色啪| 亚洲全国av大片| 国产色视频综合| 久久精品亚洲精品国产色婷小说| 国产午夜精品久久久久久| 欧美成人免费av一区二区三区| 久久久久久人人人人人| 亚洲精华国产精华精| 久久久久免费精品人妻一区二区 | 免费看十八禁软件| 岛国在线观看网站| 亚洲人成77777在线视频| 国产精品一区二区免费欧美| 久久婷婷人人爽人人干人人爱| 午夜激情福利司机影院| 一级黄色大片毛片| 国产亚洲精品久久久久5区| 看片在线看免费视频| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站| 久久香蕉精品热| 欧美成人性av电影在线观看| 国产亚洲精品综合一区在线观看 | 国产亚洲精品第一综合不卡| 午夜免费激情av| 不卡av一区二区三区| 老汉色av国产亚洲站长工具| 日本三级黄在线观看| 国产av在哪里看| 1024手机看黄色片| 又黄又爽又免费观看的视频| 女性生殖器流出的白浆| 天堂影院成人在线观看| 久久国产乱子伦精品免费另类| 99在线视频只有这里精品首页| 欧美不卡视频在线免费观看 | 国产一区二区在线av高清观看| 欧美一级a爱片免费观看看 | 男人舔女人下体高潮全视频| 无人区码免费观看不卡| 操出白浆在线播放| 国产黄片美女视频| 18禁美女被吸乳视频| 91九色精品人成在线观看| 亚洲国产欧洲综合997久久, | 欧美日本亚洲视频在线播放| 国产真实乱freesex| 精品不卡国产一区二区三区| av视频在线观看入口| 亚洲午夜理论影院| 国产精品二区激情视频| 日韩欧美在线二视频| 久久久国产欧美日韩av| 黄色片一级片一级黄色片| 亚洲精品中文字幕在线视频| 色精品久久人妻99蜜桃| 天天一区二区日本电影三级| 久久人妻福利社区极品人妻图片| 99精品在免费线老司机午夜| 欧美黑人欧美精品刺激| 夜夜躁狠狠躁天天躁| 国产蜜桃级精品一区二区三区| 国产激情久久老熟女| 国产不卡一卡二| 狂野欧美激情性xxxx| 亚洲专区中文字幕在线| 国产av又大| 精品少妇一区二区三区视频日本电影| 国产成人精品久久二区二区91| 中文字幕精品免费在线观看视频| 叶爱在线成人免费视频播放| 人人妻人人看人人澡| 久久久国产精品麻豆| 国产激情偷乱视频一区二区| 十八禁网站免费在线| 精品第一国产精品| 亚洲第一青青草原| 国产高清videossex| 一进一出好大好爽视频| 中文字幕最新亚洲高清| 一区二区日韩欧美中文字幕| 91麻豆av在线| 久久狼人影院| 中出人妻视频一区二区| 99re在线观看精品视频| 中出人妻视频一区二区| 国产精品久久久av美女十八| 国产视频内射| 久久久国产欧美日韩av| 日本撒尿小便嘘嘘汇集6| 自线自在国产av| 亚洲人成电影免费在线| ponron亚洲| 国产成人一区二区三区免费视频网站| 国产精品综合久久久久久久免费| 国产极品粉嫩免费观看在线| 自线自在国产av| 欧美日韩福利视频一区二区| 悠悠久久av| 麻豆成人午夜福利视频| 午夜两性在线视频| 国产精华一区二区三区| 国产精品亚洲一级av第二区| 久久天躁狠狠躁夜夜2o2o| 国产精品99久久99久久久不卡| 免费人成视频x8x8入口观看| 少妇熟女aⅴ在线视频| 搡老妇女老女人老熟妇| 欧美乱码精品一区二区三区| 国产精品精品国产色婷婷| 女人爽到高潮嗷嗷叫在线视频| 成人亚洲精品一区在线观看| 一级毛片高清免费大全| 欧美日韩中文字幕国产精品一区二区三区| 成人精品一区二区免费| 淫秽高清视频在线观看| 亚洲国产日韩欧美精品在线观看 | 国产精品久久久av美女十八| 99久久综合精品五月天人人| av在线播放免费不卡| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦免费观看视频1| 99热这里只有精品一区 | 亚洲黑人精品在线| 女同久久另类99精品国产91| a在线观看视频网站| 亚洲 欧美 日韩 在线 免费| 香蕉国产在线看| 免费在线观看影片大全网站| 一级作爱视频免费观看| 男人舔女人的私密视频| 亚洲成a人片在线一区二区| 成人三级黄色视频| 欧美日韩亚洲国产一区二区在线观看|