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

    Initiation and Evolution of Long-Lived Eastward-Propagating Mesoscale Convective Systems over the Second-Step Terrain along Yangtze-Huaihe River Valley※

    2022-04-02 05:51:48YuanchunZHANGJianhuaSUNRuyiYANGandRuoyunMA
    Advances in Atmospheric Sciences 2022年5期

    Yuanchun ZHANG, Jianhua SUN, Ruyi YANG, and Ruoyun MA

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

    2Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519080, China

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

    4Public Meteorological Service Center, China Meteorological Administration, Beijing 100044, China

    ABSTRACT Based on the previous statistical analysis of mesoscale convective systems (MCSs) over the second-step terrain along Yangtze-Huaihe River Valley, eight representative long-lived eastward-propagating MCSs are selected for model-based sensitivity testing to investigate the initiation and evolution of these types of MCSs as well as their impact on downstream areas.We subject each MCS to a semi-idealized (CNTL) simulation and a sensitivity (NOLH) simulation that neglects condensational heating in the formation region.The CNTL experiment reveals convection forms in the region downstream of a shortwave trough typified by persistent southwesterly winds in the low-to midtroposphere.Upon merging with other convective systems, moist convection develops into an MCS, which propagates eastward under the influence of midtropospheric westerlies, and moves out of the second-step terrain.The MCS then merges with pre-existing local convection over the plains; the merged convection reinforces the cyclonic wind perturbation into a mesoscale vortex at 850 hPa.While this vortex moves eastward to regions with local vortex at 850 hPa, another vortex at 925 hPa is also intensified.Finally,the vortices at 850 and 925 hPa merge together and develop into a mesoscale convective vortex (MCV).In contrast, MCSs fail to form and move eastward in the NOLH experiment.In the absence of eastward-propagating MCSs, moist convection and mesoscale vortices still appear in the plains, but the vortex strength and precipitation intensity are significantly weakened.It is suggested the eastward-propagating MCSs over the second-step terrain significantly impact the development and enhancement of moist convection and vortices in the downstream areas.

    Key words: mesoscale convective systems, second-step terrain, mesoscale convective vortex, numerical sensitivity simulation

    1.Introduction

    Mesoscale convective systems (MCSs) are well-known significant weather systems that produce hazardous weather, such as heavy rain, lightning, and hail (Maddox,1980; Tollerud and Collander, 1993; Zheng et al., 2013).With improvements in satellite remote sensing and radar detection, many experts explored the spatial and temporal distribution of MCSs (Maddox, 1980; Augustine and Howard,1991; Laing and Fritsch, 1997).Following Maddox (1980),numerous studies have evaluated the statistical features of MCSs by applying size and duration criteria of the cold cloud shields based on infrared (IR) satellite images in various regions (Augustine and Caracena, 1994; Laing and Fritsch, 1997; Anderson and Arritt, 1998; Mathon and Laurent, 2001; Jirak et al., 2003; Rafati et al., 2017; Zheng et al., 2008; Yang et al., 2015, Yang et al., 2020).Because of the complex terrain distribution in China, there are regional differences in the spatial and temporal features of MCSs in China (Zheng et al., 2008; Qi and Zheng, 2009).Most previous studies of MCSs in China have focused on the Tibetan Plateau (TP) and Yangtze-Huaihe River Valley(YHRV).Mesoscale convective systems (MCSs) over the TP show clear temporal variation with a peak frequency of occurrence in the afternoon.Over the Sichuan Basin, MCSs mostly occur during the nighttime, and convection over the plains can also initiate during the afternoon and early morning (Jiang and Fan, 2002; Zheng et al., 2008; Bai et al.,2011; Fu et al., 2013; Hu et al., 2016).

    However, few studies have paid attention to MCSs on the second-step terrain in China.The second-step terrain usually refers to the high mountain ranges from the Da Hinggan Mountains over northeast China (including the Yanshan Mountains, the Taihangshan Mountains, the Loess Plateau, and the Mongolian Plateau over north China and the Qinling Mountains over central China) to the Yunnan-Guizhou Plateau over southwest China (partly shown in Fig.1a).Some studies have found that severe convective systems appear over the YHRV to the east of the second-step terrain(Meng et al., 2013; Zheng et al., 2013; Zheng and Sun,2016).Previous simulation studies show that these convective systems can be traced back to convective systems over the second-step terrain (He and Zhang, 2010; Sun and Zhang, 2012; Zhang and Sun, 2017, Zhang et al., 2018).Many studies found that the initiation mechanisms of MCSs are different when the underlying surface conditions are different (Velasco and Fritsch, 1987; Trier and Parsons, 1993;Laing and Fritsch, 1997; Zipser et al., 2006; Zheng et al.,2008; Qi and Zheng, 2009, Rasmussen and Houze, 2011;Luo et al., 2014; He et al., 2018; Zhuo et al., 2012 ).Some case studies found that latent heating over the second-step terrain enhances the intensity, longevity, and eastward movement of mesoscale vortices over the eastern plains (Fu et al.,2017).In addition, the combination of the mountain-plain solenoid (MPS), due to the heating difference between the second-step terrain and its eastern plains, and the enhanced nocturnal low-level jet (LLJ) favored convection initiation east of the second-step terrain (Zhang et al., 2014, 2018).

    The above mentioned previous findings are concluded from case studies.Therefore, it is crucial to explore and summarize the common mechanisms responsible for the initiation and development of MCSs over the second-step terrain.A recent statistical study identified and tracked MCSs around the second-step terrain along the middle reaches of the YHRV (the region is shown in the solid brown box in Fig.1b) and classified these MCSs into eastward-propagating and quasi-stationary types (Yang et al., 2019, 2020).It is found that eastward-propagating MCSs, which originate from the second-step terrain, enhance the precipitation in downstream regions.Based on these studies, eight longlived eastward-propagating MCSs over the second-step terrain will be selected objectively in the present study to explore their initiation and detail their impact on downstream regions.Section 2 introduces the objective method used to select the typical MCS cases.Section 3 elaborates on the specific model design of the semi-idealized simulation.The initiation and evolution characteristics of this type of long-lived eastward-propagating MCSs are analyzed in section 4.The results of the sensitivity experiment are discussed in section 5.Finally, conclusions and discussion are presented in section 6.

    2.Case selection

    A total of 316 eastward-propagating MCSs that form over the second-step terrain are divided into four categories according to their key features (such as their initiation, trajectories, and termination, Yang et al., 2019).One of the four categories for MCSs (i.e., a set of 55 MCSs) consists of those that had the longest lifetime and made the largest contribution to warm-season heavy rainfall (the tracks of these 55 MCSs are shown in Fig.1b) in downstream regions of second-step terrain.To explore the common mechanisms of initiation and development of these long-lived MCSs and their impact on downstream regions, the correlation coefficients related to the geopotential height of the 500 hPa surface in the key region for these 55 MCSs are calculated to select the MCS cases which initiate and develop under similar synoptic circulations.The key region at 500 hPa is determined as the region covering 100°–120°E, 26.5°–36.5°N (the solid brown rectangle in Fig.1b).The correlation coefficients for these 55 long-lived eastward-propagating MCSs at their formation time in the key region are calculated using the following formula:

    Fig.1.(a) Distribution of terrain elevation higher than 500 m (units: m, color shading), d01 and d02 stand for the domains of the simulations in the text.(b) Distribution of terrain height higher than 500 m (units: m, color shading), the solid brown rectangle represents the key region used for correlation coefficient calculation, and the brown dashed rectangle outlines the region in which MCSs form (formation region) and the saturation latent heating is turned off.The purple trajectories are the tracks of 55 MCSs, which initiate over the second-step terrain and propagate eastward to the downstream regions.

    Based on calculations using the above formula, we identified eight representative MCSs among the long-lived eastward-propagating 55 MCSs, for which the correlation coefficient between any two exceeds 0.6.In addition, all eight MCSs have lifetimes longer than 15 h, with precipitation falling over regions east of 115°E (Table 1).At 500 hPa (Fig.2), case 7 is influenced by a cold vortex in the south and airflow around the periphery of the Western Pacific subtropical high (WPSH), whereas the remaining cases are characterized by a westerly trough located east of the Tibetan Plateau (TP) and a strengthened WPSH with its ridge extending westward to 110°E.Mesoscale convective systems(MCSs) form in the region of southwesterly winds associated with the trough and WPSH, where warm air advection and positive vorticity advection provide favorable environmental conditions for convective initiation.Although the synoptic weather patterns at 500 hPa are generally the same among the eight MCSs, the low-level weather patterns differ.An examination of the 850 hPa weather pattern for the eight MCSs suggests that all MCSs form in the mid-eastern mountainous region of the second-step terrain along YHRV with only a slight difference in the exact locations (figures not shown).The precipitation distribution associated with the eight MCSs extends from the Qinling Mountains to downstream areas.The tracks of the MCSs with the largest spatial influence span the middle and lower reaches of YHRV and the East China Sea.Among the eight MCSs, cases 3 and 4 produced the most accumulated precipitation (over 30 mm), while cases 7 and 8 are associated with less precipitation (Fig.3).The rain bands associated with the eight MCSs generally cover the region east of the second-step terrain to the middle and lower reaches of YHRV; therefore, they are representative of the MCSs over this region and are appropriate to use for further simulation and mechanism studies.

    Table 1.Time and location of the formation and termination of eight MCSs.

    3.Experiment configuration

    Due to the limited spatial and temporal resolution of conventional observational data, numerical simulations are used in the present study to investigate the influence of longlived eastward-propagating MCSs on precipitation over the second-step terrain along YHRV in the downstream region.The primary goal of the current study is to examine the impact of the second-step terrain on the initiation of MCSs.Some previous studies have employed semi-idealized numerical simulations to reveal the common impact of terrain on the MCSs and the associated mechanism of diurnal variation (Trier et al., 2010; Sun and Zhang, 2012; Chen et al.,2016).Specifically, semi-idealized numerical simulations are carried out using the Weather Research and Forecasting model (WRF, V3.8; Skamarock and Klemp, 2008).The initial and boundary conditions are averaged from the selected eight MCSs provided by the Climate Forecast System Reanalysis (CFSR) dataset from the National Centers for Environmental Prediction (NCEP) with a horizontal grid spacing of 0.5° × 0.5° at 6-hour time intervals.Composite initial and boundary conditions are obtained based on the following scheme.The formation time of the MCS is defined as Tn(0), with n representing the sequence number for eight real cases.The formation time is defined herein as the first time when the contiguous cold cloud shield of the MCS(TBB ≤ 221 K) reaches an area of at least 5000 km2.Since the temporal resolution of the reanalysis data is 6 h, the field at 6 hours before the formation time is denoted as Tn(?6), and that for 6 hours after the formation time is denoted as Tn(6), and so on.The simulation covers the period from 12 hours before to 24 hours after the formation time, which gives a simulation period of 36 h, so the initial conditions are given by the average field of the eight MCSs cases at T n(?12), and the boundary conditions are obtained from the composite fields of eight MCSs at Tn(?6), Tn(0),Tn(6), T n(12), T n(18), and T n(24), respectively.The simulation starts at 12 hours before the formation time of the eight MCSs (defined as t = 00 h), which is 0200 Local Standard Time (LST, LST = UTC + 8; LST is used hereafter).

    The model domain is illustrated in Fig.1a.Two oneway nested domains with a 9 km and 3 km horizontal grid spacing are used, which yields 530 lon.× 390 lat.and 1039 lon.× 796 lat.grids, respectively, in the two domains.A total of 50 sigma layers are set in the vertical with the model top at 50 hPa.For both the inner and outer domains, the WRF single-moment 6-class microphysics scheme (WSM6; Hong and Lim, 2006), the Rapid Radiative Transfer Model(RRTM) longwave radiation scheme (Mlawer et al., 1997),the Dudhia shortwave radiation scheme (Dudhia, 1989), the Noah land surface model scheme (Chen and Dudhia, 2001),the Yonsei University (YSU) planetary boundary layer scheme (Hong et al., 2006) are used.The Kain-Fritsch cumulus parameterization (Kain, 2004) is applied to the outer domain.The above setting is referred to as the CNTL experiment.

    Fig.2.The synoptic weather pattern at 500 hPa for the eight eastward-propagating MCS cases at their formation time: geopotential height (solid black contours, units: gpm).temperature (red dashed contours, units: K), wind(blue arrow, units: m s?1), and relative vorticity (color shading, units: 10–6 s–1).The red dots represent the formation locations of MCSs.Grey shading represents terrain elevations > 500 m.(a) Case 1 at 0500 LST 28 June 2000; (b) case 2 at 1700 LST 28 June 2001; (c) case 3 at 1400 LST 22 June 2002; (d) case 4 at 0200 LST 23 June 2016; (e) case 5 at 1700 LST 01 July 2003; (f) case 6 at 1700 LST 21 June 2006; (g) case 7 at 2300 LST 13 August 2010; (h) case 8 at 1700 LST 8 July 2012.

    Fig.3.Distribution of accumulated precipitation (color shading, units: mm) during the lifetimes of the eight MCSs.Grey shading represents terrain elevations > 500 m.

    Latent heating is a key factor for the organization, development, and maintenance of moist convection.It is thought to be difficult for convection to organize into MCSs if inadequate latent heating release is present after convection initiation.Therefore, a sensitivity experiment with the diabatic heating in the formation region (106°–113°E, 28°–35°N)turned off (NOLH) is performed to study the impact of terrain on long-lived eastward-propagating MCSs and precipitation over downstream areas.

    4.Results from the CNTL simulation

    4.1.Verification of CNTL results

    The composite circulation of the eight MCSs at 500 hPa at 0200 LST shows that the shortwave trough is located at the eastern edge of the Tibetan Plateau and Sichuan Basin.At the same time, the WPSH covers the region south of YHRV and South China.Southwesterly winds on the periphery of WPSH favor the eastward propagation of MCSs to the middle and lower reaches of YHRV (figures not shown).In the lower levels of the troposphere (i.e., 925 hPa),cyclonic wind shear is maintained over the lee side of the second-step terrain.The associated cyclonic curvature in the wind field continually intensifies, which finally gives rise to the formation of a mesoscale vortex over the boundaries between Henan, Anhui, and Shandong Provinces by 0200 LST the following day.The vortex intensifies and propagates eastward to downstream regions (Figs.4a, c, e, g).The synoptic circulation pattern at 500 hPa in the CNTL simulation is similar to that of the composite results from eight real cases in reanalysis data.Still, the simulated WPSH in lower latitudes is slightly stronger than that of the reanalysis data, thereby exerting control over areas further west(the simulated 5880 gpm isohypse extends west of 110°E).In addition, the simulated trough in the westerlies to the east of the TP is stronger (figure not shown).The simulated flow at 925 hPa reproduces vortex formation over the eastern edge of the second-step terrain, its development, and eastward propagation.The formation location of the vortex in the simulation agrees well with that from composite results,but the simulated vortex is stronger (Figs.4b, d, f, h).Figure 5 compares long-lived eastward-propagating trajectories between the eight simulated MCSs and eight real MCS cases.The simulated contiguous cold cloud of the MCS reaches 5000 km2at 1530 LST, which is delayed by about 1.5 h compared to the composite results (the MCS forms at 1400 LST).In general, the simulation reproduces the average locations of the eight MCSs during their formation, eastward propagation, and dissipation stages (Fig.5a).The simulated MCS initiates at the eastern edge of the second-step terrain (boundaries between Shaanxi and Hubei Provinces),about 1ofurther north of the average initiation location of the real cases.After formation, the MCS propagates eastward out of the second-step terrain and influences the southern and central regions of Henan Province, northern areas of Anhui Province, and the YHRV before finally decaying to the east into the sea.The simulated accumulated rainfall distribution from 12 hours before until 12 hours after the MCS formation is basically consistent with that from the eight real MCSs, except that the simulated precipitation maxima are slightly greater in the simulation compared to the real cases (Fig.5b).

    The above comparison between the CNTL and composite results from real cases clearly indicates that the CNTL simulation reproduces the eastward-propagating trajectories, rainbelt coverage, and background synoptic circulation relatively well during MCS evolution.Therefore, the initiation and evolution of long-lived, eastward-propagating MCSs over the second-step terrain along the YHRV are explored based on the simulation results.Additionally, sensitivity experiments are designed to further analyze the latent heating effects on downstream mesoscale systems.

    4.2.Formation and evolution of long-lived eastwardpropagating MCS

    Based on the simulation results, the formation and evolution of the MCS are characterized by four stages.Namely,the MCS formation stage (1330–1530 LST), propagation stage over the second-step terrain (1500–1900 LST), mesoscale convective vortex (MCV) formation stage (1900–2330 LST), maintenance stage (2330–0900 LST), and the dissipation stage (0830–1400 LST).

    4.2.1.MCS formation stage (1330-1530 LST)

    The 850 hPa convergence between northwesterly winds to the north of the Qingling Mountains and southeasterly winds to the southeast of the Qingling Mountains enhances the vertical motion and cloud-top heights.Convection with cloud-top TBB colder than 221 K (C1; Fig.6a) first appears over Qingling Mountains in the southeastern areas of Shanxi Province (109.5°E, 33.3°N) at 1330 LST, and the corresponding rainfall intensity is more than 20 mm h–1(Fig.6a).Meanwhile, the intensified cyclonic circulation over the eastern edge of the second-step terrain (~112°E) leads to the continuous southward advancement of northeasterly winds along the periphery of the circulation, which results in the formation of a convergence region at 1430 LST (Fig.6b).The northwesterly and southeasterly winds associated with this convergence region trigger a new area of convection to the north of C1 (Fig.6b; hereafter C2, 110°E, 33.75°N).The new area of convection, C2, gradually intensifies after initiation and merges with C1 about one hour later, leading to the formation of a mesoscale convective system (hereafter MCS1) with a cold cloud shield coverage more extensive than 5000 km2(Fig.6c; the MCS criteria in Yang et al.,2020).The rainfall rate associated with MCS1 reached 50 mm h–1.

    4.2.2.Propagation stage over the second-step terrain(1530-1900 LST)

    After formation, MCS1 intensifies and propagates eastward over the second-step terrain.In the meantime, a convergence region between the northeasterly and southeasterly winds is located to the east of MCS1 (Fig.6c).When MCS1 moves into the convergence region, it merges with preexisting isolated convection therein and further intensifies (Fig.6d).Accordingly, several precipitation centers appear with maximum rainfall intensities reaching 100 mm h–1(Fig.6d).At 1700 LST (Fig.7a), while MCS1 propagates eastward on the eastern edge of the second-step terrain, the corresponding wind convergence intensified.The convergence at the northern end of the southwesterly LLJ intensifies and promotes the development of MCS1 at 1900 LST (Fig.7c).

    Fig.4.Left column: composite relative vorticity (color shading, units: 10–5 s–1) and flow field at 925 hPa of the eight MCSs using CFSR reanalysis data at (a) 1400 LST, (c) 2000 LST, (e) 0200 LST, and(g) 0800 LST, respectively.Right column: simulated relative vorticity (color shading, units: 10–5 s–1)and flow field (right) at 925 hPa at (b) 1330 LST, (d) 2000 LST, (f) 0200 LST, and (h) 0800 LST.Grey shading represents terrain elevations >500 m.

    Fig.5.(a) The real propagation trajectories (solid blue lines) of the eight MCSs and the simulated path of the MCS in CNTL simulation (solid red line).Dots represent MCS locations at each time.Brown shading represents terrain elevations > 500 m.Purple box covers the MCS formation regions.(b) Precipitation amounts from 6 h before and 6 h after the formation of MCS averaged over the eight MCSs (red contours, units: mm) and precipitation amount in CNTL simulation (color shading, units:mm).Grey shading represents terrain elevations > 500 m.

    Fig.6.The CNTL results: wind field (wind barb, 4 m s–1; speeds with > 8 m s–1 are indicated in blue) at 850 hPa,hourly precipitation (color shading, units: mm h–1), and boundary of TBB = 221 K in CNTL simulation (purple lines).Grey shading represents terrain elevations > 500 m.

    4.2.3.MCV formation stage (1900-2300 LST)

    After 1900 LST, the intensified wind perturbation at 850 hPa gradually develops into mesoscale vortices.The initial system, MCS1, moves out of the second-step terrain and merges with local convection in downstream regions at 2000 LST (the merged convection is referred to as MCS hereafter).Meanwhile, because of upward motion appearing at~1.5 km (figures not shown), wind convergence appears over the plains at 850 hPa (114.4°E, 34.5°N), located to the northeast of the MCS.With the enhancement of upward motion after 2000 LST, the wind convergence develops into a cyclonic vortex (Fig.7e; V1).The strong wind convergence (113.5°E, 34°N) at 925 hPa first appears at ~2030 LST, 30 minutes later than its appearance at 850 hPa.The convergence gradually evolves into a new mesoscale vortex(V2) after 2130 LST (Figs.7f–h).

    Fig.7.CNTL results: wind field (units: m s–1) at 850 hPa (left column) and 925 hPa (right column), TBB (shading,units: K).The grey shading represents the terrain elevations > 500 m.

    After the MCS and its associated vortex propagate eastward out of the second-step terrain and move into the V1 region, the merged V1 intensifies into a mature mesoscale vortex (Fig.7g).With the enhancement of V1, the strong upward motion extends both upward and downward (figures not shown).The vortex center of V2 at 925 hPa shifts below the V1 center at 2230 LST (Fig.7h).Subsequently,V2 at 925 hPa combines with the merged V1 at 850 hPa.Finally, the mature MCV forms at lower levels (below 700 hPa)after 2330 LST.

    4.2.4.Maintenance stage (2300-0900 LST) and dissipation stage (0900-1400 LST)

    The mature MCV structure below 700 hPa is maintained over the YHRV until 0900 LST.However, the MCS accelerates eastward and decays during this MCV maintenance stage.Consistent with the enhanced cyclonic circulation, the precipitation related to the MCV significantly increases before reaching a maximum at 0700 LST (figures not shown), which causes heavy rainfall over the middle reaches of the Yangtze River.The MCV begins to dissipate after 0900 LST; however, cyclonic wind shear is maintained to the east of the mesoscale vortex.The MCV gradually changes its shape into a zonally-oriented elliptical structure.Moreover, the MCV begins to dissipate along with a weakened daytime LLJ, and the related precipitation consequently decreases (figure not shown).

    The above analysis divides the evolution period of the MCS and its related MCV into five stages.1) The MCS formation stage (1330–1530 LST) starts with convection (C1)triggered by the convergence between northwesterly and southeasterly wind at 850 hPa over Qingling Mountains.After propagating eastward and merging with local convection over the mountains, C1 develops into MCS1.2) The MCS1 moves eastward with its related wind convergence during the propagation stage over the second-step terrain(1530–1900 LST).3) The MCV formation stage (1900 to 2300 LST) occurs as MCS1 propagates eastward out of the eastern edge of the second-step terrain and merges with preexisting convective systems over the plains, then the merged MCS and its corresponding wind perturbation intensifies at 1930 LST.Meanwhile, the locally triggered vortex(V1) at 850 hPa and V2 at 925 hPa appear over the plains northeast of the MCS after 2230 LST.The 850 hPa vortex(V1) intensifies after the MCS-related vortex moves into the region of V1.The enhanced upward motion of V1 promotes the eastward displacement of the vortex center of V2 at 925 hPa.The combined effects of V1 and V2 initiate the formation of a mature MCV structure (below 700 hPa) over the YHRV.4) The MCV is maintained and causes heavy rainfall over the middle reaches of the Yangtze River during 2300~0900 LST before dissipating between 0900 and 1400 LST.

    4.3.Environmental conditions at convection initiation

    The above analysis demonstrates that the eastwardpropagating MCS influences the vortex and precipitation over the downstream regions.However, the above analysis focuses on the stage when the MCS has attained a relatively strong intensity (i.e., after the TBB falls below 221 K).This section will discuss and compare the environmental conditions prior to and after convection initiation.Based on the simulated composite radar reflectivity, a single cell (SC1) with reflectivity larger than 10 dBZ appears over the Qinling Mountains (109.27°E, 33.28°N) at 0930 LST(Fig.8a).Twenty minutes later, the reflectivity increases to the criterion of convection initiation (35 dBZ, May and Ballinger, 2007; Zhang et al., 2016).When the convection over upstream regions propagates eastward, SC1 strengthens through its merger with the upstream convection.The cloud-top height of SC1 increases and its TBB decreases.While the TBB above the maximum reflectivity is below 221 K at 1330 LST (Fig.8e), the convection develops into C1 (Fig.6a).The environmental conditions before and after SC1 initiation and the initiation mechanism will be further explored in the following.

    Favorable thermal and moisture conditions contribute to convection initiation.At 500 hPa, a westerly propagating trough appears east of TP three hours before SC1 initiation(0730 LST) and is maintained.The 5880 gpm isohypse extends westward to 110°E and controls the regions south of 30°N.Convection initiates under southwesterly winds ahead of the trough and on the northwestern periphery of WPSH (Fig.9b).After convection initiation, the steering winds gradually turn westerly, which favors the eastward propagation of convection.The strong southwesterly wind west of the WPSH transports abundant water vapor and warm airflow into the initiation regions near 700 hPa(Figs.9c, f), which causes the relative humidity to exceed 90%, high equivalent potential temperature, and strong instability.At lower levels (~850 hPa), the southwesterly wind rotates cyclonically over the lee side of the secondstep terrain, transporting abundant warm and moist flow(Fig.9e).The wind shear between southeasterly and southwesterly wind over the eastern edge of the second-step terrain (about 109.3°E, 33°N) provides favorable dynamic conditions for convection initiation.Temperatures increase to south of the convection initiation location while temperatures decrease to the north of the convection initiation location.The convection tends to initiate over the regions with the highest temperature gradient (Fig.9f).After convection initiation, the downstream temperature gradient is maintained; cooler temperatures over the northern regions continue to extend eastward, and the temperature over southern regions also increases.Six hours after convection initiation,the temperature gradient in the convection region is still large.

    Figure 10 shows the soundings before and after convection initiation at the center of SC1 (109.27°E, 33.27°N).The air is close to saturation, and the convective available potential energy (CAPE) reaches 901 J kg–1at 0920 LST(Fig.10a).Ten minutes later, the upward motion at low levels overcomes the convective inhibition (CIN), and the ascending air saturates upon reaching the lifting condensation level (LCL) near 830 hPa.The wind veers from low to middle levels consistent with warm advection in this layer.In addition to promoting ascent, the warm advection associated with the southerly wind component promotes atmospheric instability.With the upward extension of the saturation layer, the moist layer deepens, and the CAPE increases to 1080 J kg–1(Fig.10b).The strong upward motion initiates SC1 at 1020 LST.After convection initiation and the associated release of instability, the CAPE decreases to 332 J kg–1(figure not shown).The equivalent potential temperature increases at 0920 LST (1 h before SC1 initiation), causing the enhancement of the equivalent potential temperature gradient.The maximum gradient appears at SC1 initiation time, 1020 LST (Fig.10c).After 0940 LST, the condensation and associated latent heat release promotes the strong upward motion under a strong CAPE condition(Fig.10e).The cloud water mixing ratio rises sharply and extends to upper levels.When the cloud water mixing ratio maximum is lifted to about 650 hPa (figure not shown), the convection initiates at 1020 LST.The vertical cross-section of horizontal wind shear shows that southeasterly wind is maintained at 850 hPa before 0920 LST, then the wind gradually turns to easterly and changes to northeasterly at 0950 LST.Meanwhile, the wind field at 800 hPa turns from southerly to southeasterly.The southwesterly wind is maintained above 750 hPa.The wind speed gradually increases from 750 to 600 hPa before SC1 initiation yet decreases across this layer after SC1 initiation.The above analysis found that the wind veers with height, confirming that warm advection is present in the initiation region.The wind shear increases at middle and lower levels before SC1 initiation,conducive to convective development.After initiation, the enhanced wind shear at middle and upper levels favors the development and organization of convection (Coniglio et al., 2006; Chen et al., 2015).

    Fig.8.Radar composite reflectivity (shading, units: dBZ ) in CNTL from the 7.5 to 11.5 simulation hours.SC1 symbolizes a single-cell and C1 stands for the convective complex.The grey shading represents terrain elevations > 500 m.

    5.Sensitivity experiment results

    Based on the CNTL simulation, we chose to perform a sensitivity experiment by turning off the latent heating related to condensation.The experiment is designed to explore the impact of condensational heat release on eastward-propagating MCSs and convection and precipitation in downstream regions.The brown dashed box in Fig.1 illustrates the regions in which the latent heat of condensation is turned off in a way such that the condensational heating and evaporative cooling in the CNTL setting are substituted by the only evaporative cooling in the NOLH experiment(Figs.11a, b).

    The accumulated precipitation from MCS1 formation(1530 LST) until the end time of NOLH simulation of the two simulations is compared in Figs.12a, b.In the CNTL run, accumulated precipitation greater than 100 mm occurred along a northeast-southwest axis.The associated precipitation maxima was mainly located over the southern regions of Henan Province and northern regions of Jiangsu Province.After turning off the latent heating in the NOLH experiment, the precipitation west of 114°E decreases sharply, and the rainfall over the eastern regions still maintains an east-west orientation.Still, the precipitation maxima maintains itself over the boundaries between Shandong and Jiangsu Provinces.

    Fig.9.The CNTL results at 0730 LST (left column) and 1030 LST (right column).Panels (a) and (b) are the wind field (wind barb, 5 m s–1), geopotential height (solid blue lines, units: gpm), and temperature (red dashed lines, units:K) at 500 hPa.The grey shading represents terrain elevations > 500 m.Panels (c) and (d) show the wind field (wind barb, 5 m s–1), equivalent potential temperature (solid red line, units: K), and relative humidity (shading, units: %) at 700 hPa.The grey shading represents terrain elevations > 3000 m.Panels (e) and (f) show the wind field (wind barb,units: 5 m s–1), geopotential height (solid blue lines, units: gpm), and temperature (shading, units: K) at 850 hPa.The grey shading represents terrain elevations > 1500 m.

    When the latent heating is turned off over the formation region in the NOLH experiment, the appearance of the first radar reflectivity is delayed, and the reflectivity intensity is weaker.The convection which forms in southern regions of Shanxi Province weakens, and the merged convection does not develop into new MCS1 (figure not shown).Therefore, an MCS is not triggered over the second-step terrain and obviously cannot go on to influence the downstream regions.This result can be used to infer the effects of MCS1 on the weather of downstream regions.

    In the CNTL simulation, the MCV over downstream regions maintains itself at lower levels (850 and 925 hPa) at 2000 LST.While the latent heating is turned off in NOLH,the vortex at lower levels still forms at 2000 LST.To compare the vortex intensity in the two simulations, a key region, centered on the vortex with a 300 km × 300 km horizontal extent, is chosen for analysis.The average relative vorticity is calculated within the key region for the period from the formation to maturity of the vortex (Figs.12c–f).Without the impact of an eastward-propagating MCS, the vortex intensity in NOLH is much weaker than that in CNTL,especially during the first half-period of MCV duration(2000–0130 LST).Moreover, the average vorticity difference between the two simulations in the key region at 925 hPa is greater than that at 850 hPa (Figs.11c, d).It means that the eastward-propagating MCS greatly impacted the development and maintenance of the downstream low-level vortex, especially at the 925 hPa level.

    Fig.10.Soundings at the SC1 location in CNTL before initiation, temperature profile (red line), dew point temperature profile (red dashed line), and stratification curve (blue dashed line).Panel (a) is at 0920 LST and (b)1010 LST.The physical parameters of evolution at SC1 initiation location with time are shown in (c)–(e).Panel(c) is the equivalent temperature (units: K), (d) is diabatic heating (units: 10–3 K s–1), (e) is the vertical velocity(units: m s–1), and (f) is horizontal wind field (wind barb, 5 m s–1, colors represent wind speeds).

    Fig.11.Latent heating (units: 10–3 K s–1) distribution at 700 hPa over the formation areas in CNTL (a) and NOLH(b).The grey shading represents terrain elevations > 500 m.Averaged vorticity over the key areas of the vortex at (c)850 and (d) 925 hPa of CNTL (red lines) and NOLH (blue lines).

    Throughout the simulation period of NOLH, heavy rainfall was generally lacking west of 114°E along the YHRV,except for the dispersed heavy rainfall distribution (Figs.12a, b).This result implies that an eastward-propagating MCS directly influences the related precipitation in downstream regions.Additionally, because of the weak lowerlevel vortex east of the second-step terrain, the precipitation associated with the vortex also decreases.The MCV at 850 hPa matures at 2230 LST, causing heavy rainfall over YHRV (Fig.12c).As for rainfall in the CNTL simulation,the rainbelt covers a larger area, and the accumulated precipitation is greater (Figs.12c–e).In the later stages of the MCV, the precipitation in the two simulations mainly appears over the vortex center and its eastern section.Still,the rainfall in the sensitivity experiment is weaker and is less intense (Figs.12d–f).

    6.Conclusions and discussion

    Based on the statistical analysis of eastward-propagating MCSs over the second-step terrain along YHRV (Yang et al., 2019, 2020), this study investigates the formation and development of a typical type of long-lived eastwardpropagating MCS.We performed semi-idealized sensitivity testing using a numerical simulation from eight representative MCS cases.The conclusions are as follows.

    First, eight eastward-propagating MCSs were selected based on similarities in their synoptic weather patterns and convection initiation conditions.Next, composite meteorological fields from the eight MCSs are used as initial and boundary conditions for WRF simulation.The semi-idealized simulation successfully reproduces the formation of convection over the second-step terrain, its eastward propagation, its merger with pre-existing convection in the downstream areas,and the associated precipitation pattern.Results from the CNTL experiment indicate that the complete evolution of the MCS and its related MCV includes five stages: the MCS formation stage (1330–1530 LST), the propagation stage over the second-step terrain (1530–1900 LST), the MCV formation stage (1900–2300 LST), the maintenance stage(2300–0900 LST) and the dissipation stage (0900–1400 LST).Convection forms in regions with southwesterly winds in the low-to-midtroposphere commonly observed downstream of a 500 hPa shortwave trough.In this case, the shortwave was located east of the TP and on the northwest periphery of the WPSH, where abundant warm and moist air provides favorable moisture conditions for convection initiation.Convection initiation results from the release of CAPE triggered by the low-level convergence of southeasterly and northeasterly winds.After initiation, convection gradually propagates eastward under the influence of westerlies in the middle troposphere.During this time, moist convection develops and intensifies into MCS due to the low-level convergence and unstable stratification in the lower troposphere.

    Fig.12.Accumulated precipitation (units: mm) from the 11.5 to 36 simulation hours of CNTL (a) and NOLH (b).Wind field at 850 hPa (wind barb, 4 m s–1; speeds with > 8 m s–1 are indicated in blue) and precipitation (shading,units: mm h–1) from the 20.5 and 33 simulation hours, (c) and (e) are for CNTL, (d) and (f) are for NOLH.The grey shadings represent terrain elevations > 500 m.

    Figure 13 shows the conceptual model of the impact of an eastward-propagating MCS over the second-step terrain on the evolution of an MCV over downstream regions.While an MCS propagates out of the eastern edge of the second-step terrain and merges with the convection systems over the plains, the corresponding wind perturbation intensifies into a vortex at 850 hPa.The mesoscale vortex moves eastward, and the enhanced southwesterlies on the southeastern periphery of this vortex gradually promote enough local wind convergence to develop a local vortex(V1) at 850 hPa.This intensified local vortex merges with the leeside vortex at 925 hPa (V2) and finally develops into a mature MCV.The presence of enhanced nocturnal LLJ enables the merged convection to develop further.The MCV then intensifies and moves eastward and subsequently merges with the continuously strengthening vorticity centers on its eastern (downstream) side.At this time,the MCV reaches the mature stage with notable precipitation along its southern flank where convergence associated with the LLJ is most strongly focused.

    Fig.13.The conceptual model of the impact of an eastward-propagating MCS over the second-step terrain on the evolution of MCV over the downstream regions, including stages: eastward-propagation of MCS1 out of the secondstep terrain, merger with the local convection system, formation and maintenance of MCV.

    Results from the sensitivity experiment with diabatic heating in the formation region turned off indicate that MCSs do not form and move eastward over the second-step terrain when no diabatic heating is available.In the absence of eastward-propagating MCSs, convective and mesoscale vortices still exist in the plains to the east of the second-step terrain along YHRV, but the vortex strength and precipitation intensity weaken markedly.This result indicates that the eastward movement of these long-lived MCSs has a significant impact on the development and enhancement of convection and vortices in the downstream areas.

    Based on a semi-idealized simulation from composite fields, the present study reveals the formation, development,and impact on downstream systems of a typical type of eastward-propagating MCSs.Future studies will analyze real cases and simulate them to gain a deeper understanding of the impact of second-step terrain along YHRV on the initiation of MCSs and the relevant mechanisms related to mesoscale vortices induced by eastward-propagating MCSs.

    Acknowledgements.This article is dedicated to Prof.Fuqing ZHANG, who greatly contributed to our long-term international cooperation on mesoscale meteorology.We sincerely appreciate Prof.Fuqing ZHANG for all the suggestions, discussions, and help regarding the relationship between MCSs and mesoscale vortexes east of the second-step terrain in China.This research was supported by the National Key R&D Program of China (Grant No.2018YFC1507200) and the National Natural Science Foundation of China (Grant No.41975057).

    亚洲av日韩在线播放| 大陆偷拍与自拍| 亚洲va日本ⅴa欧美va伊人久久| 免费女性裸体啪啪无遮挡网站| 黄色 视频免费看| 美女主播在线视频| 亚洲视频免费观看视频| 99久久99久久久精品蜜桃| 午夜激情久久久久久久| 大型av网站在线播放| 欧美午夜高清在线| 欧美久久黑人一区二区| 久久免费观看电影| kizo精华| 亚洲伊人色综图| 建设人人有责人人尽责人人享有的| 亚洲欧洲日产国产| 中文字幕另类日韩欧美亚洲嫩草| 国产精品久久久人人做人人爽| 18禁观看日本| 美女国产高潮福利片在线看| 啦啦啦在线免费观看视频4| 亚洲国产欧美在线一区| 如日韩欧美国产精品一区二区三区| 国产成人精品久久二区二区91| 国产亚洲精品久久久久5区| 免费女性裸体啪啪无遮挡网站| 国产精品久久久av美女十八| 一级黄色大片毛片| 人妻久久中文字幕网| 国产精品久久久久久人妻精品电影 | 亚洲情色 制服丝袜| 国产黄频视频在线观看| 大型av网站在线播放| 日韩欧美一区视频在线观看| 久久久久久免费高清国产稀缺| 丁香六月欧美| 一本—道久久a久久精品蜜桃钙片| 男女边摸边吃奶| 涩涩av久久男人的天堂| 夜夜骑夜夜射夜夜干| 国产成人一区二区三区免费视频网站| 国产精品av久久久久免费| 搡老岳熟女国产| 在线 av 中文字幕| 黑人欧美特级aaaaaa片| 亚洲综合色网址| 免费观看a级毛片全部| 国产有黄有色有爽视频| 精品国产一区二区三区四区第35| 欧美 日韩 精品 国产| 美国免费a级毛片| 午夜激情久久久久久久| 亚洲精华国产精华精| 日韩制服丝袜自拍偷拍| 欧美日韩av久久| 大片免费播放器 马上看| 首页视频小说图片口味搜索| 丝袜人妻中文字幕| 一个人免费在线观看的高清视频| 99久久人妻综合| 国产精品1区2区在线观看. | 成人三级做爰电影| 18禁美女被吸乳视频| 国产主播在线观看一区二区| 不卡一级毛片| 国产成人免费无遮挡视频| 色尼玛亚洲综合影院| 欧美人与性动交α欧美精品济南到| 国产一区二区三区综合在线观看| 在线观看人妻少妇| 亚洲熟女毛片儿| 欧美激情久久久久久爽电影 | 99国产极品粉嫩在线观看| 亚洲专区字幕在线| 高潮久久久久久久久久久不卡| 夜夜夜夜夜久久久久| 亚洲人成77777在线视频| 精品亚洲成a人片在线观看| 日韩一区二区三区影片| 午夜91福利影院| 国产精品一区二区在线不卡| 99热网站在线观看| 欧美精品一区二区免费开放| 国产免费福利视频在线观看| 国产一区二区激情短视频| 人妻一区二区av| 欧美av亚洲av综合av国产av| 欧美+亚洲+日韩+国产| 日本黄色视频三级网站网址 | 久久精品国产亚洲av香蕉五月 | 黄色毛片三级朝国网站| 成人影院久久| 欧美另类亚洲清纯唯美| 在线观看免费午夜福利视频| 高清毛片免费观看视频网站 | 中文字幕色久视频| 满18在线观看网站| 又大又爽又粗| 一边摸一边抽搐一进一小说 | 中文亚洲av片在线观看爽 | 中亚洲国语对白在线视频| 国产亚洲精品第一综合不卡| 不卡一级毛片| 国产无遮挡羞羞视频在线观看| 久久国产精品影院| 高清欧美精品videossex| 91大片在线观看| 午夜福利在线观看吧| 男女之事视频高清在线观看| 国产免费现黄频在线看| 黑丝袜美女国产一区| 国产精品 国内视频| 美女高潮到喷水免费观看| av欧美777| 日韩制服丝袜自拍偷拍| 亚洲成人免费电影在线观看| 在线十欧美十亚洲十日本专区| 成人精品一区二区免费| 欧美中文综合在线视频| 亚洲人成电影观看| 亚洲九九香蕉| 久久天堂一区二区三区四区| 午夜福利欧美成人| 精品亚洲成a人片在线观看| 亚洲av成人不卡在线观看播放网| 日本黄色视频三级网站网址 | 一进一出抽搐动态| 老鸭窝网址在线观看| 90打野战视频偷拍视频| 人人妻人人爽人人添夜夜欢视频| 99re在线观看精品视频| 久久久久久久大尺度免费视频| 中文亚洲av片在线观看爽 | 亚洲av电影在线进入| 美女高潮喷水抽搐中文字幕| 日韩免费av在线播放| 91成人精品电影| 国产精品久久久av美女十八| 12—13女人毛片做爰片一| 国产精品99久久99久久久不卡| 国产精品成人在线| 巨乳人妻的诱惑在线观看| 亚洲国产欧美日韩在线播放| 国产精品国产高清国产av | 麻豆乱淫一区二区| 无限看片的www在线观看| 色尼玛亚洲综合影院| 欧美国产精品一级二级三级| 国产国语露脸激情在线看| 亚洲五月婷婷丁香| 久久久久久久久久久久大奶| 精品福利永久在线观看| 交换朋友夫妻互换小说| 国产成人精品久久二区二区91| 一边摸一边抽搐一进一小说 | 久久国产亚洲av麻豆专区| 国产一区二区 视频在线| 久久人妻av系列| 美女扒开内裤让男人捅视频| 精品欧美一区二区三区在线| 大片电影免费在线观看免费| 久久久久视频综合| 亚洲精品国产一区二区精华液| 亚洲成a人片在线一区二区| 久久久久国内视频| 叶爱在线成人免费视频播放| 在线观看一区二区三区激情| 久久精品国产综合久久久| 真人做人爱边吃奶动态| 国产成人免费无遮挡视频| 精品卡一卡二卡四卡免费| 女性被躁到高潮视频| 中文字幕精品免费在线观看视频| 黑人操中国人逼视频| 9191精品国产免费久久| 久久人人爽av亚洲精品天堂| 国产成人啪精品午夜网站| 国产精品电影一区二区三区 | 老司机影院毛片| 国产精品 国内视频| 亚洲色图 男人天堂 中文字幕| 成人黄色视频免费在线看| 人妻一区二区av| 久久久久久久国产电影| 国产一区有黄有色的免费视频| 久久av网站| 亚洲久久久国产精品| 大香蕉久久网| 精品国产乱码久久久久久男人| 久久99一区二区三区| 一进一出抽搐动态| 国产97色在线日韩免费| 久久人人爽av亚洲精品天堂| 又黄又粗又硬又大视频| 999久久久精品免费观看国产| 国产野战对白在线观看| 19禁男女啪啪无遮挡网站| 欧美精品人与动牲交sv欧美| 又大又爽又粗| 亚洲人成伊人成综合网2020| 淫妇啪啪啪对白视频| 老司机深夜福利视频在线观看| 婷婷成人精品国产| 国产淫语在线视频| 国产在线免费精品| 麻豆av在线久日| a在线观看视频网站| 久久香蕉激情| 午夜激情久久久久久久| 亚洲精品久久午夜乱码| 天天躁日日躁夜夜躁夜夜| 久久国产精品影院| 久久国产精品影院| 美女视频免费永久观看网站| 欧美变态另类bdsm刘玥| 极品教师在线免费播放| bbb黄色大片| 国产99久久九九免费精品| 精品国产超薄肉色丝袜足j| 国产免费av片在线观看野外av| 免费观看人在逋| 中文欧美无线码| 操出白浆在线播放| 久久国产精品人妻蜜桃| xxxhd国产人妻xxx| 成人特级黄色片久久久久久久 | 新久久久久国产一级毛片| 久久久久精品人妻al黑| 欧美精品一区二区大全| 午夜精品久久久久久毛片777| 久久精品国产99精品国产亚洲性色 | 五月天丁香电影| 女人被躁到高潮嗷嗷叫费观| 美女视频免费永久观看网站| 日本欧美视频一区| 90打野战视频偷拍视频| 欧美成狂野欧美在线观看| 久久精品国产99精品国产亚洲性色 | 黑丝袜美女国产一区| 桃红色精品国产亚洲av| 午夜激情久久久久久久| 午夜日韩欧美国产| 国产精品98久久久久久宅男小说| 日韩成人在线观看一区二区三区| 女人久久www免费人成看片| 成人手机av| 国产人伦9x9x在线观看| 在线 av 中文字幕| 亚洲精品在线观看二区| 亚洲精品粉嫩美女一区| 丝袜人妻中文字幕| 成人18禁在线播放| 水蜜桃什么品种好| 久久精品熟女亚洲av麻豆精品| 一夜夜www| 亚洲免费av在线视频| 久久精品国产综合久久久| 丝袜在线中文字幕| 国产精品二区激情视频| 久久精品亚洲av国产电影网| 色尼玛亚洲综合影院| 欧美乱码精品一区二区三区| 久久久久精品国产欧美久久久| 久久精品成人免费网站| 亚洲国产av新网站| 免费在线观看日本一区| 老司机深夜福利视频在线观看| 国产精品影院久久| 亚洲成av片中文字幕在线观看| 麻豆av在线久日| 日本av免费视频播放| 亚洲人成伊人成综合网2020| www.熟女人妻精品国产| 欧美性长视频在线观看| 热re99久久精品国产66热6| 国产高清激情床上av| 性色av乱码一区二区三区2| 免费不卡黄色视频| 久久精品国产综合久久久| 飞空精品影院首页| 亚洲久久久国产精品| 欧美日韩成人在线一区二区| 久久热在线av| 国产亚洲av高清不卡| 一个人免费看片子| 亚洲精品一二三| 日韩大码丰满熟妇| 中文字幕最新亚洲高清| av免费在线观看网站| 日韩 欧美 亚洲 中文字幕| 中亚洲国语对白在线视频| www.熟女人妻精品国产| 久久国产精品影院| 99国产极品粉嫩在线观看| 亚洲熟妇熟女久久| 久9热在线精品视频| 欧美久久黑人一区二区| 欧美日韩精品网址| 免费少妇av软件| 一个人免费在线观看的高清视频| 超色免费av| 黑人猛操日本美女一级片| 欧美 亚洲 国产 日韩一| 一边摸一边做爽爽视频免费| 他把我摸到了高潮在线观看 | 一级片'在线观看视频| 午夜91福利影院| 国产深夜福利视频在线观看| 久久久水蜜桃国产精品网| 久久性视频一级片| 国产精品免费大片| 超碰成人久久| 国产又爽黄色视频| 中文字幕人妻丝袜制服| 久久久久久亚洲精品国产蜜桃av| 热re99久久国产66热| 久久久久久人人人人人| 岛国毛片在线播放| 两人在一起打扑克的视频| 免费在线观看视频国产中文字幕亚洲| 涩涩av久久男人的天堂| 91av网站免费观看| 高清av免费在线| 国产精品久久久久成人av| videos熟女内射| 久久人妻熟女aⅴ| 首页视频小说图片口味搜索| 亚洲午夜理论影院| 乱人伦中国视频| 国产精品自产拍在线观看55亚洲 | 国产成人精品久久二区二区91| 无人区码免费观看不卡 | 午夜福利欧美成人| 91av网站免费观看| 视频区图区小说| 中文字幕精品免费在线观看视频| 欧美在线黄色| 亚洲av欧美aⅴ国产| 成人三级做爰电影| 亚洲五月色婷婷综合| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 精品国产国语对白av| 在线观看一区二区三区激情| 久久中文字幕人妻熟女| 国产亚洲欧美精品永久| 一二三四在线观看免费中文在| 天堂动漫精品| 成人国语在线视频| 激情在线观看视频在线高清 | 欧美精品高潮呻吟av久久| 日本a在线网址| www.熟女人妻精品国产| 亚洲精品国产区一区二| 欧美黑人欧美精品刺激| 国产国语露脸激情在线看| 无限看片的www在线观看| 国产高清激情床上av| 日本av免费视频播放| 在线观看免费高清a一片| 人人妻,人人澡人人爽秒播| 久久99一区二区三区| 久久久国产一区二区| 亚洲欧洲精品一区二区精品久久久| 黑人巨大精品欧美一区二区mp4| 久久久久久免费高清国产稀缺| 人妻一区二区av| 老司机午夜十八禁免费视频| 老熟妇乱子伦视频在线观看| 19禁男女啪啪无遮挡网站| 欧美成人免费av一区二区三区 | 男女高潮啪啪啪动态图| 亚洲国产成人一精品久久久| 欧美精品人与动牲交sv欧美| 免费一级毛片在线播放高清视频 | 女同久久另类99精品国产91| 日韩欧美三级三区| 高清毛片免费观看视频网站 | 午夜福利视频在线观看免费| 精品福利永久在线观看| 欧美大码av| 老汉色∧v一级毛片| 午夜福利,免费看| 欧美日韩一级在线毛片| 精品久久久精品久久久| 久久久欧美国产精品| 午夜福利视频精品| 欧美+亚洲+日韩+国产| 69av精品久久久久久 | 亚洲欧美色中文字幕在线| 极品教师在线免费播放| 国产97色在线日韩免费| 久热爱精品视频在线9| 国产欧美日韩精品亚洲av| 亚洲国产av影院在线观看| 国产老妇伦熟女老妇高清| 国产一卡二卡三卡精品| 免费在线观看黄色视频的| 欧美 亚洲 国产 日韩一| 人妻一区二区av| 国产成人欧美在线观看 | 无限看片的www在线观看| 亚洲,欧美精品.| 国产单亲对白刺激| 国产高清视频在线播放一区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品国产色婷婷电影| 大型黄色视频在线免费观看| 交换朋友夫妻互换小说| 久久青草综合色| 男女无遮挡免费网站观看| 别揉我奶头~嗯~啊~动态视频| 欧美 亚洲 国产 日韩一| 日韩 欧美 亚洲 中文字幕| 午夜精品久久久久久毛片777| 91精品国产国语对白视频| 免费看a级黄色片| 母亲3免费完整高清在线观看| 人人妻人人添人人爽欧美一区卜| 黄色 视频免费看| 久久人妻福利社区极品人妻图片| 久久久国产一区二区| 久久久久久久大尺度免费视频| 在线看a的网站| 久久人人爽av亚洲精品天堂| 亚洲国产毛片av蜜桃av| 如日韩欧美国产精品一区二区三区| 亚洲午夜理论影院| 宅男免费午夜| 日韩大片免费观看网站| 日本黄色视频三级网站网址 | 国产av又大| 黑人猛操日本美女一级片| 999久久久国产精品视频| 亚洲av欧美aⅴ国产| 午夜福利视频精品| 午夜福利一区二区在线看| 老司机午夜十八禁免费视频| 日韩中文字幕欧美一区二区| 91av网站免费观看| 脱女人内裤的视频| 精品人妻在线不人妻| 久久久精品国产亚洲av高清涩受| 亚洲欧美日韩另类电影网站| 极品教师在线免费播放| 精品少妇久久久久久888优播| 久久亚洲精品不卡| av片东京热男人的天堂| 国产成人欧美| 一级,二级,三级黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | av免费在线观看网站| 久久国产精品大桥未久av| 搡老乐熟女国产| 日韩成人在线观看一区二区三区| 欧美日韩一级在线毛片| 九色亚洲精品在线播放| 国产在线视频一区二区| 成人三级做爰电影| 18禁观看日本| 国产成人影院久久av| 午夜福利影视在线免费观看| 欧美午夜高清在线| 国产成人一区二区三区免费视频网站| 亚洲欧美精品综合一区二区三区| 黑丝袜美女国产一区| 久久人妻av系列| 亚洲第一av免费看| 99热国产这里只有精品6| 欧美成人午夜精品| 男女下面插进去视频免费观看| 午夜福利视频精品| 国产精品免费一区二区三区在线 | 亚洲国产精品一区二区三区在线| 黄色毛片三级朝国网站| 国产不卡av网站在线观看| 欧美一级毛片孕妇| 91字幕亚洲| 青草久久国产| 激情在线观看视频在线高清 | 免费高清在线观看日韩| 啦啦啦在线免费观看视频4| 国产高清激情床上av| av片东京热男人的天堂| 欧美亚洲日本最大视频资源| 日韩大码丰满熟妇| 9热在线视频观看99| 国产免费现黄频在线看| 一级毛片女人18水好多| 宅男免费午夜| 黑人猛操日本美女一级片| 国产在线观看jvid| 精品国产一区二区久久| 悠悠久久av| 久久午夜亚洲精品久久| 国产视频一区二区在线看| 亚洲av欧美aⅴ国产| av不卡在线播放| 亚洲精品粉嫩美女一区| 50天的宝宝边吃奶边哭怎么回事| 精品人妻在线不人妻| 日本撒尿小便嘘嘘汇集6| 精品国产亚洲在线| 电影成人av| av不卡在线播放| 国产一区二区激情短视频| 亚洲精品国产色婷婷电影| 欧美av亚洲av综合av国产av| 操出白浆在线播放| 另类精品久久| 人成视频在线观看免费观看| 日日爽夜夜爽网站| 美女高潮喷水抽搐中文字幕| 老司机福利观看| 在线观看一区二区三区激情| 色婷婷av一区二区三区视频| 欧美+亚洲+日韩+国产| 亚洲一码二码三码区别大吗| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产av影院在线观看| tube8黄色片| www.精华液| 91麻豆精品激情在线观看国产 | 亚洲天堂av无毛| 欧美+亚洲+日韩+国产| 99热网站在线观看| 午夜福利在线观看吧| 五月天丁香电影| 少妇裸体淫交视频免费看高清 | www.熟女人妻精品国产| 丁香欧美五月| 日韩中文字幕视频在线看片| 国产男靠女视频免费网站| 两人在一起打扑克的视频| 精品亚洲成a人片在线观看| 两人在一起打扑克的视频| 国产精品香港三级国产av潘金莲| 国产高清videossex| 91精品三级在线观看| 精品少妇内射三级| 人人妻人人澡人人爽人人夜夜| 国产成+人综合+亚洲专区| 大型av网站在线播放| 一级黄色大片毛片| 麻豆av在线久日| 欧美日韩黄片免| 欧美日韩视频精品一区| 国产精品久久久久久精品古装| 免费在线观看黄色视频的| 丝瓜视频免费看黄片| av又黄又爽大尺度在线免费看| 日韩精品免费视频一区二区三区| 亚洲精品国产色婷婷电影| 99精品欧美一区二区三区四区| 美女福利国产在线| 天天添夜夜摸| www日本在线高清视频| 亚洲七黄色美女视频| 激情在线观看视频在线高清 | 午夜福利在线观看吧| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片| 国产精品久久久久久精品古装| 别揉我奶头~嗯~啊~动态视频| 欧美日韩亚洲高清精品| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷成人精品国产| 欧美日韩一级在线毛片| 午夜免费成人在线视频| 免费不卡黄色视频| 在线亚洲精品国产二区图片欧美| 午夜精品久久久久久毛片777| 每晚都被弄得嗷嗷叫到高潮| 国产熟女午夜一区二区三区| av免费在线观看网站| 夜夜爽天天搞| 黄色成人免费大全| 亚洲av国产av综合av卡| 国产精品免费视频内射| 国产成人影院久久av| 老司机靠b影院| 亚洲一卡2卡3卡4卡5卡精品中文| 久久精品国产a三级三级三级| 一本色道久久久久久精品综合| av在线播放免费不卡| 熟女少妇亚洲综合色aaa.| 国产主播在线观看一区二区| www.精华液| 久久毛片免费看一区二区三区| 中文字幕最新亚洲高清| 高清视频免费观看一区二区| 中文字幕av电影在线播放| 国产xxxxx性猛交| 大片电影免费在线观看免费| 男女无遮挡免费网站观看| 熟女少妇亚洲综合色aaa.| 在线永久观看黄色视频| 国产精品 欧美亚洲| 国产av一区二区精品久久| 一级毛片电影观看| 亚洲精品一二三| 人人妻人人添人人爽欧美一区卜| 黄片小视频在线播放| 亚洲精品av麻豆狂野| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 午夜福利在线免费观看网站| 亚洲av片天天在线观看| 国产av精品麻豆| 欧美 亚洲 国产 日韩一| 人成视频在线观看免费观看| 日韩一区二区三区影片| 波多野结衣一区麻豆| 999久久久国产精品视频|