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

    On the Mechanism of a Terrain-Influenced Snow Burst Event during Midwinter in Northeast China

    2021-04-20 00:42:10NaLIBaofengJIAOLingkunRANXinyongSHENandYanbinQI
    Advances in Atmospheric Sciences 2021年5期

    Na LI, Baofeng JIAO*,, Lingkun RAN, Xinyong SHEN, and Yanbin QI

    1Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

    2Key Laboratory of Meteorological Disaster, Ministry of Education/Joint International Research Laboratory of Climate and Environment Change/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science and Technology, Nanjing 210044, China

    3Jilin Weather Modification Office, Changchun 130062, China

    ABSTRACT Short-duration snow bursts with heavy snow represent one type of hazardous weather in winter which can be easily missed by the winter weather warnings but often results in great hazards. In this paper, the mechanism for the occurrence of such events was investigated with the aid of a localized terrain-influenced snow burst event in Northeast China. The snow burst was produced by an eastward-moving cold-frontal snowband which encountered the downstream complex terrain of the Changbai Mountains and intensified. To ascertain the role of orography on the snow burst, numerical experiments,together with a parallel sensitivity experiment removing Changbai Mountains, were performed to attempt to distinguish the contributions of cold-frontal system and orographic effects to produce the heavy snow. Diagnosis showed that without the influence of Changbai Mountains, the release of conditional instability (CI) and inertial instability (II) within a weak frontogenetical environment was responsible for the snowband maintenance. Orographic effects played important roles in enhancing the snowband and increasing the snowfall intensities. The enhancement mechanism was related to the interactions of the cold-frontal snowband and the topography. On the one hand, orographic frontogenesis and persistent ascent, created by orographic gravity waves over the terrain, greatly enhanced the orographic lifting. The intensification of the lifting promoted the release of CI and thus enhanced the snowfall. On the other hand, pre-existing orographic instabilities were released due to the passing of the cold-frontal snowband, which could also serve to intensify the snowband over terrain and thus increase the snowfall.

    Key words: terrain-influenced snow burst, snowband, instability, lifting

    1. Introduction

    Heavy snow represents one type of major catastrophic weather event that occurs at middle and high latitudes during the cold season. Heavy snow events and their accompanying high winds, persistent low temperatures, and low visibilities typically have severe impacts on ecosystems, socioeconomic factors, human lives, and property (Changnon and Changnon, 2006; Wang et al., 2010). Consequently, extensive research has been conducted on heavy snow during the past century (e.g., Bennetts and Hoskins, 1979; Sanders and Bosart, 1985; Ulbrich, 2001).

    Heavy snow is frequent in regions of China during wintertime, such as the Tibetan Plateau, Xinjiang, Northeast China, Shandong Peninsula, and even low-latitude South China. Research on heavy snow events in China began in the 1970s, and early work mainly involved synoptic-scale studies. For example, Zhu and Shou (1994) studied the formation of a secondary circulation along the fronts of winter heavy snow events. Tao and Wei (2008) researched the continuous occurrence of freezing rain and snow in southern China from mid-January to early February 2008 and discovered that this persistent precipitation resulted from the activity of a cut-off low from Europe and the transport of abundant warm and moist air from South China Sea, forming a melting layer in Southern China. Since the late 20th century, increasing attention has been paid to the mesoscale features of snowfall with the development of numerical models and observation techniques. Asuma et al. (1997) focused on the properties of winter monsoon clouds using dual-polarization Doppler radar data, and Wang and Cheng (2000,2001) indicated that the snowstorm of January 1996 developed under symmetric instability (SI) conditions based on mesoscale model simulations.

    Nevertheless, previous studies focused mainly on snowfall events with long durations and significant accumulation. In recent years, snow squalls or snow bursts that are characterized by short durations but heavy snow received more attention. The concern of such events stems from the fact that they usually do not meet the warning criteria but frequently create hazardous weather. As one type of a linear convection event, snow squalls share some common characteristics with warm season squall lines and may lead to extremely low visibility as well as hazardous driving conditions (Pettegrew et al., 2009). Market et al. (2002) presented a climatological analysis of thundersnow and showed that most thundersnow events in U.S. are associated with a transient mid-latitude cyclone, whereas the remaining events correspond to orographic forcing, lake effect snow,or, on rare occasions, fronts. Liu and Moore (2004) studied the climatic characteristics of lake effect snowstorms in southern Ontario, Canada, during 1992?99, and noted that lake effect snow occurred in the presence of steep lapse rates and subsequent CI. Pettegrew et al. (2009) researched a winter convection case in the Midwestern United States that produced 4 cm of snow in a short time and created nearwhiteout conditions; the shallow forced convection was related to frontogenesis and convective instability in a sheared environment. Milrad et al. (2014) compared two snow bursts that were not associated with lake effect or mid-latitude cyclones, thus providing a good comparison of different types of instabilities which exist in unsaturated and saturated environments.

    The above studies greatly enhance our understanding of both the short duration and the persistent heavy snow events, which are beneficial for improving predictive skill.Nevertheless, to the best of our knowledge, within the literature published in China, research involving frontal snowsqualls or snow-burst events in the cold season is still limited despite the fact that hazards induced by them are comparably great, such as the case we will study in this paper. A snow burst event occurred in Northeast China on 26–27 January 2017. During this event, more than 12-mm snow fell in a very short period of time without an adequate forecast or advance warning. As a result, 12 multicar collisions were reported on several highways around Harbin (news from https://www.sohu.com/), claiming the lives of seven people and resulting in significant economic losses. This strongly justifies the great importance of explicitly understanding the dynamic mechanism responsible for the occurrence of snow bursts, which we believe will help to strengthen the earlywarning ability of such short-duration snow events in China.

    Observations shown in Fig. 1 show a snow burst that was associated with a cold-frontal snowband, which formed along the foothills of the Dahingganling Mountains and then moved eastward, before encountering the high and complex terrain of the Changbai Mountains. Heavy snow was produced on the windward side and over the Changbai Mountains, shown by the distributions of the three snowfall centers during the event in Fig. 1: Jilin center in Jilin Province(10 mm, red dot, 43.8°N, 126.5°E); Benxi center in Liaoning Province (12.3 mm, purple dot, 41.3°N, 123.7°E; and centered on the border of Liaoning and Jilin Province(denoted as “BLJ”, 6.2 mm, yellow dot). The question then arises: How was the snow burst produced in those snowfall centers under the influence of complex terrain?

    The importance of orography upon the distribution of snowfall has long been recognized and explored. Adhikari and Liu (2019) analyzed the global geographical distribution of thundersnow events using four years of observations from the Global Precipitation Measuring Mission Ku-band radar. It was found that mountainous regions such as the Himalayas, Tibet, the Andes, and the Zagros have highestfrequency of thundersnow events, but they also occur in continental and coastal regions.

    Fig. 1. Observed cumulative liquid-equivalent precipitation(color dots, units: mm) over 24 hours ending at 0000 UTC on 27 January and the terrain height (shaded, units: m) in Northeast China. The numbers “1”, “2”, “3” and “4” respectively indicate Heilongjiang Province, Jilin Province,Liaoning Province and Inner Mongolia. Dotted arrow indicates the direction of the band motion. Red dot in Jilin province,purple dot in Liaoning province and yellow dot on the border of Jilin and Liaoning province (BLJ) respectively indicate the three snowfall centers of the snow burst.

    Extensive studies have been conducted concerning the orographic enhancement of summer rainfall, however investigations regarding the orographic enhancement of snowfall during cold seasons are much less, especially ones that focus upon the explicit mechanism by which–the orography increases the snowfall. However, the orographic enhancement of both summer rainfall and cold season snowfall share many similar physical processes such as the orographic lifting, orographic frontogenesis and convergence,and the “seeder-feeder” mechanism. Schumacher et al.(2010) analyzed two different kinds of snowbands downstream of elevated terrain initiated by both surface and elevated frontogenesis; notably, the major bands were maintained by moist conditional symmetric instability (CSI) or dry SI and IIs on the right exit region of a jet. Schumacher et al. (2015), a follow up to Schumacher et al. (2010), studied the effect of upstream terrain as it relates to the occurrence and intensity of the bands. Nevertheless, most of the studies have emphasized the passive effects of orography to weather systems while the possible interactive influence of passing weather systems upon orographic processes, which may also play a role to enhance the precipitation, is mostly ignored.

    In this paper, the mechanism that induces the snow burst on 26?27 January 2017 in Northeast China is investigated using the output of a high-resolution numerical experiment and a parallel sensitivity experiment that removes the downstream Changbai Mountains. The purpose of the sensitivity experiment is to separate the processes associated with the cold front and those associated with the orography in order to ascertain: (1) to what extent that the orography enhances the snowfall to produce the snow burst; (2) the process that supports the snowband development without terrain; and (3) the net orographic effects and their role on the enhancement of snowfall over complex terrain. In particular, the interactions between the cold-frontal snowband and orography is emphasized.

    In this paper, section 2 describes the large-scale environmental settings of the snow burst which occurred on 26?27 January 2017 in Northeast China that is central to this study.Section 3 introduces the details and configurations of the high-resolution numerical simulation and the nature of the parallel sensitivity experiment of the snow burst, and, furthermore validates the control experiment and analyzes the extent to which that the orography enhances the snowfall. Section 4 analyzes the output data which sheds light upon the physical mechanisms responsible for the snow burst associated with the interaction of the cold-frontal snowband and orography. Finally, section 5 provides a summary of the work.

    2. Large-scale configurations of the 26?27 January 2017 snow burst event

    The distribution of the snowfall amount during the snow burst on 26?27 January 2017 (UTC) in Northeast China is presented in Fig. 1. The snowfall amount was extracted from conventional surface observations in the form of liquid-equivalent precipitation; the observations were obtained from the National Meteorological Center (NMC)of the Chinese Meteorological Administration (CMA)(Zhou et al., 2018). From Fig.1, one can see a banded structure to the snowfall, during this event, across Jilin and Liaoning Provinces; although the width of the band as seen from the 24-h cumulative liquid-equivalent precipitation is larger than its real width due to the eastward movement of the band (dotted arrow). In addition, it is noted that the banded structure of snowfall in Heilongjiang Province around the Xiaohingganling Mountains is not as evident as that observed in Jilin and Liaoning Provinces because snowfall in the former area was mostly produced by a surface cyclone which can be seen in Fig. 2.

    Figure 2 presents the large-scale environment of the snowband, where the gridded surface and mandatory level analyses were taken from the Global Forecasting System of the National Centers for Environmental Prediction(NCEP/GFS), which has a horizontal grid spacing of 0.5×0.5. Here, the conditions observed at 0600 UTC, 26 January, 2017 were applied in order to analyze the synoptic configuration. This time is coincident with the period that the snowband experienced evident development over Changbai Mountains and produced heavy snow. As shown in Fig. 2, the mid-troposphere at 500 hPa in Northeast China was controlled by an intensifying, eastward-moving trough.Located over eastern Inner Mongolia (see number “4” in Fig. 1), the trough continuously transported cold air (the cold advection area enclosed by the thick red solid line in Fig. 2a) into Northeast China, forming a northeast-southwest-oriented cold front (dotted blue line). Above the frontal boundary, in the upper troposphere at 200 hPa, a remarkable upper-level jet streak (ULJ) was detected with a maximum wind speed exceeding 60 m s. This front and ULJ created a favorable large-scale ascent mechanism, which was the basic synoptic forcing for the precipitation. Consequently, snow occurred in the warm-sector of the cold front, which aligned with the exit region of the ULJ.

    At the surface, Northeast China was under the influence of an extratropical cyclone, whose center was located over Heilongjiang Province (Fig. 2b, see number “1” in Fig. 1), and the cold front extended equatorward from the low-pressure center. A comparison between Fig. 1 and Fig. 2 indicates that the snowfall over the Xiaohingganling Mountains in northeastern Heilongjiang Province was caused by the above-mentioned cyclone, whereas the snowfall in Jilin and Liaoning Provinces (respectively number“2” and “3” in Fig. 1) was due to the cold front. In this paper, we focus predominantly on the snowband along the cold front.

    Fig. 2. (a) Geopotential height (solid blue lines, units: 10 gpm) and temperature (dashed red lines, units: °C) at 500 hPa and wind speeds greater than 40 m s–1 at 200 hPa; (b) sea level pressure (solid lines, units: hPa, interval: 2.5 hPa) and wind vectors with wind speeds (colored shading) greater than 12 m s–1 at 850 hPa; (c) equivalent potential temperature (solid lines, units: K, interval: 2 K) and water vapor flux divergence [shaded areas, units: 10–5 g (hPa m2 s) –1] at 850 hPa analyzed from GFS at 0600 UTC on 26 January 2017; (d) TBB observed by Himawari satellite at 0600 UTC 27 January 2017. The thick black line in (a) is the 500-hPa trough; the dotted blue lines in (a)–(d) are the cold front; the area enclosed by the red solid line in (a) is the cold advection area. The red triangles indicate the ten stations with the largest snowfall totals from 0600 to 1200 UTC on 26 January 2017 to denote the snowband regions.

    The low level, 850 hPa, analysis reveals that the cold front exhibited strong cyclonic wind shear over Northeast China (Fig. 2b). The low-level jet (LLJ) formed by the southwesterly wind in the warm sector of the cold front stretched over the eastern part of Northeast China with a maximum wind speed reaching 27 m s. This LLJ transported adequate heat and water vapor to Northeast China and provided for strong convergence upon collision with the cold, dry northwesterly flow from high latitudes behind the front. As shown in Fig. 2c, water vapor strongly converged along the cold front in the warm sector, with the largest water vapor flux convergence associated with the largest amount of snowfall (indicated by red triangles which denote stations with the largest snowfall totals observed during 0600?1200 UTC on 26 January 2017).

    The satellite imagery (Fig. 2d) at 0600 UTC 27 January 2017 also indicated a convective banded structure. One can find two relatively isolated convective bands, which are moving toward the Northeast China and corresponded with the dual snow-bands mentioned next in Fig. 3a.

    From the above synoptic scale analysis, the snow burst in Northeast China on 26 January 2017 (UTC) was mainly caused by an eastward-moving cold front. The along-front snowband occurred primarily in the warm sector of the cold front, where there was a strong LLJ and considerable moisture convergence. The snowfall centers (red, purple and yellow dots in Fig. 1, see introduction) tended to occur on the windward side of the Changbai Mountains. To explore the details of the combined effects of the front and complex terrain to produce the snow burst and thus reveal the occurrence mechanisms of the snow burst, high-resolution numerical experiments (including control and sensitivity experiments) of this case are conducted in the next section.

    3. Numerical experiments

    3.1. Model configurations

    Fig. 3. Cumulative liquid-equivalent precipitation (units: mm) from observations (first line), the CTL experiment(second line) and TRNP experiment (third line): (a, d, g) 0000?1200 UTC, (b, e, h) 0000?0600 UTC, and (c, f, i)0600?1200 UTC on 26 January 2017. The gray shaded areas in the simulation are terrain heights. The observations were obtained by interpolating the station data to the set of 0.1°×0.1° grids.

    Version 3.8.1 of the Advanced Research Weather Research and Forecasting Model (WRF-ARW; Skamarock et al., 2008) was used to simulate the 26?27 January 2017 snowburst event at convection-appropriate resolutions. The model contained only one domain, which included all of Northeast China and the surrounding regions (34°?54°N,110°?138°E). The model center was located at 45°N,124°E, and the horizontal grid spacing was 4 km. The number of horizontal grid points was 641×561, and the terrain-following vertical coordinate, sigma, had 51 vertical levels from the surface to 50 hPa, with more vertical levels near the surface and fewer levels aloft (Laprise, 1992). The Community Atmospheric Model (CAM) scheme (Collins et al.,2004, NCAR Tech Note) was employed for longwave and shortwave radiation; in addition, the Pleim-Xiu Land Surface Model (PX LSM) (Pleim and Xiu, 1995; Xiu and Pleim, 2001) and the Asymmetrical Convective Model, version 2 (ACM2) planetary boundary layer (PBL) (Pleim,2007) were used. The precipitation process in the model was represented by the Morrison double-moment scheme(Morrison et al., 2009). The hourly fifth global climate reanalysis produced by European Centre for Medium-Range Weather Forecasts (ERA5/ECMWF) with a horizontal grid spacing of 0.25°×0.25° and a vertical resolution of 37 levels was used to initialize the model and produce lateral boundary conditions at 3-h intervals. The experiment was initialized at 1200 UTC on 25 January 2017 and integrated for 24 h. History files were generated once an hour.

    On the basis of the above control simulation with full terrain (denoted as "CTL"), one idealized experiment is further conducted (denoted as "TRNP"). All configurations in the TRNP experiment are the same as the CTL experiment except that all the terrain height greater than 200 m in the Changbai Mountains are set to be 200 m to remove the effects of Changbai Mountains in TRNP (200 m is the mean elevation of the plains in Northeast China). Through this simple sensitivity experiment, one can basically separate the component of snowfall that was attributed solely to the cold-front and thus ascertain to what extent the orography enhanced the snowfall. Additionally, the physical processes associated with the cold front and the orography can be distinguished and roughly separated.

    3.2. Validation of the control simulation

    Figure 3 illustrates the simulated output of cumulative liquid-equivalent precipitation for both the CTL and TRNP experiments for comparative purposes. Since little precipitation fell after 1200 UTC on 26 January within Northeast China, only the precipitation from 0000 to 1200 UTC is shown. Figures 3a?c were obtained by interpolating the station observations into a set of 0.1°×0.1° grids. Due to the interpolation, the location of BLJ center in Fig. 3a is a little different from that in Fig. 1. Figures 3a?f can be used to validate the performance of the simulation while Figs. 3d?i can be used to analyze the enhancement of the snowfall due to the complex terrain.

    Through a comparison of Figs. 3a and 3d, the CTL simulation reproduced the observed snowfall very well, including the detailed banded structure of the precipitation and the three snowfall centers within the snowband. As shown in Fig. 3d, both the Jilin center and the Benxi center were encompassed by narrow small-scale bands (approximately 300 km long and 30 km wide) within the overall snowband; the Jilin center even exhibited a double-banded structure. Notably,the BLJ center was quite evident in the simulation results,as it exhibited a cellular precipitation structure that extended beyond China. The much larger extent of BLJ in the simulation compared to that in the observation was likely due to the sparseness of observation stations in the mountains and the lack of observations outside of China (Fig. 1).

    In addition to the spatial distribution of the snowband,its temporal evolution was also effectively simulated, as shown in Figs. 3b?c and Figs. 3e?f. During the observational period, snowfall was most abundant during 0600?1200 UTC; the circumstances were the same in the simulation, although the simulated precipitation in the Benxi center appeared earlier than the observed precipitation. Nevertheless, Fig. 3 still validates the accuracy of the model in simulating the snowband. The solid performance of the model implies that the physical processes associated with the snow band were simulated reliably, especially regarding both the mesoscale and small-scale processes that cannot be resolved by reanalysis data.

    3.3. Enhancement of snowfall by orography

    To evaluate the role of orography in producing the snow burst, snowfall amounts and distributions from the CTL and TRNP experiments are compared. From Figs. 3d?f and 3g?i, the snowfall area produced by the cold-frontal snowband does not change much after removing the downstream terrain of the Changbai Mountains. However, the snowfall intensity is significantly influenced. In Fig. 3g, all three snowfall centers in Jilin and Liaoning provinces are greatly reduced, especially the Benxi Center in Liaoning Province. Comparison between Figs. 3e?f and 3h?i indicates that the 0600?1200 UTC snowfall contributes to the difference between Figs. 3d and 3g, a time when the snowband was just becoming active over Changbai Mountains.This indicates that the orography of the Changbai Mountains plays an important role in inducing the snow burst or creating heavy snowfall centers.

    To quantitatively evaluate to what extent orography effects enhance the snowfall, Table 1 calculates the numbers of grids within the snowband region for different snowfall intensities. The results of Table 1 are consistent with Fig. 3d?i in that the snowfall area greater than 3 mm (12 h)does not appreciably change in the CTL and TRNP experiments. However, when comes to the > 6 mm and > 9 mm snowfall areas, the orographic effects are quite evident.They contribute 49.7% to the 12-h snowfall greater than 6 mm compared to 26.8% of the cold front, and contributes 66.1% to the 12-h snowfall greater than 9 mm compared to 7.4% of the cold front.

    In the following section, we will use the TRNP experiment to analyze the mechanism that supports the snowband development associated with the cold front. Then, by comparing CTL and TRNP experiments, the orographic effects or the interaction between the cold-frontal snowband and orography that intensifies the snowband can be more clearly observed.

    4. Development mechanism of snowband over the Changbai Mountains

    4.1. Overview of the snowband

    As radar observations were not available to view the snowband development for this case, the simulated composite radar reflectivities are analyzed instead to observe the development process of the snowband. Figure 4 presents the simulated composite radar reflectivity from both CTL and TRNP experiments during the entire snowfall process superimposed upon the terrain around Northeast China. Two types of snowfall structures are found in both experiments:an arc-shaped cellular snow region along the northeastern border of Heilongjiang Province related to the surface cyclone there (Fig. 2b), and a linear, northeast-southwest-oriented snowband that formed along the front that was later influenced by the downstream complex terrain of the Changbai Mountains. The linear snowband, which resulted in considerable snowfall in Jilin and Liaoning Provinces (as shown in Figs. 3d and 3f), is our focus.

    Figures 4a?h clearly display how the snowband was initialized and organized as well as how it developed while moving eastward and climbing over the mountains. The snowband was initiated from a narrow convection line along the foothills of the upstream terrain of the Dahingganling Mountains at approximately 2100 UTC on 25 January. In the first few hours (Figs. 4a?c), these convective lines moved southeastward slowly and produced weak snowfall. This process was observed to be the initial stage of snowband develop-ment, when the upstream terrain influenced the band evolution. At 0400 UTC (Fig. 4d), the snowband experienced major development as it approached the Changbai Mountains. Both spatial dimensions, along and normal to the band, grew perceptibly. Since 0400 UTC (Figs. 4d?h), the snowband evidently became organized along the cold front and gradually passed the Changbai Mountains. As shown in Fig. 4f, at 0600 UTC, two double-banded structures existed within the snowband and were distributed in Jilin and Liaoning Provinces respectively; these smaller-scale lines contributed to the large total amount of snowfall during this snow burst.

    Table 1. Number of grids within the snowband region (36°?46°N, 120°?130°E) for different snowfall intensities.

    Fig. 4. Simulated composite radar reflectivity (colored shading, units: dBZ) and terrain height (gray shading, units:m) on 26 January 2017: (a?h) the CTL experiment; (i?p) the TRNP experiment. The red line in Fig. 4f indicates the line for the following vertical cross sections.

    Compared to CTL, the snowband in TRNP presents a similar development process. Figures 4i?k reflect that the downstream Changbai Mountains do not show evident influence on the initial development of the snowband, which is consistent with the snowfall distributions in Figs. 3e and 3h.After 0400 UTC, it is observed that the effects of the modified terrain to the snowband distribution is small (Figs. 4l?p and Figs. 4d?h). This indicated that orographic effects do not significantly influence the development of the snowband;the mechanism of the initial development mainly comes from the cold-frontal system. However, the area with radar reflectivities greater than 20 dBZ within the snowband is largely reduced when the Changbai Mountains are removed,indicating a constructive role of the complex terrain in enhancing the snowfall intensities (Figs. 4e?f and 4m?n).

    4.2. Snowband dynamic structure

    Before addressing the mechanisms that support the development and intensification of the snowband, Fig. 5 first presents the dynamic structures of the band in both the CTL and TRNP experiments at 0600 UTC, which could help us to more clearly understand its development. As in Fig. 4f,the line, which the vertical-sections go through, crosses the Benxi and BLJ snowfall centers.

    In the CTL experiment, as in Figs. 5a?c, an evident transverse vertical circulation associated with the ULJ can be found in the upper levels (shaded arrows). The low-level frontal circulation (thick black curved arrow) was largely influenced by the ULJ transverse circulation. The downward motions (Figs. 5a?b) of the frontal circulation occurred mostly along the saturated equivalent potential temperature surfaces and were confined to the frontal zone(blue dotted lines). The widespread upward motion in the prefrontal areas was combined with the ascending branch of the ULJ (Figs. 5b?c), thus providing a large-scale ascent mechanism for the snowband. In addition to this primary circulation related to the jet-front system, another prevalent flow pattern was present in the form of alternatively distributed rising and sinking motions (Fig. 5c) over the terrain, in the warm sector of the cold front. These motions were generated by the orographic gravity waves when air flows over the mountains. The orographic gravity waves are mixed with the upper-level gravity waves produced by ULJ, forming deep ascending and descending motion over the terrain.

    Fig. 5. Vertical cross-sections of (a, d) the saturated equivalent potential temperature (black solid line, units: K) and simulated radar reflectivity (shaded areas, units: dBZ); (b, e) the total wind speed (shaded areas, units: m s–1); and (c,e) the vertical velocity (shaded areas, units: 10–2 m s–1) along the line shown in Fig. 4f at 0600 UTC on 26 January 2017. The left column comes from CTL experiment, while the right column comes from TRNP experiment. The thick solid lines with 10-dBZ and 20-dBZ radar reflectivities depict the location of the snowband. The symbols “L1”,“L2”, “L3” indicate convective cells. "ULJ" is the upper-level jet. The blue dotted lines indicate the frontal zone. The black curved arrows are the frontal circulation, and the shaded arrows are the ULJ transverse vertical circulation. The red stars indicate locations of ascentassociated with the orographic waves.

    The snowband (CTL experiment), indicated by its high reflectivity values, was located at the leading edge of the front corresponding to a wind shift at low levels. At 0600 UTC, three distinct convective cells can be seen. Each cell,according to its location relative to the front and the topography, had a different dynamic structure. The cell (denoted as “L1”) at the back of the snowband was above the lowlevel front and featured a slanted ascent embedded within the frontal circulation. Figure 4f shows stratiform precipitation at the back of the snowband corresponding to L1. The cell in the area where the surface front is confronted with topography (denoted as “L2”) includes two different types of ascent: the western part of the cell near the frontier of the cold front displayed a slanted ascent, while the eastern part over the terrain exhibited a nearly-upright ascent (note that the term "upright" here is a relative concept since the vertical velocity increased by a factor of 10). The cell (denoted as “L3”) in front of the snowband was entirely over the terrain with alternatively distributed upward and downward motions.

    It is noted that the cells L2 and L3 in Figs. 5a, c, e actually correspond to the smaller-scale convective lines within the snowband in Fig. 4f and they respectively induce the large snowfall of the Benxi and BLJ centers. These cells also appear in the TRNP experiment but with weaker intensities. Upon comparing Figs. 5a, c, e and 5b, d, f the difference of the snowband dynamic structure between CTL and TRNP is mainly attributed to differential vertical velocities.Without the influence of Changbai Mountains, only the frontal circulation is present around the snowband, characterized by slantwise ascent. The alternatively distributed upward and downward motions, due to the orographic gravity waves, are not present in the TRNP snowband.

    4.3. Snowband maintenance and intensification mechanism over terrain

    For banded precipitation related to frontal zones, the CSI associated with frontogenesis-based lifting may be the predominant mechanism for explaining the initiation, organization, and development of the band (Bennetts and Hoskins,1979). The conditions for CSI require the moist potential vorticity (MPV) to be negative when the atmosphere is inertially and conditionally stable (reviewed by Schultz and Schumacher, 1999). In some cases, when the saturation condition is not satisfied, dry SI and inertial instability (II) can also appear, thereby generating a sloped ascent and initiating a banding structure, as in the cases discussed by Schultz and Knox (2007) and Schumacher et al. (2010). The conditions for II require the absolute vorticity to be negative,while the conditions for dry SI require the atmosphere to be both inertially and gravitationally stable and for the potential vorticity (PV) to be negative.

    In addition, Schultz and Knox (2007) revealed the simultaneous presence of CI, II and SI in real cases, a phenomenon also discussed by Schumacher et al. (2010). Following Emanuel (1980), Jascourt et al. (1988) defined the coexistence of CI and CSI as convective-symmetric instability,although the mesoscale circulation related to this situation remains unclear. Xu (1986) discussed the interaction between CI and SI in a prefrontal environment and proposed two mechanisms for the development of rainbands based on the sequence of CI and SI: "upscaling" and "downscaling". More recently, Morcrette and Browning (2006) proposed a release of SI caused by ΔM-adjustment, which could also lead to slanted circulation and precipitation bands near fronts.

    In this subsection, these mechanisms associated with instabilities and corresponding liftings for banded convection are examined for the current snowband. The criteria used to evaluate the instabilities are as follows:

    where (θ)=dθ/d

    z

    , (θ)=dθ/d

    z

    , θis the equivalent potential temperature, θ is the potential temperature, ζis the absolute vertical vorticity, R H denotes the relative humidity, PV is potential vorticity defined by PV=(ω·?θ)/ρ,MPVis the moist potential vorticity defined by MPV=(ω·?θ)/ρ , ρ is the air density, ωis the threedimensional vorticity vector. In addition, the total wind was used to calculate ζ, PV and MPV (Novak et al., 2004,2006; Schultz and Knox, 2007).4.3.1.

    Maintenance mechanism of the snowband

    Figure 6 presents horizontal maps of the different types of instabilities and frontogenesis at 0600 UTC in both the CTL and TRNP experiments. The right panel of Fig. 6 comes from TRNP experiment for the purpose of examining the mechanism responsible for the maintenance of the snowband. The TRNP experiment is analyzed first, because the main process in maintaining snowband development can be deduced in the absence of orographic effects. Conversely, the orographic effects that enhance the snowfall can also be ascertained by extracting the maintenance mechanism from the CTL experiment.

    Fig. 6. (a, e) CI (green shading) and CSI (yellow shading) areas at a model level of η=0.97; (b, f) same as (a, e) but at a model level of η=0.88; (c, g) II (contours) and negative PV (gray shading) areas at a model level of η =0.97; (d, h) same as (c, g) but at a model level of η=0.79 at 0600 UTC on 26 January 2017. The left column comes from the CTL experiment, while the right column comes from the TRNP experiment.The black solid lines are 10-dBZ and 20-dBZ radar reflectivity contours denoting the snowband region. The blue box encloses the II area within the band.

    In Figs. 6a?d, all types of instabilities listed in Eq. (1)appeared within and around the TRNP snowband region(denoted by the 10-dBZ radar reflectivity contour). As shown in Fig. 6a, CI (green shaded areas) was the most widely distributed but only in the near-surface layer, with a small corresponding convective available potential energy(CAPE; not shown). CSI (yellow shaded areas), which is usually considered to be the mechanism responsible for banded convection, also appear, are distributed along the snowband.Dry SI with negative PV and positive absolute vorticity(gray shaded areas) was only present at the near-surface level within and in front of the band (Fig. 6c). Compared to the SI, as shown in Figs. 6c?d, II was quite evident in this case; in the near-surface layer, II was distributed predominantly at the southeast part and in front of the convective region (Fig. 6c), while in the low troposphere, one of the II belts (blue box in Fig. 6d) appeared within the band. All these instabilities revealed by Fig. 6 are possible candidates for the development of the snowband and they are embedded in a weak frontogenetical environment (Fig. 7a) which provides a basic lifting mechanism.

    Fig. 7. Frontogenesis in the near-surface level at η=0.97 of (a)the CTL experiment and (b) TRNP experiment at 0600 UTC on 26 January 2017. The black solid lines are 10-dBZ and 20-dBZ radar reflectivity contours denoting the snowband region.

    The presence of all the instabilities does not mean that they all play significant roles to support the development of the snowband, although they are all embedded in a frontogenetical environment. Figure 8 and Figs. 9a?b illustrate the vertical structure of the instabilities and associated frontogenesis around the TRNP snowband corresponding to the dynamic structures of the snowband in Figs. 5d?f. In Fig. 8,corresponding to the convective cells L2 and L3, CI, CSI,and dry SI all appear under the 2-km height, with CI (blue box) the most evident. In Fig. 9c, L2 has a frontogenesis region in the near-surface layer to release CI while L3 has an elevated frontogenesis region at about 1?2 km above ground level. The two frontogenesis regions contribute to release CI to support the development of L2 and L3 and thus support the development of the snowband. In the higher levels, the circumstance is simpler because only II(orange box) appear in those levels. Distributed at the base of L1, the elevated II and corresponding frontogenesis (Fig.9a) contribute to the slanted circulations in L1 which aligns along the cold frontal boundary (Figs. 5d?f).

    The above analysis indicates that the main mechanism that supports the development of the snowband without the influence of terrain, is the release of CI and II in a weak frontogenetical environment. Next, we analyze the processes that enhance the snowband or snowfall in the presence of complex terrain.

    Fig. 8. Vertical cross-sections of (a) the negative absolute vorticity (thin solid lines), negative PV (blue shaded areas) and(b) the CI areas with (θse)z <0 (red solid lines) and CSI areas with MPV <0, (θse)z >0 and RH >90% (blue shaded areas) at 0600 UTC on 26 January 2017. The thick solid lines are the 10-dBZ and 20-dBZ simulated radar reflectivity contours; the gray lines are potential temperatures in (a), and θse in (b). The orange box indicates an elevated II area, and the blue box indicates the near-surface CI area. The cross-sections are along the same lines as those in Fig. 5.

    Fig. 9. Vertical cross-sections of frontogenesis at (a, c) 0400 UTC and (b, d) 0600 UTC on 26 January 2017.The left column comes from TRNP experiment while the right column comes from CTL experiment. The blue boxes indicate the orographic frontogenesis area before the snowband encounters the downstream terrain. The cross-sections are along the same lines as those in Fig. 8.

    4.3.2.

    Mechanism of snowband intensification by terrain

    By comparing the dynamic structures (Figs. 5c and 5f)and distributions of instabilities (Figs. 8a?b and Figs. 10c?d) around the snowband, three aspects become evident with regard to how orographic effects can intensify the snowfall.The first comes from the instabilities as seen in Figs. 6a?d and 6f?i in which regions of CSI, II and even dry SI are all enlarged especially over Changbai Mountains. The second comes from the frontogenesis as in Figs. 7 and 9, while the third comes from the intensification of lifting due to orographic upslope flow and also orographic gravity waves as indicated in section 4.2.

    Figure 10 shows the vertical distributions of the instabilities from 0500 UTC to 0900 UTC along the same line as that in Fig. 5. In this process, the dissipating L1 disappeared at 0800 UTC, accompanied by the release of an II region at a height of approximately 2 km. This process is similar to the TRNP snowband in Fig. 8.

    L2, which produces the large snowfall in Benxi center(Figs. 3d and 4f), first experienced intensification between 0500 UTC and 0600 UTC and then gradually weakened on the windward side of the mountains. Similar to L2 in TRNP experiment, the intensification of L2 was induced mainly by the CI at the low levels, as there was no evidence of either II or CSI appearing within the L2 region (Figs. 10a?d).However, a comparison of Figs. 10d and 8b or Figs. 6a and 6e fails to show intensification of CI in the presence of terrain. Therefore, it is quite likely that the intensification of a lifting mechanism promotes the release of CI and thus makes L2 intensify, inducing the large snowfall in Benxi Center.

    The increase of frontogenetical lifting can be seen in Figs. 7a?b and Figs. 9a?d. This process is easily understood since the blocking of terrain to the windward flow increases the flow convergence and deformation which are two basic driving conditions forfrontogenesis. It is noted that before the snowband encounters the downstream terrain, relatively strong orographic frontogenesis already existed on its windward side (Figs. 9c). It is when the coldfrontal snowband approaches this region of orographic frontogenesis that it begins to trigger the release of CI around the cold front to support the development of L2 and L3 (Fig.9d). Apart from the frontogenetical lifting, the orographic lifting also plays an important role to release CI.

    As indicated in section 4.2, the atmosphere over the terrain was characterized by alternatively upward and downward motions, revealing four evident localized regions of ascent (indicated by the four red stars in Fig. 5c) distributed in the upslope, downslope and valley areas of the local terrain. These ascending motions did not exhibit a propagating characteristic and remained in the same place even before the snowband encounters the Changbai Mountains,which can be seen in Fig. 11. When the snowband passed through these areas, the lifting was considerably intensified,such as for L2. This gravity-wave lifting is much stronger than the lifting induced by orographic windward slope or the frontogenesis. Moreover, when the snowband passed through the descending area beneath the ascending area, the lifting was greatly inhibited, which is why L2 began to weaken after 0700 UTC, even though it was still on the windward of the mountains. The decrease of L2 over the windward region implies that the gravity-wave lifting was likely the main dynamic mechanism that influenced the intensification of L2, compared to the frontogenesis and orographic physical lifting.

    Fig. 10. Vertical cross-sections (left column) of the negative absolute vorticity (thin solid line), negative PV (blue shaded areas) and (right column) of the CI areas with (θse)z <0(red solid lines) and CSI areas with MPV <0 , (θse)z >0 and RH >90% (blue shaded areas) at (a?b) 0500 UTC, (c?d) 0600 UTC, (e?f)0700 UTC, (g?h) 0800 UTC and (i?j) 0900 UTC on 26 January 2017. The thick solid lines are the 10-dBZ and 20-dBZ simulated radar reflectivity contours; the gray lines in the left column are potential temperatures, which are θse. The blue dotted lines indicate the leading edge of the front. The crosssections are along the same lines as those in Fig. 4f.

    Fig. 11. Vertical cross-sections of the vertical velocity (shaded areas, units: 10?1 m s?1) at (a) 0300 UTC and (b) 0400 UTC on 26 January 2017. The cross-sections are along the same lines as those in Figs. 4c and 4h The arrows are vertical circulation along the section The thick solid lines with 10-dBZ and 20-dBZ radar reflectivities depict the location of the snowband.

    The development process of L3, which produced larger snowfall in BLJ center over the elevated terrain provides evidence for another mechanism by which the terrain influenced the snowband. As shown in Figs. 10a?b, at 0500 UTC, when L3 was triggered before L2, several smallscale IIs and CSIs were distributed in the area. In the horizontal maps, there is also a thin layer of CI that is not reflected in the vertical structure (Fig. 9a). At 0700 UTC, L3 developed into a mature convective cell, and the unstable areas within it were considerably reduced, implying the release of these instabilities. A notable contrast can be seen from a comparison of Figs. 10a?b and Figs. 10i?j, which delineate the instability areas before and after the snowband passed, respectively. At 0900 UTC, considerable II and CSI were reduced compared to the levels at 0500 UTC in the path of the snowband between 124.2°E and126.4°E.

    It is difficult to clarify how these instabilities released because much of the cold frontal circulation, slanted or upright convection of the snowband, and orographic circulation were mixed. However, it can be inferred that the ascending area associated with the cold front and the accompanying frontogenesis played important roles when the band passed over the terrain. Before the snowband and the cold front reached the terrain, orographic frontogenesis (Fig. 9b)and ascent (Fig. 11) had already occurred over the terrain.These processes did not release the instabilities; however,after the snowband and the cold front passed, these instabilities were generally not observed. A possible justification of the above process is from Schumache et al. (2010). They noted that the absolute vorticity and potential vorticity may be altered by the convection which causes a neutralization of negative absolute vorticity and release of inertial instability. This process was quite similar to the scenario described above.

    From the above analysis, we can identify two processes that were important to the intensification of the snowband over terrain: the release of CI due to the orographic frontogenesis and gravity-wave lifting and the release of orographic instabilities (including CI, II, and CSI over the terrain) due to the passing of the cold front and snowband over the terrain. The two processes associated with the interaction of a cold-frontal snowband and complex terrain contributed to the snow burst in Benxi and BLJ regions. Finally, it is noted that the above analyses are basically from the vertical sections along the Benxi and BLJ centers in Liaoning provinces. In fact, another set of cross sections are also analyzed and most of the main processes that are associated with the snowband maintenance and intensification over terrain are similar.

    5. Summary

    This paper studied a terrain-influenced snow burst event moved across Northeast China with a gale and brief snowfall during the daytime on 26 January 2017. The events were characterized by a northeast-southwest-oriented convective band along a cold front with several embedded small-scale convective lines. It represents one typical type of winter snow event with short-duration and heavy snowfall rates in China. These events are frequently ignored by the early warnings of disastrous weather but often induce great hazards due to the very low visibilities and dangerous driving conditions. Therefore, investigations on the occurrence mechanisms of such events are of great importance.

    Initially triggered along the foothills of the high-elevation Dahingganling Mountains, the snowband that produced the snow burst moved eastward with the cold front and encountered another region of high terrain, the Changbai Mountains. Heavy snow was produced within the snowband over the complex terrain, which motivates the investigation into the effects of orography to the snow burst.

    The WRF-ARW model was used to reproduce this event, with simulations largely consistent with the observations. Based on the control simulation, a sensitivity experiment that removed the downstream topography of the Changbai Mountains was conducted to separate the contributions of the cold-frontal system and the topography upon the heavy snowfall. A quantitative analysis suggests that the downstream topography plays a significant role on the local enhancement in the snowfall.

    Simulations with high temporal and spatial resolutions were then applied to reveal maintenance mechanisms of the snowband and also its enhancement under the influence of terrain. The corresponding results are summarized into the conceptual models in Fig. 12.

    As in Fig. 12, the instabilities (including CI, SI, CSI and II), which are considered to be possible mechanisms for banded convection, coexisted around the snowband.Without the influence of downstream terrain, it was discovered that the release of II and CI may have contributed to snowband maintenance. The release of II in a weak frontogenesis environment was responsible for the banded convection over the cold front and caused a slanted ascent. Meanwhile, the release of CI with frontogenesis as the lifting mechanism, triggered and supported the convective line that ran ahead of the cold front.

    Fig. 12. Sketch maps of the mechanisms for the occurrence of the snow burst over the Changbai Mountains: (a) maintenance of snowband without influence of mountains; (b) snowband enhancement associated with the interaction of cold front and topography to produce the heavy snow.

    When the downstream terrain was present, two processes are found to be responsible for the intensification of snowfall. On the one hand, the release of CI associated with the cold front due to mesoscale frontogenesis still occurred.This process was greatly enhanced by orographic frontogenesis and gravity waves over terrain which appeared to intensify the orographic lifting. On the other hand, several instabilities (including CI, II, and CSI) simultaneously coexisted over the terrain. When the cold front and the snowband passed over the terrain, these instabilities were mostly released, which was considered to be another mechanism by which the terrain influenced the development of the snowband.

    As a single-event study, this paper provides a new insight that the terrain influences the precipitating band passing over it and thus produces heavy snowfall. It helps us to understand the occurrence of a snow burst event in Northeast China, to understand the importance of the Changbai Mountains during a snow burst, and to provide scientific evidence and insight with the intent of issuing early weather warnings to avoid the accidents induced by such snow burst events.

    Finally, several problems remain that are worth further investigation: (1) What effect did the upstream terrain have on the snowband development? (2) How did the various instabilities form over the terrain? (3) What are the microphysical processes associated with the interaction of the coldfrontal snowband and topography? These questions will be studied in future papers.

    This study was supported by the Strategic Pilot Science and Technology Special Program of the Chinese Academy of Sciences (Grant No. XDA17010105), the National Key Research and Development Project (Grant No.2018YFC1507104), The Key Scientific and Technology Research and Development Program of Jilin Province (Grant No.20180201035SF), and the National Natural Science Foundation of China (Grant Nos. 41875056, 41575065 and 41790471,42075013), and Chongqing Technology Innovation and Application Development Special Key Project (Grant No. cstc2019jscxtjsbX0007).

    国产在线观看jvid| 高清黄色对白视频在线免费看| 国产精品成人在线| 男女国产视频网站| 狂野欧美激情性bbbbbb| 精品一品国产午夜福利视频| 亚洲av欧美aⅴ国产| 一区在线观看完整版| 国产精品 国内视频| 69av精品久久久久久 | 国产精品av久久久久免费| 欧美成狂野欧美在线观看| 日本猛色少妇xxxxx猛交久久| 欧美乱码精品一区二区三区| 伊人久久大香线蕉亚洲五| 在线精品无人区一区二区三| 香蕉国产在线看| 正在播放国产对白刺激| 国产精品国产三级国产专区5o| av欧美777| 成年美女黄网站色视频大全免费| 免费不卡黄色视频| 无限看片的www在线观看| 老熟女久久久| 亚洲人成77777在线视频| av福利片在线| 国产av国产精品国产| 亚洲精品在线美女| 老汉色av国产亚洲站长工具| 男人操女人黄网站| 国产精品一区二区在线不卡| 精品国产超薄肉色丝袜足j| 国产xxxxx性猛交| 人妻一区二区av| 纵有疾风起免费观看全集完整版| 免费av中文字幕在线| 色老头精品视频在线观看| 欧美av亚洲av综合av国产av| 免费高清在线观看日韩| 十分钟在线观看高清视频www| 在线观看舔阴道视频| 日韩制服丝袜自拍偷拍| 精品国产一区二区三区久久久樱花| 亚洲国产av影院在线观看| 精品国产乱子伦一区二区三区 | 国产欧美亚洲国产| 男人操女人黄网站| 亚洲av片天天在线观看| 国产亚洲av片在线观看秒播厂| 国产一区二区激情短视频 | av在线老鸭窝| 99re6热这里在线精品视频| 十八禁人妻一区二区| 国产片内射在线| 99国产精品一区二区三区| 久久久久国产精品人妻一区二区| 国产精品1区2区在线观看. | 日本a在线网址| 窝窝影院91人妻| 搡老乐熟女国产| 亚洲国产中文字幕在线视频| 黄色片一级片一级黄色片| 又大又爽又粗| 日日爽夜夜爽网站| 精品亚洲乱码少妇综合久久| 男女下面插进去视频免费观看| 日韩中文字幕视频在线看片| 国产精品1区2区在线观看. | 又黄又粗又硬又大视频| 97在线人人人人妻| 精品久久久久久久毛片微露脸 | 国产亚洲一区二区精品| 亚洲av美国av| 人人妻人人澡人人看| 真人做人爱边吃奶动态| 国产在线一区二区三区精| 久久精品人人爽人人爽视色| 精品国产乱码久久久久久小说| 欧美大码av| 亚洲熟女毛片儿| 高清视频免费观看一区二区| 少妇人妻久久综合中文| 精品少妇内射三级| 乱人伦中国视频| 成在线人永久免费视频| 国产精品免费大片| 不卡一级毛片| 黄片小视频在线播放| 91国产中文字幕| 丰满人妻熟妇乱又伦精品不卡| 免费日韩欧美在线观看| 久久久久久人人人人人| 女人被躁到高潮嗷嗷叫费观| 午夜两性在线视频| 欧美亚洲日本最大视频资源| 极品少妇高潮喷水抽搐| 亚洲专区中文字幕在线| 成年动漫av网址| 亚洲成av片中文字幕在线观看| 天堂中文最新版在线下载| 欧美日韩视频精品一区| 午夜久久久在线观看| 国产男女内射视频| 啦啦啦啦在线视频资源| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片| 中文字幕精品免费在线观看视频| 国产精品国产av在线观看| 这个男人来自地球电影免费观看| 亚洲色图综合在线观看| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| 啦啦啦中文免费视频观看日本| 午夜精品国产一区二区电影| 狠狠狠狠99中文字幕| 亚洲国产中文字幕在线视频| 韩国高清视频一区二区三区| 黄色视频,在线免费观看| 欧美黄色淫秽网站| 国产精品久久久久成人av| 国产精品麻豆人妻色哟哟久久| 亚洲久久久国产精品| 久久久久久久久久久久大奶| 亚洲欧洲精品一区二区精品久久久| 国产成人啪精品午夜网站| 黑人欧美特级aaaaaa片| 国产成人a∨麻豆精品| 日本五十路高清| 精品免费久久久久久久清纯 | 国产一区二区在线观看av| 亚洲专区字幕在线| 国产一卡二卡三卡精品| 中文字幕人妻熟女乱码| 国产男人的电影天堂91| 国产高清国产精品国产三级| 成人亚洲精品一区在线观看| 久久性视频一级片| 青春草亚洲视频在线观看| 日韩免费高清中文字幕av| 国产高清视频在线播放一区 | 亚洲欧美成人综合另类久久久| 日本av手机在线免费观看| 99热全是精品| 1024香蕉在线观看| 久久久久久久久久久久大奶| 国产精品免费大片| 老司机亚洲免费影院| 亚洲伊人色综图| 男人添女人高潮全过程视频| 大香蕉久久网| a在线观看视频网站| 女人精品久久久久毛片| 丰满人妻熟妇乱又伦精品不卡| 美国免费a级毛片| 老熟女久久久| 啦啦啦 在线观看视频| 桃红色精品国产亚洲av| 狂野欧美激情性xxxx| 国产1区2区3区精品| 不卡一级毛片| 午夜福利,免费看| 在线观看免费午夜福利视频| 国产精品99久久99久久久不卡| 欧美老熟妇乱子伦牲交| 美女扒开内裤让男人捅视频| 午夜两性在线视频| 国产伦人伦偷精品视频| 国产人伦9x9x在线观看| 老熟妇仑乱视频hdxx| 天天添夜夜摸| 亚洲情色 制服丝袜| 婷婷成人精品国产| 永久免费av网站大全| 久久午夜综合久久蜜桃| 美女视频免费永久观看网站| 男女之事视频高清在线观看| 日日夜夜操网爽| 亚洲国产毛片av蜜桃av| 一区在线观看完整版| 咕卡用的链子| 欧美老熟妇乱子伦牲交| 亚洲第一欧美日韩一区二区三区 | 69精品国产乱码久久久| 成年美女黄网站色视频大全免费| 十八禁网站免费在线| 男女之事视频高清在线观看| 久久天堂一区二区三区四区| 国产一级毛片在线| 亚洲欧美激情在线| 精品福利观看| 亚洲午夜精品一区,二区,三区| 下体分泌物呈黄色| 午夜精品久久久久久毛片777| 久久久欧美国产精品| 欧美97在线视频| 日韩电影二区| 国产成人系列免费观看| 男女午夜视频在线观看| av天堂在线播放| 波多野结衣av一区二区av| 美女午夜性视频免费| 国产亚洲av片在线观看秒播厂| 一区二区日韩欧美中文字幕| 午夜福利乱码中文字幕| 午夜免费观看性视频| 国产在线视频一区二区| 丝袜美腿诱惑在线| 69av精品久久久久久 | 久久国产精品大桥未久av| videosex国产| 国产精品九九99| 国产一区二区激情短视频 | 男女国产视频网站| 黑人欧美特级aaaaaa片| 老司机在亚洲福利影院| 国产精品1区2区在线观看. | 黄网站色视频无遮挡免费观看| 国产日韩欧美视频二区| 在线观看www视频免费| 超碰97精品在线观看| 999久久久国产精品视频| 亚洲精品一区蜜桃| 国产深夜福利视频在线观看| 久久精品国产a三级三级三级| 老熟女久久久| 日韩 欧美 亚洲 中文字幕| 他把我摸到了高潮在线观看 | av视频免费观看在线观看| 黄片小视频在线播放| 精品国产乱码久久久久久男人| 国产免费福利视频在线观看| 亚洲九九香蕉| videosex国产| 亚洲av日韩在线播放| 亚洲精品中文字幕一二三四区 | 交换朋友夫妻互换小说| 麻豆av在线久日| 欧美亚洲 丝袜 人妻 在线| 久久久久久亚洲精品国产蜜桃av| 免费久久久久久久精品成人欧美视频| 制服人妻中文乱码| 黄色毛片三级朝国网站| 在线观看免费高清a一片| 欧美精品啪啪一区二区三区 | 亚洲成人免费av在线播放| 搡老岳熟女国产| kizo精华| 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品久久久久5区| av超薄肉色丝袜交足视频| 两性夫妻黄色片| 亚洲成国产人片在线观看| 成人影院久久| 亚洲情色 制服丝袜| 久久综合国产亚洲精品| 一级毛片电影观看| 日日摸夜夜添夜夜添小说| 美女国产高潮福利片在线看| 午夜免费成人在线视频| 国产日韩欧美在线精品| 久久天堂一区二区三区四区| 一区二区三区精品91| 一进一出抽搐动态| 青春草亚洲视频在线观看| 97精品久久久久久久久久精品| 亚洲精品久久成人aⅴ小说| 黄色视频,在线免费观看| 免费观看a级毛片全部| 成人影院久久| 黄网站色视频无遮挡免费观看| 无遮挡黄片免费观看| 久久久国产精品麻豆| 亚洲中文av在线| 91国产中文字幕| 可以免费在线观看a视频的电影网站| 欧美大码av| 国产亚洲精品一区二区www | 国产三级黄色录像| 交换朋友夫妻互换小说| 亚洲人成电影免费在线| 午夜成年电影在线免费观看| 中文字幕高清在线视频| 手机成人av网站| 久久久久久人人人人人| 日韩电影二区| 色94色欧美一区二区| 桃红色精品国产亚洲av| 十分钟在线观看高清视频www| xxxhd国产人妻xxx| 真人做人爱边吃奶动态| 97人妻天天添夜夜摸| 午夜精品国产一区二区电影| 精品视频人人做人人爽| 俄罗斯特黄特色一大片| 久久九九热精品免费| 婷婷丁香在线五月| 亚洲成人手机| 亚洲中文字幕日韩| 国产精品久久久av美女十八| 在线天堂中文资源库| 美女大奶头黄色视频| 亚洲国产av影院在线观看| 亚洲欧美日韩高清在线视频 | 精品卡一卡二卡四卡免费| 啦啦啦啦在线视频资源| 久久天堂一区二区三区四区| 亚洲欧美一区二区三区黑人| 少妇人妻久久综合中文| 精品国内亚洲2022精品成人 | 亚洲精品国产精品久久久不卡| 18在线观看网站| 亚洲国产精品一区三区| 曰老女人黄片| h视频一区二区三区| av网站免费在线观看视频| 91麻豆av在线| 水蜜桃什么品种好| 国产又爽黄色视频| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费视频网站a站| 亚洲国产精品一区二区三区在线| 三级毛片av免费| 亚洲色图综合在线观看| 在线观看www视频免费| 久久亚洲精品不卡| 国产av精品麻豆| 亚洲国产av新网站| 美女福利国产在线| 天堂俺去俺来也www色官网| 国产精品久久久人人做人人爽| 在线看a的网站| 大片电影免费在线观看免费| 亚洲欧美精品综合一区二区三区| 少妇被粗大的猛进出69影院| av网站免费在线观看视频| 日本猛色少妇xxxxx猛交久久| 王馨瑶露胸无遮挡在线观看| 久久久久久人人人人人| 十八禁网站免费在线| 国产欧美日韩一区二区精品| xxxhd国产人妻xxx| 久久久久久久国产电影| 视频区欧美日本亚洲| 超色免费av| 成人影院久久| 三上悠亚av全集在线观看| 亚洲中文av在线| 精品久久久久久电影网| 9191精品国产免费久久| 狠狠精品人妻久久久久久综合| 亚洲视频免费观看视频| 人成视频在线观看免费观看| 老司机午夜十八禁免费视频| 亚洲精品av麻豆狂野| 国产深夜福利视频在线观看| 一个人免费在线观看的高清视频 | 中文字幕最新亚洲高清| 亚洲五月色婷婷综合| 国产不卡av网站在线观看| 亚洲精品自拍成人| 2018国产大陆天天弄谢| 欧美日韩av久久| 搡老岳熟女国产| av超薄肉色丝袜交足视频| 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 人妻一区二区av| av超薄肉色丝袜交足视频| av网站在线播放免费| 女性被躁到高潮视频| 国产精品久久久av美女十八| 久久精品国产亚洲av高清一级| 国产亚洲欧美精品永久| 老司机影院成人| 午夜福利视频在线观看免费| 久久人妻熟女aⅴ| 亚洲国产av影院在线观看| 精品高清国产在线一区| 国产野战对白在线观看| 免费女性裸体啪啪无遮挡网站| 99热网站在线观看| 精品卡一卡二卡四卡免费| 美女午夜性视频免费| 天天操日日干夜夜撸| 高清av免费在线| 成年人午夜在线观看视频| 日韩电影二区| 久久久久久久久久久久大奶| 丰满迷人的少妇在线观看| 热re99久久精品国产66热6| 免费人妻精品一区二区三区视频| 久久久久精品国产欧美久久久 | 人人妻,人人澡人人爽秒播| 美女脱内裤让男人舔精品视频| 成人三级做爰电影| 久久久久久免费高清国产稀缺| 午夜福利,免费看| 国产国语露脸激情在线看| 欧美乱码精品一区二区三区| 久久国产精品男人的天堂亚洲| 大型av网站在线播放| 欧美另类一区| 妹子高潮喷水视频| 亚洲第一青青草原| 宅男免费午夜| 亚洲精品自拍成人| 精品少妇久久久久久888优播| av线在线观看网站| 国产有黄有色有爽视频| 国产亚洲欧美精品永久| 久久国产精品影院| 亚洲精品美女久久久久99蜜臀| 热99re8久久精品国产| 亚洲欧美一区二区三区久久| 在线观看一区二区三区激情| 母亲3免费完整高清在线观看| 两个人看的免费小视频| 不卡一级毛片| 久久久国产精品麻豆| 美女扒开内裤让男人捅视频| 国产精品一二三区在线看| 91成年电影在线观看| 高清av免费在线| 国产高清视频在线播放一区 | 久久热在线av| 久久精品国产a三级三级三级| 久久99一区二区三区| 人人妻人人添人人爽欧美一区卜| av天堂在线播放| 日韩大码丰满熟妇| 中文字幕高清在线视频| 国产主播在线观看一区二区| 成人黄色视频免费在线看| 一区二区三区乱码不卡18| 蜜桃在线观看..| 极品人妻少妇av视频| 国产精品1区2区在线观看. | 99香蕉大伊视频| 亚洲午夜精品一区,二区,三区| 宅男免费午夜| 亚洲第一欧美日韩一区二区三区 | 亚洲av电影在线观看一区二区三区| 超碰97精品在线观看| 久久人妻福利社区极品人妻图片| 国产男女内射视频| 亚洲专区中文字幕在线| 日本a在线网址| 成人影院久久| 97在线人人人人妻| 久久久久久久久久久久大奶| 久久99热这里只频精品6学生| 久久久国产欧美日韩av| 男人操女人黄网站| 岛国在线观看网站| 一个人免费在线观看的高清视频 | av在线播放精品| 欧美av亚洲av综合av国产av| 水蜜桃什么品种好| 久久人人爽人人片av| 国产精品一区二区在线观看99| 国产精品一区二区免费欧美 | 大香蕉久久成人网| 精品高清国产在线一区| av欧美777| 国产欧美日韩综合在线一区二区| 在线十欧美十亚洲十日本专区| 美女高潮喷水抽搐中文字幕| 成人国语在线视频| 精品少妇一区二区三区视频日本电影| 女性生殖器流出的白浆| www.精华液| 999久久久精品免费观看国产| 大型av网站在线播放| 国产亚洲午夜精品一区二区久久| av电影中文网址| 久久毛片免费看一区二区三区| 美女大奶头黄色视频| 久久久久久久久免费视频了| 91成人精品电影| 两性夫妻黄色片| 亚洲一区中文字幕在线| 岛国毛片在线播放| 久久久久精品国产欧美久久久 | 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 建设人人有责人人尽责人人享有的| 国产福利在线免费观看视频| 男男h啪啪无遮挡| 亚洲久久久国产精品| 亚洲国产av影院在线观看| 婷婷成人精品国产| 亚洲精华国产精华精| 久久国产亚洲av麻豆专区| 亚洲视频免费观看视频| 18禁裸乳无遮挡动漫免费视频| 黄色 视频免费看| 亚洲激情五月婷婷啪啪| 亚洲自偷自拍图片 自拍| 90打野战视频偷拍视频| 窝窝影院91人妻| 午夜福利,免费看| 日本猛色少妇xxxxx猛交久久| 乱人伦中国视频| 一区二区三区精品91| 美女国产高潮福利片在线看| 久久女婷五月综合色啪小说| 青春草亚洲视频在线观看| av在线老鸭窝| 最黄视频免费看| 久久人人爽人人片av| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频| 亚洲视频免费观看视频| 少妇的丰满在线观看| 免费在线观看影片大全网站| 一级片免费观看大全| 在线观看免费高清a一片| 国产一卡二卡三卡精品| 大码成人一级视频| 午夜精品国产一区二区电影| 丝袜人妻中文字幕| 色婷婷久久久亚洲欧美| 成年人免费黄色播放视频| av一本久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩一级在线毛片| 午夜两性在线视频| 男女免费视频国产| 精品卡一卡二卡四卡免费| 高清视频免费观看一区二区| www.精华液| 亚洲国产毛片av蜜桃av| av天堂久久9| 亚洲精品国产av蜜桃| 欧美亚洲 丝袜 人妻 在线| 国产老妇伦熟女老妇高清| 国产欧美亚洲国产| 中国国产av一级| 国产日韩欧美亚洲二区| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| 狠狠狠狠99中文字幕| 欧美乱码精品一区二区三区| www.av在线官网国产| 久久久久久久精品精品| 亚洲,欧美精品.| 男女国产视频网站| av欧美777| 亚洲第一欧美日韩一区二区三区 | 国产精品一区二区在线不卡| 精品国产一区二区三区四区第35| 十分钟在线观看高清视频www| 国产在线一区二区三区精| 日韩制服骚丝袜av| 国产成人欧美| 国产精品成人在线| 一级黄色大片毛片| cao死你这个sao货| 亚洲自偷自拍图片 自拍| 90打野战视频偷拍视频| 亚洲 欧美一区二区三区| 午夜福利,免费看| 久久午夜综合久久蜜桃| 久久久精品免费免费高清| 大片电影免费在线观看免费| 午夜影院在线不卡| 久久久久久亚洲精品国产蜜桃av| 精品人妻1区二区| 视频区图区小说| av线在线观看网站| 日本撒尿小便嘘嘘汇集6| 国产在线一区二区三区精| 黄色片一级片一级黄色片| 在线观看www视频免费| 丝袜人妻中文字幕| 又紧又爽又黄一区二区| 亚洲av欧美aⅴ国产| 国产精品久久久久久精品电影小说| 中文字幕另类日韩欧美亚洲嫩草| 一本久久精品| 国产成人精品无人区| 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| 天天躁夜夜躁狠狠躁躁| 99久久国产精品久久久| 亚洲伊人色综图| 国产精品秋霞免费鲁丝片| 成人国产一区最新在线观看| 操出白浆在线播放| 在线亚洲精品国产二区图片欧美| 国产欧美日韩一区二区精品| 亚洲国产欧美日韩在线播放| 亚洲色图 男人天堂 中文字幕| 热re99久久精品国产66热6| 大香蕉久久成人网| 动漫黄色视频在线观看| 成人手机av| 婷婷色av中文字幕| 国产色视频综合| 另类精品久久| 99香蕉大伊视频| 嫩草影视91久久| 日本撒尿小便嘘嘘汇集6| 精品亚洲成国产av| 免费不卡黄色视频| 少妇精品久久久久久久| 亚洲第一欧美日韩一区二区三区 | 黄片播放在线免费| 一二三四社区在线视频社区8| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久久久99蜜臀| 精品国内亚洲2022精品成人 | 少妇人妻久久综合中文| 新久久久久国产一级毛片| 首页视频小说图片口味搜索|