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

    A Case Study of the Initiation of Parallel Convective Lines Back-Building from the South Side of a Mei-yu Front over Complex Terrain

    2021-04-20 00:41:36QiweiWANGYiZHANGKefengZHUZheminTANandMingXUE
    Advances in Atmospheric Sciences 2021年5期

    Qiwei WANG, Yi ZHANG, Kefeng ZHU, Zhemin TAN*, and Ming XUE,2

    1Key Laboratory of Mesoscale Severe Weather and School of Atmospheric Sciences,Nanjing University, Nanjing 210023, China

    2Center for Analysis and Prediction of Storms and School of Meteorology,University of Oklahoma, Norman, OK 73072, USA

    ABSTRACT Parallel back-building convective lines are often observed extending to the southwest of some mesoscale convective systems (MCSs) embedded in the mei-yu front in China. The convective lines with echo training behavior can quickly develop into a stronger convective group of echoes, resulting in locally heavy rainfall within the mei-yu front rainband. The initiation mechanism of the back-building convective lines is still unclear and is studied based on high-resolution numerical simulation of a case that occurred during 27?28 June 2013. In the present case, the new convection along the convective lines was found to be forced by nonuniform interaction between the cold outflow associated with the mei-yu front MCSs and the warm southerly airflow on the south side of the mei-yu front, which both are modified by local terrain. The mei-yu front MCSs evolved from the western to the eastern side of a basin surrounded by several mesoscale mountains and induced cold outflow centered over the eastern part of the basin. The strong southwest airflow ahead of the mei-yu front passed the Nanling Mountains and impacted the cold outflow within the basin. The nonuniform interaction led to the first stage of parallel convective line formation, in which the low mountains along the boundary of the two airflows enhanced the heterogeneity of their interaction. Subsequently, the convective group quickly developed from the first stage convective lines resulted in apparent precipitation cooling that enhanced the cold outflow and made the cold outflow a sharp southward windshift. The enhanced cold outflow pushed the warm southerly airflow southward and impacted the mountains on the southeast side of the basin, where the roughly parallel mountain valleys or gaps play a controlling role in a second stage formation of parallel convective lines.

    Key words: back-building convective line, convective initiation, complex terrain, mei-yu front

    1. Introduction

    Mesoscale convective systems (MCSs) are often observed accompanying back-building convective lines which support their organization and maintenance for a long time and induce heavy precipitation (Schumacher and Johnson, 2005; Luo et al., 2013; Wang et al., 2014; Duffourg et al., 2016; Ding et al., 2018; Lee et al., 2018; Huang et al.,2019; Tsuguti et al., 2019). The back-building process consists of the periodic appearance of new cell upstream merging with the preexisting line of convection (Bluestein and Jain, 1985; Bluestein et al., 1987) or “echo training” of convective cells, and is primarily due to the interaction between the MCS cold outflow and the surrounding airflow.

    The characteristics of the MCSs and the associated back-building convective lines are different in different synoptic environments. The back-building MCSs east of the Rocky Mountains in the United States have a quasi-stratiform area of convection that produces a region of stratiform rain on the rear side of the MCS, and the back-building process is generally dependent on lifting provided by storm-generated cold pools (Schumacher and Johnson,2005). Over the northwestern Mediterranean, back-building quasi-stationary MCSs often form in an environment of moist and conditionally unstable air, which is carried by a southwesterly to southeasterly low-level marine flow(Ducrocq et al., 2013). The low-level cold pool under the MCS can lift the marine flow (Ducrocq et al., 2008) or modify the low-level circulation to enhance convergence forcing locally (Duffourg et al., 2016) to induce back-building processes. When this moist and conditionally unstable flow impinges on some of the mountains along the western Mediterranean, orographic lifting has been proposed as a mechanism for triggering heavy precipitation events involving back-building MCSs (Buzzi et al., 1998; Rotunno and Ferretti, 2001; Ducrocq et al., 2008).

    In east China, the MCSs associated with a quasi-stationary or slow-moving front during the mei-yu season (also known as the Changma season in Korea or the Bai-yu season in Japan) are essential systems for producing heavy rainfall events over the Yangtze-Huai Rivers basin (Ding, 1992;Chen et al., 1998; Zhang and Tan, 2009; Luo et al., 2013).Some MCSs embedded with quasi-linear back-building convection are often observed along the mei-yu front. Luo et al.(2013) investigated a back-building MCS along a mei-yu front during the midnight-to-morning hours and found two scales of the convective organization during the development of the MCS: one is the east- to northeastward “echo training” of convective cells along individual rainbands, and the other is the southeastward “band training” of the rainbands (the movement of the parallel rainbands) along the quasi-linear MCS. MCSs with a similar convective organization as in Luo et al. (2013) are also found in extreme rainfall events occurred over coastal South China during the pre-summer rainy season (e.g., Wang et al., 2014; Liu et al.,2018), while the mesoscale mountains near the coastal line sometimes can play a role in the convective initiation and organization.

    The back-building MCS in Luo et al. (2013) occurred in an area of flat terrain, while many mountains with a mean height of ~1000 m can be found over the Yangtze River-Huaihe River basin (Wang and Tan, 2006), which could have an effect on the back-building MCSs along the mei-yu front. Sun et al. (2015) found the mountains favor the formation and maintenance of the back-building pattern of a meiyu front MCS, but the role of the mountains was not revealed in detail. For the mei-yu front MCS over complex terrain, it is still not clear how the mountains modify the cold outflow under the MCS and how the mountains alter the back-building process of the convective lines associated with the MCS. Moreover, the mountains on the south side of the mei-yu front can modify the upstream airflow that approaches the mei-yu front, which can affect the evolution of the back-building MCS. This study will investigate a typical case involving back-building MCSs over complex terrain during the mei-yu season, with the primary aim being to focus on the role of the complex terrain in the formation of the back-building convective lines associated with the mei-yu front MCSs.

    The rest of this paper is organized as follows: Section 2 gives an overview of the case. The numerical model and experiment design for the simulation of the case are introduced in section 3. The results are discussed in detail in section 4. Finally, a summary is given in section 5, including a conceptual model of the role of the complex terrain in the formation of the back-building convective lines.

    2. Overview of the case

    This case took place over a basin area in Jiangxi Province of China, which is surrounded by the Wuyi,Yushan, Wugong, Jiuling and Mufu, Dabie, and Huangshan Mountains (Fig. 1a). The basin area is enlarged in Fig. 1b, and an important gap between the Yushan and Wuyi Mountains is denoted with an arrow. Our focus are the two stages of convective line formation, and the associated back-building processes that occurred primarily within two zones, which are denoted as two dashed rectangles in Fig. 1b, respectively.

    The NCEP 0.5°×0.5° reanalysis data (National Centers for Environmental Prediction, 2003) are used to portray the environmental conditions for the case. At 1800 UTC 27(0200 LST 28), June 2013 (Fig. 1c), the mei-yu front region characterized by weak temperature gradients and low wind speeds covered the basin area. The leading edge of the meiyu front, identified by the sharp contrast of the near-surface horizontal wind vectors at 10 m above ground level along the south boundary of the weak gradients of 2-m temperature, was located on the north side of the Nanling and Wuyi Mountains. To the west of the basin, there was a low pressure center embedded within the mei-yu front. The circulation of the low showed interactions with the southerly airflow ahead of the mei-yu front and also with the Jiuling and Mufu Mountains, where two cloud clusters labeled with A and B formed. The cloud clusters in Fig. 1 are outlined by the ?40°C isotherm of the infrared brightness temperature(TBB) retrieved from Multi-functional Transport Satellite(MTSAT).

    Over the next 6 hours, the cloud cluster B moved eastward and quickly developed within the basin area (Fig. 1d).The cloud cluster A also developed and propagated over the Jiuling and Mufu Mountains. The southerly airflow ahead of the mei-yu front became stronger at this time, and both of the two cloud clusters extended southwestward to the first back-building zone, showing more or less back-building features. The back-building features are hypothesized to be a result of interaction between the stronger southerly airflow and the cold outflow associated with the precipitation of the two cloud clusters.

    Fig. 1. Panels (a) and (b) show the terrain features. The terrain height is shaded in grayscale. A basin and surrounded mountains of Nanling, Wugong, Yushan, Wuyi, Jiuling, Mufu, Dabie, and Huangshan are denoted in (a). The basin area outlined with a solid box in (a) is enlarged in (b) with two zones of back-building and the gap between Yushan and Wuyi Mountains denoted. Panels (c?f) show the mean sea level pressure (blue contours, hPa), 2-m temperature(red contours, °C), 10-m horizontal winds (vectors, m s?1), and the infrared brightness temperature (TBB) at –40°C(thick green contours) at (c) 1800 UTC 27 June 2013, (d) 0000, (e) 0600, and (f) 1200 UTC 28 June 2013. The fields of pressure, temperature, and wind are plotted from the NCEP 0.25°×0.25° reanalysis data. The TBB data are obtained from the satellite data archive for research and education at http://weather.is.kochi-u.ac.jp/sat/GAME/. The evolution of two focused cloud clusters (labeled with “A” and “B”) in (c) is shown in Fig. 2. The bold dashed line denotes the south leading edge of the mei-yu front. The first zone and the second zone of back-building are also shown in (d?e) and (e?f), respectively. See the text for the details.

    Six hours later (Fig. 1e), cloud cluster B weakened and moved out of the basin area, and cloud cluster A also weakened and separated into cloud cluster segments. To the southwest of cloud cluster A, a new cloud cluster C formed,which also showed a back-building character on its southwest side within the first back-building zone. The cloud clusters C and A exhibited linear shapes roughly parallel to each other, and their back-building development appeared to push the leading edge of the mei-yu front southward.

    After another 6 hours (Fig. 1f), cloud clusters C and A also moved out of the basin area and dissipated. Meantime,a new cloud cluster D covering the second back-building zone, which originated roughly from where C quickly grew,and a new cluster E east of D developed rapidly along the Yushan and Wuyi Mountains, and further pushed the leading edge of the mei-yu front southward.

    The detailed evolution of the cloud clusters mentioned above is shown in Fig. 2 in hourly maps of TBB. Note the line features of the cloud clusters A and C in Fig. 2d, which can be identified from the lower values of TBB and the wave pattern of the south sides of the two cloud clusters.This is the first stage (stage I) of back-building starting from the first back-building zone shown in Fig. 1b. The development of cloud cluster D also shows line features of backbuilding (Fig. 2j), and the back-building development of the cloud cluster E east to D is much clearer (Figs. 2h-j). The second stage (stage II) of back-building starts from the second back-building zone shown in Fig. 1b and is also associated with the interaction between the cold outflow in the basin area and the southerly airflow ahead of the mei-yu front. At the same time, the mountains should play an important role as the leading edge of the mei-yu front moved southward.

    3. Numerical model and experiment design

    Version 5 of the Advanced Regional Prediction System (ARPS, Xue et al., 2000, 2001, 2003), a non-hydrostatic atmospheric model, is used to perform the simulations of the case. The model domain has 651 and 579 horizontal grid points in X and Y directions with a 1-km grid spacing, which is nested in a 3-km X–Y grid of 483 × 403 horizontal grid points (Fig. 3). A generalized terrain-following coordinate is used in vertical, and 53 stretched levels are defined with the grid spacing increasing from about 20 m near the ground to about 800 m near the model top at about 20 km. The model terrain and land surface characteristics on the 3- and 1-km grids are derived from the 30-arc seconds (approximately 1 kilometer) global elevation and land-use data set from the U.S. Geological Survey.

    The model is initiated at 1200 UTC 27 June 2013 and is integrated for 36 hours until 0000 UTC 29 June 2013. For the 3-km grid, the lateral boundary conditions (LBCs) are generated from GFS global 0.5-degree analysis fields at 6-hour intervals obtained from the Historical Unidata Internet Data Distribution (IDD) Gridded Model Data (National Centers for Environmental Prediction, 2003). The initial condition for the 3-km grid is created using the ARPS three-dimensional variational (3DVAR) data analysis scheme (Gao et al., 2004) based on the GFS global 0.5-degree analysis fields and the conventional radiosonde and surface observations shown in Fig. 3. The nested 1-km grid gets its initial condition from the interpolation of the 3-km analysis and its LBCs from the 3-km forecasts at 10-min intervals. The ARPS is used in its full physics mode, as in Wang and Xue(2012).

    In this case, the cold outflow is generally confined within the basin surrounded by the Mufu-Jiuling and Huangshan Mountains on the north side, and several low mountains are scattered among the first back-building zone where stage I convective lines formed. It is not clear whether the low mountains play a role in the formation of stage I convective lines. Neither is it clear whether the Mufu-Jiuling and Huangshan Mountains play a role in altering the strength of the cold outflow within the basin and, thus, the back-building processes associated with the cold outflow interacting with the southerly flow along the leading edge of the meiyu front during the two stages of convective line formation.To evaluate the two points in question, sensitivity experiments with the mountains parallel to the leading edge of stage I convective lines removed (NoMidMts) and with the Mufu-Jiuling and Huangshan Mountains on the north side of the basin removed (NoNorthMts) were designed, respectively (Table 1).

    4. Results

    4.1. The simulated convective lines and cold outflow

    Figure 4 shows the observed (left, obtained from radar observations at Nanchang, Ji’an, and Shangrao as shown in Fig. 3) and the simulated (right) composite radar reflectivity (RFC) at 0300 and 0900 UTC 28 June 2013. The lines of the cloud cluster in Figs. 2d and 2j primarily correspond to the lines of RFC shown as dashed lines in Figs. 4a and 4c,respectively. The line features of the RFC at the two stages are generally reproduced in the simulation (dashed lines in Figs. 4b and 4d), except the simulated RFC is slightly greater, and the position of the convective lines is slightly to the south of the observed convective lines. The simulation also reproduces the broad stratiform region that follows the convective lines at the two stages (Figs. 4a?b and 4c?d).The convective line corresponding to the cloud cluster E in Fig. 2j was initiated later in the simulation and did not separate from the other convective lines (Fig. 4d), as observed(Fig. 4c), which will be discussed next.

    Fig. 2. Evolution of TBB (shaded) in the solid box region in Fig. 1a at (a) 0000, (b)0100, (c) 0200, (d) 0300, (e) 0400, (f) 0500, (g) 0600, (h) 0700, (i) 0800, and (j) 0900 UTC 28 June 2013. The focused cloud clusters are labeled as A, B, C, D, and E. The black dashed lines denote the positions of the cloud clusters. The red dashed lines in(d) and (j) show the line features of the cloud clusters. The TBB data are obtained from the satellite data archive for research and education at http://weather.is.kochiu.ac.jp/sat/GAME/.

    Fig. 3. The 3-km model domain nested with a 1-km subdomain (solid box). The small triangles and circles mark the stations from conventional radiosonde and surface networks,respectively. The small squares mark radar locations of Nanchang (NC), Ji’an (JA), and Shangrao (SR), with large circles indicating their maximum radar ranges. The dashed box denotes the horizontal plotting region in later figures. The terrain elevation is shaded in grayscale with mountainous regions denoted as in Fig. 1a.

    Table 1. Experiment design.

    The two stages of back-building convective line formation are associated with the interaction between the cold outflow in the basin and the southerly airflow ahead of the meiyu front, as discussed in section 2. The cold outflow linking with precipitation in the basin is shown in Fig. 5. The precipitation data used in Figs. 5a and 5b is the hourly, 0.1° × 0.1°,gauge-satellite merged precipitation product, which can be obtained from the National Meteorological Information Center of China. (available at http://data.cma.cn/data/detail/data-Code/SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRI D_0.10/keywords/CMORPH.html). The station temperature in Figs. 5c and 5d is taken from the observations of the automatic weather stations (AWS).

    The centers of accumulated precipitation from 0000 UTC 27 to 0300 UTC 28 June 2013 are mainly over the Jiuling Mountains and the western part of the basin between the Wuyi and Huangshan Mountains (Fig. 5a), roughly corresponding to the centers of cold outflow when stage I convective lines formed (Fig. 5c). The development of stage I convective lines induced heavy precipitation over the next 6 hours, and the accumulated precipitation by 0900 UTC 28 June 2013 becomes larger, spreading southward and extending to the southwestern part of the basin between the Wugong and Yushan Mountains (Fig. 5b). The cold outflow centers correlate less with the centers of the accumulated precipitation (Figs. 5b and 5d) due to the cold outflow that responded more directly to the recent 6-hour precipitation. The cold outflow advanced towards the Yushan and Wuyi Mountains and pushed the leading edge of the mei-yu front (stronger gradient of temperature) southeastward over the next 6 hours (Figs. 5c and 5d). It is conjectured that the Yushan and Wuyi Mountains may have modified the interaction between the cold outflow and warm southerly airflow along the leading edge of the mei-yu front, which should connect with the formation of stage II back-building convective lines.

    The simulated cold outflow in the basin is weaker than shown in the observations (Figs. 5c?d and 5e?f), while the cold outflow centers over the western part of the basin between the Wuyi and Huangshan Mountains (Fig. 5e) and approaching the Yushan and Wuyi Mountains (Fig. 5f) are generally captured in the simulation. The simulated weaker cold outflow is partially associated with initial errors and model uncertainties, which starts at 1200 UTC 27 June 2013 and failed to produce the observed long-time accumulated precipitation, as in Figs. 5a, b. On the other hand, were the simulation to have begun earlier, the prediction of the convective line formation would not be as good due to the uncertainty associated with the initiation and evolution of the convection in the model (not shown).

    Fig. 4. Observed (left) and simulated (right) composite radar reflectivity at (a, b) 0300 and (c, d) 0900 UTC 28 June 2013. The dashed lines denote the line features of the composite radar reflectivity.

    The warm southerly airflow near the surface is also poorly simulated in terms of temperature (Figs. 5c?d and 5e?f). At 0300 UTC 28 June 2013, the simulated boundary of the warm southerly airflow is diffuse near the surface(Fig. 5e). This may be partially associated with the imperfect representation of the mountain and valley circulations to the lee side of the Nanling Mountains due to shorter integration of the model. By 0900 UTC, the simulated near-surface boundary of the southerly airflow becomes clearer and more consistent with that of the observations (Figs. 5d and 5f). While the boundary differs locally near the gap between the Yushan and Wuyi Mountains, it is the warm southerly airflow in the observations (Fig. 5d) or the cold outflow in the simulation (Fig. 5f) that intrudes into the gap. The cold outflow intrusion into the gap in the simulation results in a later formation of the convective line than that in the observations (the right convective line in Figs. 4c and 4d).

    Note that the orientation of the stratiform region in Fig. 4 is roughly west-east at 0300 UTC, but changes to southwestnortheastward by 0900 UTC 28 June 2013. Figure 2 suggests that the stratiform region partially comes from the decaying of the former convective groups along the mei-yu front.When new convective regions generated on the south side of the stratiform region, processes associated with the detrainment of air from the convective regions and the transition of the weakening convective regions to stratiform regions should have played an important role as discussed in Luo et al. (2010). Figure 6 shows the soundings at Nanchang within the stratiform region at two times. Mid-level westerly winds are seen between 700 and 500 hPa, and the mean wind direction between the two levels primarily indicates the orientation of the stratiform region at the two times,respectively. Generally, the orientation of the stratiform region is directly linked with mid-level flow associated with the front. The motion of the convective lines shows a similar orientation to that of the stratiform region due to the mid-level flow, while the orientation of the convective lines or cell motion closely follows the direction of the southwest flow that impinges on the mei-yu front.

    Fig. 5. Accumulated precipitation from 0000 UTC 27 June to (a) 0300 and (b) 0900 UTC 28 June 2013, and observed near-surface temperature at (c) 0300 and (d) 0900 UTC 28 June 2013. The plots of simulated near-surface temperature at 0300 and 0900 UTC 28 June 2013 are shown in (e) and (f), respectively. The red star in (a) denotes the location of Nanchang (NC), where the soundings are extracted from model forecasts and are shown in Fig. 6. The observed precipitation data are from the gauge-satellite merged precipitation data obtained from the National Meteorological Information Center of China (available at http://data.cma.cn/data/detail/dataCode/SEVP_CLI_CHN_MERGE_CMP_PRE_HOUR_GRID_0.10/keywords/CMORPH.html). The observed temperature data were from the observations of the automatic weather stations (AWS) shown as dots in (c).

    Fig. 6. Skew-T plots of soundings extracted from the model simulation at th e red star in Fig. 5a, at 0300 (red), and 0900(blue) UTC 28 June 2013.

    Below 900 hPa, southeast winds dominate (Fig. 6), indicating that the cold air flows from the cold outflow center to the sounding location at Nanchang. As the leading edge of the front advanced southward between 0300 and 0900 UTC 28 June 2013, the main body of the stratiform region gradually deviated from Nanchang, leaving drier and slightly warmer air in the mid-and-lower levels at Nanchang, but the near-surface cold airflow can still reach Nanchang, making the near-surface air at 0900 UTC colder than that at 0300 UTC (Fig. 6).

    Overall, the convective lines, the cold outflow within the basin, and the warm southerly airflow ahead of the meiyu front are reasonably simulated with acceptable correspondence to the observations, despite differences in strength and position details. This establishes the credibility of the model simulation to be further analyzed with a certain degree of confidence in the next sections.

    4.2. Formation of convective lines during stage I

    Low-level horizontal cross-sections at heights of 250 m and 500 m mean sea level (MSL) are plotted in Fig. 7 during the formation of stage I convective lines. The displays clearly show the interaction between the cold outflow and the warm southerly airflow and their modification by the local mountains. The convective lines along A1B1-A5B5 in Figs. 7c?f correspond to those shown as dashed lines in Fig. 4b.

    Within 30 minutes prior to the appearance of the line features of the convective group (Figs. 7a?d), the cold outflow centers at 250 m and 500 m MSL are located over the western parts of the basin, and between the Wuyi and Huangshan Mountains, similar to the cold outflow centers near the surface (Fig. 5e). The boundary between the cold outflow and the warm southerly airflow at 250 m MSL is quite consistent with that near the surface (Figs. 5e and 7c) and is in a quasi-steady state (Figs. 7a and 7c). Within the basin, the airflow at 250 m MSL is generally colder than that near the surface because of the temperature lapse rate and less accumulated precipitation and cooling near the surface due to a shorter integration of the model.

    At 0230 UTC (Fig. 7a), the warm southerly airflow passes the Nanling, Wugong, and Yushan Mountains and converges with the cold outflow roughly along a line from the Jiuling Mountains to the north of the Yushan Mountains(AB in Figs. 7c?f). At 500 m MSL (Fig. 7b), the airflow pattern is similar, while the southerly airflow is stronger than that at 250 m MSL, resulting in clearer convergence forcing where it meets the cold outflow. The strongest convergence forcing is found downstream of the Wugong Mountains at both of the two levels due to strong lee side convergence (Figs. 7a?b).

    The airflows at the two levels maintain a quasi-steady pattern over the next 30 min (Figs. 7c?d), except that the southerly airflow becomes slightly warmer and stronger. The stronger southerly airflow strengthens the convergence forcing downstream of the Wugong Mountains, supporting the maintenance of convection along A1B1 on the warm air side of the convergence zone. The convection along A1B1 on the cold air side displays some features of back-building(Figs. 7c and 7e). On the warm air side, two antecedent convective elements advected toward and tended to merge with the cold air side convective group where back-building features appear. This causes the convective group to take on a linear appearance along A2B2 and A3B3. Where the convective lines A4B4 and A5B5 extend southwestward from the convective group, back-building features show most clearly.

    It is easy to understand that the western part of the convective group shifts to the south of the cold-warm airflow boundary due to the convection forming on the warm air side since back-building is generally associated with coldwarm airflow interaction along its boundary. In this case,the interaction between the cold and warm airflow is nonuniform. One possible reason for this is that this interaction was modified by the local mountains, as shown in the vertical cross-section roughly along their boundary (Fig. 8a).Also, the interaction was likely affected by the convection on the warm air side. The detailed interaction and the associated back-building processes are shown in vertical cross-sections through the convective lines A1B1-A5B5 in Figs. 8b?f, respectively.

    Figure 8a shows the vertical cross-section along AB in Fig. 7c, five updraft branches can be identified where the convective lines A1B1?A5B5 intersect with AB. The warm(cold) air tends to be deeper on the left (right) side of each of the terrain profiles along AB, indicating modification of the cold-warm airflow interaction by the local terrain.Along each of the five convective lines (Figs. 8b?f), successive updrafts are found to be tilting over the cold air side,which is associated with the general structure of the mei-yu front. New convective updrafts feeding the mei-yu front updraft continuously generate near the cold-warm airflow boundary at low levels, clearly demonstrating a back-building structure.

    Fig. 7. Temperature (shaded, °C), wind vectors, and convergence (blue contours, with an interval of 2×10–3 s–1) at(left) 250 m and (right) 500 m MSL, with the times at (a, b) 0230 and (c, d) 0300 UTC 28 June 2013. The red contours show the 45 dBZ composite radar reflectivity at the corresponding times. The lines in (c) and (d) indicate the position of the vertical cross-sections shown in Fig. 8. The dashed box regions in (c) and (d) are enlarged in (e)and (f), respectively.

    Upstream of the slantwise updrafts along the convective line A1B1 (Fig. 8b), convection forced by the convergence on the lee side of the Wugong Mountains (Figs. 7c?f)produces an apparent downdraft, which enhances the lowlevel convergence that contributes to the formation of the back-building new updraft. Along the convective line A2B2(Fig. 8c), upstream convection is induced by low mountains, which is weak and has no obvious effects on the downstream back-building new updraft. In contrast, the upstream convection along the convective line A3B3 (Fig. 8d) is strongly amplified by higher mountains. The amplified convection can produce a downdraft that strengthens the lowlevel convergence supporting the formation of the back-building new updraft as that along the convective line A1B1. The back-building features along the convective lines A4B4(Fig. 8e) and A5B5 (Fig. 8f) are clearer since there is no upstream convection. For some of the convective lines encountering the mountains near the cold-warm airflow boundary, the local mountains play a role in the formation of the back-building new updraft, e.g., through upstream lifting (e.g., Figs. 8d, f) or lee side downslope wind convergence (Fig. 8f).

    Fig. 8. Vertical cross sections through the lines (a) AB, (b) A1B1, (c) A2B2, (d) A3B3, (e) A4B4, and (f) A5B5 in Figs. 7c?f. Shown are temperature (shaded, °C), wind vectors, vertical velocity (blue contours, with contour levels at 0.5, 1, 2, and 5 m s?1), and 45 dBZ radar reflectivity (red contours).

    4.3. Formation of convective lines during stage II

    The processes associated with the formation of the convective lines during stage II are shown in Fig. 9. The development of the convective lines of stage I results in apparent precipitation-forced latent cooling, and the centers of the cold outflow have shifted to the southeastern part of the basin along the Yushan and Wuyi Mountains by 0830 UTC 28 June 2013 (Figs. 9a and 9b). Note the gap between the Yushan and Wuyi Mountains and the valleys among the mountain branches on the northwest side of the Yushan Mountains. The cold outflow first runs into the gap and the valleys, and strong convergence between the cold outflow and the warm southerly airflow is found along the leading edge,which is more clearly shown at 500 m MSL (Fig. 9b). The convergence boundary is followed by several back-building convective lines, indicating that the forcing for back-building comes from the cold-warm airflow interaction associated with terrain modification at low levels.

    Fig. 9. The same as in Fig. 7, but for times at (a, b) 0830 and (c, d) 0900 UTC 28 June 2013. The lines in (c) and (d)indicate the position of the vertical cross-sections shown in Fig. 10. The dashed box regions in (c) and (d) are enlarged in (e) and (f), respectively.

    The airflows described above maintain a similar interaction pattern for the following 30 min (Figs. 9c?d). During this period, the cold outflow intrusion is quick and strong and results in a stage of clearer back-building at 0900 UTC(Figs. 9c and 9d). The convective group extends southwestward roughly along line CD over the Yushan Mountains.Along the line, the convergence boundary shows a clear wave pattern (Fig. 9f), suggesting modification by the local terrain. Five convective lines are primarily identified. The convective line C1D1 is along the Yushan Mountains and the lee side of the Nanling Mountains. The convective lines C2D2?C4D4 are roughly along the major valleys on the northwest side of the Yushan Mountains, and the convective line C5D5 primarily follows the gap between the Yushan and Wuyi Mountains. Vertical cross-sections along these lines shown in Fig. 10 illustrate the processes resulting in the forcing that supports the back-building.

    Warm air alternating with cold air along the line CD show nonuniform interaction between the warm southerly airflow and the cold outflow (Fig. 10a), and the associated updrafts can be found roughly where the horizontal gradient of temperature is higher. The updrafts corresponding to the selected convective lines C2D2?C5D5 are near where they intersect with the line CD (Figs. 9c and 9d). For the convective line C1D1 that extends to the lee side of the Nanling Mountains (Fig. 10b), the airflow interaction at low levels and the associated convergence forcing are not strong due to the mountain-reduced southwest airflow on the lee side. Instead of being surface-based, the new updraft is due to the interaction between the mid-level (3?5 km MSL) stronger southerly airflow and the leading convection (Fig. 10b).The back-building new updrafts for convective lines C2D2?C5D5 are all surface-based (Figs. 10c?f), but the sources of forcing are different in detail depending on the cold-warm airflow boundary related to the topography.When the cold outflow has not reached the mountain top(Fig. 10d), the cold-warm airflow convergence provides the forcing for back-building on the near-basin side of the mountain, and the existing convection gains support from the lifting of the cold outflow by the mountain. As the cold outflow just surmounts the mountain top (Fig. 10c), the leading updraft is associated with the convergence between the warm air and a thin near-surface layer of cold downslope outflow and the lifting of the warm southerly airflow by the mountain above. When the leading convection passes the top of the mountain (Fig. 10e), the mountain downslope can enhance the cold outflow associated with the downdraft of the convection, which will strengthen the cold-warm airflow interaction, and thus the leading updraft. Along the gap between the Yushan and Wuyi Mountains (Fig. 10f), the processes are similar to that along convective line C4D4(Fig. 10e), except the topographic effects can be neglected since the terrain along the gap is much lower.

    4.4. The role of the mountains within and surrounding the basin

    For the formation of the convective lines during stage I,the low mountains along the leading edge of the cold outflow are within the basin. They are discontinuous (Figs. 7c and 8a), and the importance of these low mountains is unclear since some studies show that back-building convective lines can form over flat terrain (e.g., Luo et al., 2013).For the convective line formation of stage II, the mountains ahead of the cold outflow are on the southeast side of the basin and are continuous (Figs. 9c and 10a). The mountains play an important role, as was discussed in section 4.3. The experiment NoMidMts investigates the role of the low mountains within the basin in the formation of stage I convective lines, and the experiment NoNorthMts will examine the role of the mountains on the north side of the basin that probably changes the strength of cold outflow and subsequently the back-building processes.

    As the control experiment (CNTL), both of the two experiments are initiated at 1200 UTC 27 June 2013 and are integrated until 0000 UTC 29 June 2013. The modified terrain for each of the two experiments is smoothed to avoid integration instability associated with terrain removal that could result in a very sharp terrain profile. The smoothed terrain for the two experiments maintain the primary characteristics of the terrain in the experiment CNTL where they are not modified (Figs. 11 and 12), whereas the modified terrain should alter the processes of the convective line formation associated with different terrain conditions.

    In the experiment NoMidMts (Fig. 11), the convective group shows line features. Still, the convective line number decreases since the interaction between the cold outflow and the warm southerly airflow becomes less ununiform without being disturbed by the low mountains (Figs. 11a and 11b).In the absence of the Wugong Mountains and the enhanced convergence forcing on its lee side, the convective group is primarily along the cold-warm airflow boundary instead of shifting its western part to the warm side of the boundary as in the experiment CNTL (Figs. 11a and 7c). Some parts of the Yushan and Wuyi Mountains are also removed in the experiment NoMidMts. This makes the cold outflow spreads toward the Yushan and Wuyi Mountains more quickly than that in the experiment CNTL (Figs. 11c?d and 9c?d).

    The experiment NoNorthMts shows that the cold outflow over the original basin area becomes broader (Fig. 12),merging with the cold air that should be between the north side of the mei-yu front and the Mufu-Jiuling and Huangshan Mountains in the experiment CNTL (Fig. 1b). The merger more or less enhances the cold outflow over the basin area, whose low-level leading edge by 0300 UTC 28 June 2013 is roughly at y = 300 km (Figs. 12a?b) that is south to the leading edge of the cold outflow in the experiment CNTL (Figs. 7c?f). The convective group in the experiment NoNorthMts forms and develops further southward and later results in a quicker advance of the cold outflow into the Yushan and Wuyi Mountains than that in the experiment CNTL (Figs. 12c?d and 9c?d).

    Fig. 10. The same as in Fig. 8, but for vertical cross sections through the lines (a) CD, (b) C1D1, (c) C2D2, (d)C3D3, (e) C4D4, and (f) C5D5 in Figs. 9c?f.

    Fig. 11. Fields shown are the same as in Fig. 7 but for experiment with the mountains within the basin removed(NoMidMts) at times of (a, b) 0300, and (c, d) 0900 UTC 28 June 2013.

    To illustrate the terrain effects on the low-level forcing that contribute to the formation of the convective lines during stage I, the time evolution of key relevant fields along AB in Fig. 7c at 250 m MSL is compared between the CNTL case and the NoMidMts and NoNorthMts experiments in Fig. 13. In the CNTL case (Figs. 13a and 13b),more low-level convergence bands (slash-filled blue contours in Fig. 13a) coupling with lifting forcing (slash-filled blue contours in Fig. 13b) form due to ununiform coldwarm airflow interaction. The terrain around x = 100 km induces warm airflow convergence on the left side and cold outflow convergence on the right side. This results in convergence forcing supporting the development and maintenance of two of the convective lines (dash-filled red contours in Figs. 13a and 13b) later. The right shifting of the convective lines with time clearly shows an eastward “band training” organization of the convective lines, as summarized in Luo et al. (2013). In the NoMidMts experiment (Figs. 13c and 13d), bands of convergence forcing cannot form around x = 100 km where the terrain is removed, and the adjacent convective lines are more intermittent instead of continuous.The cold outflow becomes stronger in the NoNorthMts experiment (Figs. 13e and 13f), and it is clear that the terrain along AB in Fig. 7c still helps to maintain warm airflow convergence on its left side and cold outflow convergence lifting on its right side. This confirms that the low mountains along the cold-warm airflow boundary increase the non-uniformity of the cold-warm airflow interaction by inducing long-lasting low-level convergence forcing that favors the formation of back-building convective lines.

    The time evolution of the same fields along CD in Fig. 9d at 500 m MSL is shown in Fig. 14 to illustrate the low-level forcing associated with terrain effects that contribute to the formation of the convective lines during stage II. In the CNTL case (Figs. 14a and 14b), the convergence lifting forcing starts where the cold outflow meets the warm southerly airflow. The wavy boundary of cold-warm airflow interaction is primarily controlled by the mountain tops and valleys along CD (Fig. 10a). Several convective lines originate from some of the convergence forcing, and their maintenance along CD indicates continuous back-building that supports the formation of some of the convective lines shown in Figs. 9c?d. Note not all the convective lines in Figs. 9c?d have corresponding RFC field in Figs. 14a?b because the line CD is roughly through the cold-warm airflow boundary instead of crossing the leading parts of all the convective lines. As with the organization of the convective lines during stage I, the northeastward “band training” along CD is clear during this stage. Similar features exist in the NoMid-Mts (Figs. 14c and 14d) and NoNorthMts (Figs. 14e and 14f) experiments. However, the cold-warm airflow interaction near the line CD occurs about 1 hour earlier due to quicker development of the stronger cold outflow in the two experiments, as discussed above. The terrain plays a controlling role in the convective line formation during stage II.

    Fig. 12. Fields shown are the same as in Fig. 7 but for experiment with the mountains on the north side of the basin removed (NoNorthMts) at times of (a, b) 0300, and (c, d) 0900 UTC 28 June 2013.

    5. Summary and conclusions

    The formation of parallel convective lines along the south side of the mei-yu front during 28 June 2013 in China is studied through numerical simulations at 1 km horizontal grid spacing using a mesoscale numerical model. Two stages of convective line formation observed with back-building features are reasonably simulated.

    The two stages of convective line formation are associated with the interaction between the warm southerly airflow ahead of the mei-yu front and the cold outflow induced by the mei-yu front mesoscale convective systems (MCSs)evolving over a basin surrounded by the Wuyi, Yushan,Wugong, Jiuling, Mufu, Dabie, and Huangshan Mountains.The mountains within and on the southeast side of the basin play important roles during the two stages of convective line formation, respectively, which is summarized in a conceptual model, as illustrated in Fig. 15.

    Fig. 13. Time evolution of (a) temperature (shaded) and convergence (blue contours, > 0.5×10?3 s?1 filled with slashes), or (b) cross-line horizontal wind speed (shaded) and vertical velocity (blue contours, > 0.1 m s?1 filled with slashes) during stage I along line AB in Fig. 7c at 0.25 km MSL. The same as (a) and (b), respectively, (c) and (d)are for experiment NoMidMts, and (e) and (f) are for experiment NoNorthMts. The 45 dBZ composite radar reflectivity shown as red contours filled with dashed lines outline the convective lines.

    During stage I (Fig. 15a), the warm southerly airflow passes the Nanling, Wugong, and Yushan Mountains and meets the cold outflow within the basin. A convection region initiates from the cold-warm air boundary (or the leading edge of the mei-yu front). It is followed by a main stratiform region within a broad cloud sector covering the basin.The cold outflow within the basin is primarily due to precipitation-forced cooling associated with MCSs that propagated along the mei-yu front from the Mufu-Jiuling Mountains while it was prevented by the Wuyi and Huangshan Mountains on the east side of the basin earlier. The cold outflow center is between the Wuyi and Huangshan Mountains,where the accumulated precipitation is larger. The low-level cold outflow spreads from the cold outflow center and interacts with the warm southerly airflow nonuniformly along the leading edge of the mei-yu front. The low mountains along the leading edge of the mei-yu front can increase the non-uniformity of the interaction. Back-building processes continuously occur where the interaction is stronger along the quasi-steady leading edge of the mei-yu front and supports the formation and maintenance of the convective lines.Also, low-level convergence forcing due to interaction between the mountains and the warm southerly airflow induces several convective pieces ahead of the mei-yu front,especially the convergence forcing to the lee side of the Wugong Mountains. The convective elements are found in line with and approaching some of the convective lines. The associated downdrafts were shown to have played a role in enhancing the back-building processes near the leading edge of the mei-yu front. For each of the convective lines,the cell motion is along the convective line. In contrast, the motion of the convective lines is primarily consistent with the mid-level west-eastward flow associated with the front.

    Fig. 14. The same as Fig. 13, but during stage II, along line CD in Figs. 9c at 0.5 km MSL. The initiation of convective lines in (c?d) and (e?f) is about 1 hour earlier than that in (a?b).

    The development and eastward propagation of the convective lines formed during stage I results in new precipitation-forced cooling that strengthens the cold outflow. The cold outflow center becomes broad and extends southwestward along the southeast side of the basin, resulting in a sharp southwest shift of the convection region, the stratiform region, and the cloud sector (Fig. 15b). The strong cold outflow intrudes into the valleys on the northwest side of the Yushan Mountains and the gap between the Yushan and Wuyi Mountains and pushes the warm airflow southward. Convective lines form along the main valleys and the gap, with the convective line along the gap being stronger since the gap is broader than the valleys, and the intrusion of the cold outflow into the gap is more intense. The gap plays a role in accelerating the cold outflow and thus to enhance its interaction with the warm southerly airflow. The valleys can induce additional forcing through confluence and uplifting or downslope flow to enhance the cold-warm airflow interaction. These result in a controlling role of the terrain in the formation of the convective lines during stage II.During this stage, the cell motion is also along each of the convective lines, and the motion of the convective lines is still primarily consistent with the mid-level flow associated with the front, which has changed to be southwest-northeastward.

    Fig. 15. The conceptual model of convective line formation during (a) stage I and (b) stage II. The warm southerly airflow boundary or the south leading edge of the mei-yu front, the cold outflow center, and the cold outflow boundary are shown as red bold, blue, and blue-dashed lines, respectively. The main convection region, the main stratiform region, and the cloud boundary are outlined with red-dashed, orange, and green-dashed lines, respectively.Red ellipses embedded with small ovals indicate the convective lines, with a wider one denoting a stronger convective line. Small red ellipses represent the convection ahead of the mei-yu front. The motion of the convective cells along the convective line and the motion of the convective lines primarily following the mid-level flow are denoted with bold arrows. Terrain height is shaded, and the main mountains are labeled in (a). Important low-level wind vectors are shown, and the black ellipse in (a) indicate a sustained low-level convergence on the lee side of the Wugong Mountains.

    This conceptual model contains a broad stratiform region and is different from the schematic diagram of the quasi-linear-shaped MCS in Luo et al. (2013), in which the stratiform rainfall is insignificant. The conceptual model appears to be similar to that of the MCS climatology with training line and adjoining stratiform in North America (Schumacher and Johnson, 2005) in terms of the stratiform region, but differs in terms of the causes for the formation of the convective lines and the motion of the echo training.The echo training along the convective lines is similar to that of the back-building MCS in Schumacher and Johnson(2005) and Luo et al. (2013). More importantly, the mountains play significant roles in the conceptual model. Sensitivity experiments confirm the modification role of the low mountains within the basin during stage I and the controlling role of the Yushan and Wuyi Mountains on the southeast side of the basin during stage II of the convective line formation. Knowledge of the effects of the mountains on the convective line formation is helpful for the understanding and prediction of heavy precipitation events over the basin region during the mei-yu season in China.

    This work was primarily supported by the National Key R&D Program of China (Grant No.2017YFC1501603), the Major Research Plan of the National Natural Science Foundation of China (Grant No. 91937301), and the National Natural Science Foundation of China (Grant Nos.41775054, 41375061 and 40705019). The NCEP reanalysis data used to reproduce the weather background and to initialize the simulations were downloaded from the NCAR Data Archive. The infrared brightness temperature data used to show the cloud clusters were obtained from the satellite data archive for research and education at http://weather.is.kochi-u.ac.jp/sat/GAME/. The observed precipitation data were obtained from the National Meteorological Information Center of China (available at http://data.cma.cn/data/detail/dataCode/SEVP_CLI_CHN_MERGE_CMP_P RE_HOUR_GRID_0.10/keywords/CMORPH.html). The observed station temperature data were from the observations of the automatic weather stations (AWS). We appreciate help from Dr.Dexian FANG to provide the radar observations of Nanchang, Ji’an,and Shangrao to verify the simulations. We are grateful to the High-Performance Computing Center (HPCC) of Nanjing University for doing the numerical calculations in this paper on its blade cluster system. We thank the two anonymous reviewers,whose comments improved the manuscript.

    热99re8久久精品国产| 国产亚洲欧美98| 免费看光身美女| 亚洲精品日韩av片在线观看| 国产毛片a区久久久久| 国产 一区 欧美 日韩| 色播亚洲综合网| 99久久人妻综合| 51国产日韩欧美| av在线播放精品| 亚州av有码| 小蜜桃在线观看免费完整版高清| 一个人看视频在线观看www免费| 美女 人体艺术 gogo| 成人漫画全彩无遮挡| 久久欧美精品欧美久久欧美| 男女视频在线观看网站免费| 国产不卡一卡二| 国产精品一区二区在线观看99 | 丝袜喷水一区| 干丝袜人妻中文字幕| 一本精品99久久精品77| 日本一本二区三区精品| 亚洲精品国产av成人精品| 少妇高潮的动态图| 欧美日本视频| 国产成人a区在线观看| 日本黄色片子视频| 老司机影院成人| 久久久久九九精品影院| 一级毛片我不卡| 岛国在线免费视频观看| 丰满人妻一区二区三区视频av| 能在线免费观看的黄片| 欧美日韩乱码在线| 久久久精品94久久精品| 男女边吃奶边做爰视频| 激情 狠狠 欧美| 一级二级三级毛片免费看| 亚洲乱码一区二区免费版| 中文在线观看免费www的网站| 日韩一区二区三区影片| 精品人妻视频免费看| 久久99精品国语久久久| 亚洲精品日韩av片在线观看| 国产不卡一卡二| 亚洲av男天堂| 亚洲无线在线观看| 激情 狠狠 欧美| 女人十人毛片免费观看3o分钟| av卡一久久| 亚洲最大成人手机在线| 男插女下体视频免费在线播放| 一本久久精品| 亚洲人与动物交配视频| 久久人人爽人人片av| 免费看av在线观看网站| 黄色视频,在线免费观看| 不卡一级毛片| 一边摸一边抽搐一进一小说| 国产黄片视频在线免费观看| 熟女人妻精品中文字幕| 中文在线观看免费www的网站| av在线亚洲专区| 一区二区三区高清视频在线| 少妇丰满av| 丰满乱子伦码专区| 偷拍熟女少妇极品色| 免费搜索国产男女视频| 久久人人爽人人爽人人片va| 啦啦啦韩国在线观看视频| 国产亚洲av片在线观看秒播厂 | 亚洲美女搞黄在线观看| 亚洲av电影不卡..在线观看| 天天躁日日操中文字幕| 小蜜桃在线观看免费完整版高清| 波多野结衣巨乳人妻| 久久精品国产亚洲av涩爱 | a级一级毛片免费在线观看| 成人av在线播放网站| 国产探花极品一区二区| 午夜视频国产福利| 国产探花极品一区二区| 99视频精品全部免费 在线| 99久久久亚洲精品蜜臀av| 精品熟女少妇av免费看| 身体一侧抽搐| 亚洲真实伦在线观看| 久久精品91蜜桃| 亚洲乱码一区二区免费版| 97超视频在线观看视频| 淫秽高清视频在线观看| 国产精品女同一区二区软件| 中文精品一卡2卡3卡4更新| 91狼人影院| 在线天堂最新版资源| 国内揄拍国产精品人妻在线| 校园人妻丝袜中文字幕| 亚洲成人久久爱视频| 亚洲天堂国产精品一区在线| 亚洲,欧美,日韩| 日韩在线高清观看一区二区三区| 亚洲人成网站高清观看| 国产精品蜜桃在线观看 | 狂野欧美白嫩少妇大欣赏| 女人被狂操c到高潮| 国产 一区 欧美 日韩| 久久久久久久久中文| 国产精品人妻久久久久久| 啦啦啦观看免费观看视频高清| av在线老鸭窝| 国产伦一二天堂av在线观看| 国产毛片a区久久久久| 欧美日韩一区二区视频在线观看视频在线 | 日韩av不卡免费在线播放| 亚洲乱码一区二区免费版| 男女做爰动态图高潮gif福利片| 日韩人妻高清精品专区| 日产精品乱码卡一卡2卡三| 九九久久精品国产亚洲av麻豆| 国产一区二区激情短视频| 岛国在线免费视频观看| 成人性生交大片免费视频hd| 又爽又黄无遮挡网站| 亚洲高清免费不卡视频| 日本-黄色视频高清免费观看| 精品一区二区三区人妻视频| 亚洲激情五月婷婷啪啪| 日韩制服骚丝袜av| 97在线视频观看| 国产精品爽爽va在线观看网站| 少妇被粗大猛烈的视频| 亚洲欧美日韩高清在线视频| 久久久欧美国产精品| 国产精品.久久久| 在线观看av片永久免费下载| 九九久久精品国产亚洲av麻豆| 色哟哟·www| 日韩强制内射视频| 久久韩国三级中文字幕| 韩国av在线不卡| 国产成人aa在线观看| av在线蜜桃| 国产伦在线观看视频一区| 亚洲人成网站在线播放欧美日韩| 国产成人精品久久久久久| 国产女主播在线喷水免费视频网站 | 一进一出抽搐gif免费好疼| av免费观看日本| 国产精品福利在线免费观看| 欧美成人精品欧美一级黄| 国产成人a区在线观看| 国产精品.久久久| 国产精品一区二区性色av| 国产精品久久久久久精品电影| 97超视频在线观看视频| 超碰av人人做人人爽久久| 精品午夜福利在线看| 亚洲欧美日韩东京热| 国产老妇女一区| 免费看a级黄色片| 日韩三级伦理在线观看| 少妇裸体淫交视频免费看高清| 好男人在线观看高清免费视频| 校园春色视频在线观看| 国产在线男女| 免费无遮挡裸体视频| 色哟哟哟哟哟哟| 国产真实伦视频高清在线观看| 亚洲国产精品成人久久小说 | 天堂av国产一区二区熟女人妻| 内射极品少妇av片p| 亚洲欧美精品自产自拍| 久久精品国产亚洲av香蕉五月| 久久综合国产亚洲精品| 久久久国产成人免费| 国产午夜福利久久久久久| 国产精品一区二区在线观看99 | 欧美一区二区亚洲| 国产亚洲av嫩草精品影院| 国内久久婷婷六月综合欲色啪| 女同久久另类99精品国产91| 国产黄色小视频在线观看| 精品99又大又爽又粗少妇毛片| 小说图片视频综合网站| 我要搜黄色片| 大又大粗又爽又黄少妇毛片口| 91aial.com中文字幕在线观看| 一级毛片我不卡| 人人妻人人看人人澡| 26uuu在线亚洲综合色| 久久久色成人| 国产黄片美女视频| 综合色av麻豆| 美女高潮的动态| 国产精品人妻久久久久久| 成人毛片60女人毛片免费| 日韩欧美精品v在线| 国产在线精品亚洲第一网站| 国产午夜精品论理片| 亚洲第一区二区三区不卡| 日韩人妻高清精品专区| 亚洲无线在线观看| 最近的中文字幕免费完整| 国内精品久久久久精免费| 一个人看的www免费观看视频| 久久久成人免费电影| 嫩草影院入口| 久久精品久久久久久久性| 中文资源天堂在线| 波多野结衣高清无吗| 舔av片在线| 成人性生交大片免费视频hd| 成人美女网站在线观看视频| 日韩欧美 国产精品| 欧美一区二区国产精品久久精品| 免费在线观看成人毛片| 特级一级黄色大片| 亚洲精品日韩av片在线观看| 欧美zozozo另类| 桃色一区二区三区在线观看| 又黄又爽又刺激的免费视频.| 在线a可以看的网站| 亚洲美女视频黄频| 69av精品久久久久久| 免费看av在线观看网站| 免费看光身美女| 亚洲成人精品中文字幕电影| 在线免费观看的www视频| 日韩一区二区视频免费看| 观看免费一级毛片| 国产成人freesex在线| 变态另类成人亚洲欧美熟女| 91麻豆精品激情在线观看国产| 午夜老司机福利剧场| 色尼玛亚洲综合影院| 网址你懂的国产日韩在线| 日本色播在线视频| 最近中文字幕高清免费大全6| 级片在线观看| 女同久久另类99精品国产91| 菩萨蛮人人尽说江南好唐韦庄 | 一区二区三区高清视频在线| 变态另类丝袜制服| 99热6这里只有精品| 天天一区二区日本电影三级| 最新中文字幕久久久久| 欧美激情在线99| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 中文字幕熟女人妻在线| 秋霞在线观看毛片| 在线播放无遮挡| 日本色播在线视频| 日日摸夜夜添夜夜爱| 亚洲精品456在线播放app| 老熟妇乱子伦视频在线观看| 亚州av有码| 亚洲人与动物交配视频| 国产成人freesex在线| 又粗又爽又猛毛片免费看| 免费av毛片视频| 在线观看66精品国产| 你懂的网址亚洲精品在线观看 | 全区人妻精品视频| 欧美bdsm另类| 精品免费久久久久久久清纯| 最近2019中文字幕mv第一页| 欧美最黄视频在线播放免费| 成年女人看的毛片在线观看| 精品午夜福利在线看| 欧美在线一区亚洲| 免费看日本二区| 日韩,欧美,国产一区二区三区 | 国产精品一区二区在线观看99 | 啦啦啦啦在线视频资源| 国产乱人偷精品视频| 国产精品一区二区在线观看99 | 亚洲av一区综合| 精品国产三级普通话版| 美女内射精品一级片tv| 国产色爽女视频免费观看| 日韩一区二区视频免费看| 欧美成人免费av一区二区三区| 久久亚洲国产成人精品v| 只有这里有精品99| 狂野欧美激情性xxxx在线观看| 久99久视频精品免费| 久久精品夜色国产| 嫩草影院新地址| 熟女人妻精品中文字幕| 天堂√8在线中文| 性插视频无遮挡在线免费观看| 亚洲第一电影网av| 三级男女做爰猛烈吃奶摸视频| 久久精品夜夜夜夜夜久久蜜豆| 嫩草影院精品99| 欧美成人精品欧美一级黄| 国产精品一区二区三区四区免费观看| 99国产极品粉嫩在线观看| 丰满人妻一区二区三区视频av| 国产亚洲精品久久久com| 国产综合懂色| 亚洲人成网站高清观看| 欧美成人精品欧美一级黄| 国产精品1区2区在线观看.| 亚洲综合色惰| 最近2019中文字幕mv第一页| 亚洲国产日韩欧美精品在线观看| 国产精品无大码| 午夜亚洲福利在线播放| 成人鲁丝片一二三区免费| 少妇高潮的动态图| 久久这里只有精品中国| 国产亚洲精品久久久com| 天堂网av新在线| 日本-黄色视频高清免费观看| 日本成人三级电影网站| 欧美精品国产亚洲| 国产成人福利小说| 六月丁香七月| 啦啦啦观看免费观看视频高清| 亚洲一区二区三区色噜噜| 色综合站精品国产| 国产亚洲欧美98| 少妇高潮的动态图| 久久精品夜夜夜夜夜久久蜜豆| 一边亲一边摸免费视频| 91av网一区二区| 国产精品乱码一区二三区的特点| 久久久欧美国产精品| 久久久精品大字幕| 国产女主播在线喷水免费视频网站 | 九色成人免费人妻av| 成人性生交大片免费视频hd| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 亚洲精品久久国产高清桃花| 99精品在免费线老司机午夜| 又粗又爽又猛毛片免费看| 久久综合国产亚洲精品| 久99久视频精品免费| 女人十人毛片免费观看3o分钟| 女人被狂操c到高潮| 一级黄色大片毛片| 成年女人永久免费观看视频| 免费看av在线观看网站| 日本成人三级电影网站| 又粗又硬又长又爽又黄的视频 | 亚洲欧美中文字幕日韩二区| 大又大粗又爽又黄少妇毛片口| 日产精品乱码卡一卡2卡三| 日日干狠狠操夜夜爽| 亚洲av不卡在线观看| 白带黄色成豆腐渣| 亚洲av免费在线观看| 一本久久中文字幕| 青春草视频在线免费观看| 国产精品综合久久久久久久免费| 久久久久久国产a免费观看| 欧美性猛交╳xxx乱大交人| 男女边吃奶边做爰视频| 少妇熟女欧美另类| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品国产三级国产av玫瑰| 人妻系列 视频| 国产亚洲精品久久久久久毛片| 中文字幕熟女人妻在线| 免费看日本二区| 欧美区成人在线视频| 精华霜和精华液先用哪个| 亚洲精品久久久久久婷婷小说 | 美女脱内裤让男人舔精品视频 | 日韩成人av中文字幕在线观看| 亚洲婷婷狠狠爱综合网| 麻豆一二三区av精品| 国产不卡一卡二| 亚洲成人久久爱视频| 不卡视频在线观看欧美| 91麻豆精品激情在线观看国产| 中文字幕av成人在线电影| 中国国产av一级| 国产又黄又爽又无遮挡在线| 国产精品福利在线免费观看| 国产成人影院久久av| a级一级毛片免费在线观看| 深夜精品福利| 午夜福利成人在线免费观看| 人人妻人人澡欧美一区二区| av女优亚洲男人天堂| 久久久久免费精品人妻一区二区| 国产精品,欧美在线| 亚洲精品456在线播放app| av在线蜜桃| 又粗又硬又长又爽又黄的视频 | avwww免费| 国产精品久久久久久av不卡| 免费一级毛片在线播放高清视频| 日本色播在线视频| 久99久视频精品免费| 国产精品一区二区三区四区免费观看| 久久中文看片网| 日本撒尿小便嘘嘘汇集6| 你懂的网址亚洲精品在线观看 | 亚洲成人精品中文字幕电影| 欧美又色又爽又黄视频| 一区二区三区四区激情视频 | 日韩一本色道免费dvd| 亚洲精品影视一区二区三区av| 成人欧美大片| 99久久精品热视频| 麻豆成人午夜福利视频| 白带黄色成豆腐渣| 国产精品日韩av在线免费观看| 久久婷婷人人爽人人干人人爱| 韩国av在线不卡| 成人高潮视频无遮挡免费网站| 国产黄a三级三级三级人| 观看免费一级毛片| 精品久久久久久久久亚洲| 日韩亚洲欧美综合| 99久久九九国产精品国产免费| 国产一区二区三区在线臀色熟女| 国产一级毛片在线| 99久久无色码亚洲精品果冻| 国内久久婷婷六月综合欲色啪| 亚洲精品乱码久久久久久按摩| 亚洲va在线va天堂va国产| 91狼人影院| 变态另类丝袜制服| 级片在线观看| 久久精品国产自在天天线| a级毛色黄片| 身体一侧抽搐| 一夜夜www| 国产av麻豆久久久久久久| 男女下面进入的视频免费午夜| 人妻少妇偷人精品九色| 亚洲性久久影院| 亚洲国产欧洲综合997久久,| 中出人妻视频一区二区| 亚洲av中文av极速乱| 久久久久网色| 久久久久久久久久久免费av| 欧美三级亚洲精品| 一进一出抽搐gif免费好疼| 中文字幕av在线有码专区| 国产黄色视频一区二区在线观看 | 少妇猛男粗大的猛烈进出视频 | 国产成年人精品一区二区| 久久精品久久久久久久性| 人妻夜夜爽99麻豆av| 18禁在线无遮挡免费观看视频| 欧美激情国产日韩精品一区| 国产精品日韩av在线免费观看| 乱系列少妇在线播放| 国产精品嫩草影院av在线观看| 尾随美女入室| 只有这里有精品99| 国产精品无大码| 亚洲精品国产av成人精品| 亚洲成av人片在线播放无| 色吧在线观看| 免费黄网站久久成人精品| 自拍偷自拍亚洲精品老妇| 亚洲成人久久爱视频| 色哟哟哟哟哟哟| 欧美成人a在线观看| 免费看光身美女| 欧美日韩乱码在线| 欧美高清成人免费视频www| 日韩三级伦理在线观看| 别揉我奶头 嗯啊视频| 男女做爰动态图高潮gif福利片| 如何舔出高潮| 精品人妻视频免费看| 热99re8久久精品国产| 欧美三级亚洲精品| 亚洲国产高清在线一区二区三| 久久亚洲精品不卡| 亚洲欧美精品综合久久99| 男女啪啪激烈高潮av片| 国产淫片久久久久久久久| 亚洲中文字幕日韩| 欧美xxxx黑人xx丫x性爽| 亚洲经典国产精华液单| 成人午夜高清在线视频| 亚洲欧美成人精品一区二区| 色哟哟·www| 久久久久久久亚洲中文字幕| 亚洲乱码一区二区免费版| 成人毛片a级毛片在线播放| 欧美日韩乱码在线| 97超碰精品成人国产| av又黄又爽大尺度在线免费看 | 观看美女的网站| av专区在线播放| 国产精品国产三级国产av玫瑰| 中文字幕免费在线视频6| av在线蜜桃| 国产一区二区三区av在线 | 成人欧美大片| 亚洲精品成人久久久久久| 久久精品国产亚洲av香蕉五月| 日本撒尿小便嘘嘘汇集6| 男的添女的下面高潮视频| 99在线视频只有这里精品首页| 亚洲第一区二区三区不卡| 精品久久久久久久末码| 99久久九九国产精品国产免费| 床上黄色一级片| 亚洲欧美日韩高清在线视频| 久久精品国产亚洲av香蕉五月| 亚洲美女搞黄在线观看| 婷婷六月久久综合丁香| 91精品国产九色| 久久九九热精品免费| 色5月婷婷丁香| 久久综合国产亚洲精品| 国产精品人妻久久久久久| 亚洲无线在线观看| 欧美变态另类bdsm刘玥| 黑人高潮一二区| 夜夜夜夜夜久久久久| 麻豆成人午夜福利视频| 一级黄色大片毛片| 中国美女看黄片| 插逼视频在线观看| 春色校园在线视频观看| 夜夜看夜夜爽夜夜摸| 日韩高清综合在线| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 岛国毛片在线播放| 夜夜爽天天搞| 伦理电影大哥的女人| 毛片一级片免费看久久久久| 国语自产精品视频在线第100页| 男女视频在线观看网站免费| 婷婷精品国产亚洲av| 一边摸一边抽搐一进一小说| 成年女人看的毛片在线观看| 国内精品宾馆在线| 日韩欧美在线乱码| 伦精品一区二区三区| 亚洲精品国产成人久久av| 国产综合懂色| 免费av毛片视频| 中文在线观看免费www的网站| 亚洲自偷自拍三级| 美女被艹到高潮喷水动态| 少妇丰满av| av在线天堂中文字幕| 亚洲第一电影网av| 国产91av在线免费观看| 亚洲欧美精品专区久久| 久久草成人影院| 午夜a级毛片| 国内久久婷婷六月综合欲色啪| 国产蜜桃级精品一区二区三区| 日韩欧美国产在线观看| 丰满人妻一区二区三区视频av| 一本精品99久久精品77| 欧美+日韩+精品| 只有这里有精品99| 久久草成人影院| 看非洲黑人一级黄片| 成人亚洲精品av一区二区| 国产高清激情床上av| 中文在线观看免费www的网站| 久久中文看片网| 精品国内亚洲2022精品成人| 男女边吃奶边做爰视频| 精品久久久久久久人妻蜜臀av| 亚洲最大成人av| 免费看光身美女| 欧美激情在线99| .国产精品久久| 欧美成人一区二区免费高清观看| h日本视频在线播放| 日韩欧美精品免费久久| 国产成人精品久久久久久| 成人性生交大片免费视频hd| 日日啪夜夜撸| 久久精品国产99精品国产亚洲性色| 3wmmmm亚洲av在线观看| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩无卡精品| 亚洲成人精品中文字幕电影| 在线免费观看不下载黄p国产| 免费看a级黄色片| 特大巨黑吊av在线直播| 热99在线观看视频| 国产精品永久免费网站| 久久婷婷人人爽人人干人人爱| 99热这里只有是精品50| 国产一区二区亚洲精品在线观看| 日韩在线高清观看一区二区三区| 亚洲av第一区精品v没综合| 久久久久久久久久久免费av| 毛片女人毛片| 亚洲精品日韩在线中文字幕 | 色5月婷婷丁香| 我要看日韩黄色一级片| 麻豆成人av视频| 欧美日韩乱码在线| 日韩,欧美,国产一区二区三区 | 亚洲精品影视一区二区三区av| 三级男女做爰猛烈吃奶摸视频| 国产亚洲精品久久久com| 大型黄色视频在线免费观看| 国产亚洲欧美98| 日韩精品有码人妻一区| 亚洲美女搞黄在线观看| 亚洲综合色惰|