• <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).

    日韩大码丰满熟妇| 99国产极品粉嫩在线观看| 中文字幕人妻丝袜一区二区| av超薄肉色丝袜交足视频| 国产精品秋霞免费鲁丝片| 91大片在线观看| 成年版毛片免费区| 国产又爽黄色视频| 狠狠婷婷综合久久久久久88av| 亚洲熟女精品中文字幕| 亚洲在线自拍视频| 老司机靠b影院| 久久中文看片网| 一级片免费观看大全| 日本vs欧美在线观看视频| 亚洲 欧美一区二区三区| 国产乱人伦免费视频| 国产xxxxx性猛交| 久久人妻福利社区极品人妻图片| 久久精品国产99精品国产亚洲性色 | 热re99久久精品国产66热6| 久久亚洲精品不卡| 中国美女看黄片| 国产激情久久老熟女| 欧美性长视频在线观看| 国产野战对白在线观看| 亚洲av日韩在线播放| 18禁裸乳无遮挡动漫免费视频| 精品国产国语对白av| 成人国产一区最新在线观看| 亚洲人成77777在线视频| 操出白浆在线播放| 王馨瑶露胸无遮挡在线观看| 免费黄频网站在线观看国产| 99国产精品一区二区蜜桃av | 91麻豆精品激情在线观看国产 | 欧美国产精品一级二级三级| 国产视频一区二区在线看| 国产av一区二区精品久久| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图| 一级a爱视频在线免费观看| avwww免费| 亚洲aⅴ乱码一区二区在线播放 | 午夜两性在线视频| 日韩三级视频一区二区三区| x7x7x7水蜜桃| 又黄又粗又硬又大视频| 人妻 亚洲 视频| 中亚洲国语对白在线视频| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| 中文字幕人妻丝袜一区二区| 国产麻豆69| 99国产精品99久久久久| 国产在线观看jvid| 色精品久久人妻99蜜桃| 久久人人爽av亚洲精品天堂| 12—13女人毛片做爰片一| 啦啦啦 在线观看视频| 麻豆国产av国片精品| 淫妇啪啪啪对白视频| 国产亚洲一区二区精品| 美女高潮到喷水免费观看| 日韩有码中文字幕| 咕卡用的链子| 国产男靠女视频免费网站| 黄色视频,在线免费观看| 久久这里只有精品19| 热re99久久国产66热| 婷婷精品国产亚洲av在线 | 国产精品99久久99久久久不卡| 99国产综合亚洲精品| 国产精品一区二区在线观看99| 国产黄色免费在线视频| 国产精品久久视频播放| 亚洲一码二码三码区别大吗| 亚洲五月婷婷丁香| 成人黄色视频免费在线看| 精品久久久久久电影网| 国产成人欧美在线观看 | 久热这里只有精品99| 亚洲av日韩精品久久久久久密| 夜夜爽天天搞| 国产极品粉嫩免费观看在线| 亚洲国产中文字幕在线视频| 国产aⅴ精品一区二区三区波| 91在线观看av| 国产99白浆流出| 国产成人系列免费观看| 99在线人妻在线中文字幕 | 欧美亚洲日本最大视频资源| 俄罗斯特黄特色一大片| √禁漫天堂资源中文www| 又紧又爽又黄一区二区| 天堂√8在线中文| 国产国语露脸激情在线看| 王馨瑶露胸无遮挡在线观看| 久热这里只有精品99| 18禁裸乳无遮挡免费网站照片 | 国产免费av片在线观看野外av| 日韩大码丰满熟妇| 91精品国产国语对白视频| 成人国产一区最新在线观看| 亚洲一码二码三码区别大吗| 黄片播放在线免费| 日韩制服丝袜自拍偷拍| 亚洲av熟女| 国产成人啪精品午夜网站| 久热爱精品视频在线9| 亚洲,欧美精品.| 男人操女人黄网站| 午夜福利影视在线免费观看| 日韩欧美一区二区三区在线观看 | 激情在线观看视频在线高清 | 真人做人爱边吃奶动态| 午夜福利视频在线观看免费| 男女下面插进去视频免费观看| 男女高潮啪啪啪动态图| 亚洲全国av大片| 国产精品av久久久久免费| 欧美一级毛片孕妇| 身体一侧抽搐| 日本a在线网址| 国产成人啪精品午夜网站| 国产亚洲精品久久久久久毛片 | 日韩免费高清中文字幕av| 久久狼人影院| av福利片在线| 欧美 亚洲 国产 日韩一| av在线播放免费不卡| 丝袜人妻中文字幕| 一级毛片高清免费大全| 在线观看免费视频网站a站| 黄片大片在线免费观看| 欧美乱码精品一区二区三区| 国产免费av片在线观看野外av| 中文欧美无线码| 黄色毛片三级朝国网站| 成年动漫av网址| 男人的好看免费观看在线视频 | 老熟女久久久| a级毛片在线看网站| 人人妻人人添人人爽欧美一区卜| 国产一区有黄有色的免费视频| 老鸭窝网址在线观看| 欧美久久黑人一区二区| 国产亚洲精品第一综合不卡| 欧美日韩成人在线一区二区| 亚洲专区字幕在线| 黄片小视频在线播放| 国产成人av教育| 欧美 亚洲 国产 日韩一| 国产免费av片在线观看野外av| 久久久国产成人免费| netflix在线观看网站| 久久狼人影院| 国产成人av教育| 欧美精品啪啪一区二区三区| 国产精品二区激情视频| 人妻久久中文字幕网| 亚洲三区欧美一区| 日韩欧美免费精品| 久久国产精品人妻蜜桃| 亚洲国产毛片av蜜桃av| 国产精华一区二区三区| 国产免费男女视频| 国产一卡二卡三卡精品| 欧美日韩福利视频一区二区| 欧美精品亚洲一区二区| 亚洲精品国产精品久久久不卡| 国产97色在线日韩免费| 国产无遮挡羞羞视频在线观看| 亚洲人成77777在线视频| 欧美另类亚洲清纯唯美| 999久久久国产精品视频| 狠狠婷婷综合久久久久久88av| 国产高清国产精品国产三级| 国产精品久久久av美女十八| 亚洲一卡2卡3卡4卡5卡精品中文| 精品福利观看| 黑丝袜美女国产一区| 欧美国产精品一级二级三级| videos熟女内射| 一级,二级,三级黄色视频| 亚洲 欧美一区二区三区| 日本欧美视频一区| 麻豆国产av国片精品| 免费在线观看视频国产中文字幕亚洲| 亚洲,欧美精品.| 国产亚洲欧美98| 久热爱精品视频在线9| 久久久久久久午夜电影 | 手机成人av网站| 十八禁网站免费在线| 国产精品亚洲av一区麻豆| 一边摸一边做爽爽视频免费| 亚洲av日韩在线播放| 看免费av毛片| 国产不卡av网站在线观看| 99久久国产精品久久久| 九色亚洲精品在线播放| 成人亚洲精品一区在线观看| 黑人欧美特级aaaaaa片| 妹子高潮喷水视频| 欧美老熟妇乱子伦牲交| 国产欧美日韩综合在线一区二区| 国产精品久久电影中文字幕 | 女人高潮潮喷娇喘18禁视频| 国产精品偷伦视频观看了| 久久精品国产综合久久久| av在线播放免费不卡| 欧美日韩视频精品一区| 777久久人妻少妇嫩草av网站| 岛国毛片在线播放| 久久久久精品国产欧美久久久| 成人国产一区最新在线观看| 欧美成狂野欧美在线观看| 欧美乱色亚洲激情| 久久香蕉国产精品| 久久人人爽av亚洲精品天堂| 成年动漫av网址| 久9热在线精品视频| 久久精品亚洲精品国产色婷小说| 大香蕉久久网| 丰满迷人的少妇在线观看| 日韩成人在线观看一区二区三区| 99热国产这里只有精品6| 激情在线观看视频在线高清 | 激情在线观看视频在线高清 | 老司机在亚洲福利影院| 国产有黄有色有爽视频| 亚洲欧美日韩高清在线视频| 一区二区日韩欧美中文字幕| 精品国产一区二区久久| 悠悠久久av| 一级片免费观看大全| 精品卡一卡二卡四卡免费| 欧美激情 高清一区二区三区| 在线av久久热| 国产高清视频在线播放一区| 亚洲av欧美aⅴ国产| 亚洲中文日韩欧美视频| 亚洲久久久国产精品| av网站在线播放免费| 久久久久久久国产电影| 青草久久国产| 下体分泌物呈黄色| 51午夜福利影视在线观看| 啦啦啦视频在线资源免费观看| 久热这里只有精品99| 操美女的视频在线观看| 999精品在线视频| 交换朋友夫妻互换小说| 久久精品亚洲熟妇少妇任你| 国产精品av久久久久免费| 人妻丰满熟妇av一区二区三区 | 日本vs欧美在线观看视频| 精品国产超薄肉色丝袜足j| 满18在线观看网站| 日韩免费高清中文字幕av| 久久久水蜜桃国产精品网| 国产成人啪精品午夜网站| 新久久久久国产一级毛片| 香蕉国产在线看| 9热在线视频观看99| 90打野战视频偷拍视频| 91老司机精品| 国产一卡二卡三卡精品| 国内久久婷婷六月综合欲色啪| 91成人精品电影| 成人黄色视频免费在线看| av免费在线观看网站| 极品教师在线免费播放| 日日摸夜夜添夜夜添小说| 亚洲片人在线观看| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 最新的欧美精品一区二区| 欧美性长视频在线观看| 久久精品人人爽人人爽视色| 窝窝影院91人妻| 一进一出好大好爽视频| 亚洲国产中文字幕在线视频| 国产成人系列免费观看| 女人被狂操c到高潮| 亚洲专区中文字幕在线| 欧美日韩福利视频一区二区| 999精品在线视频| 久久精品91无色码中文字幕| 97人妻天天添夜夜摸| 免费女性裸体啪啪无遮挡网站| 精品高清国产在线一区| 国产欧美亚洲国产| 国产成人av激情在线播放| 日本欧美视频一区| 精品欧美一区二区三区在线| 国产精品久久久av美女十八| 中文字幕色久视频| tube8黄色片| 女警被强在线播放| 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| av免费在线观看网站| 精品免费久久久久久久清纯 | 久久天躁狠狠躁夜夜2o2o| 不卡av一区二区三区| 日韩视频一区二区在线观看| 亚洲国产欧美网| 国产欧美亚洲国产| 丁香六月欧美| 日本精品一区二区三区蜜桃| 精品一区二区三区av网在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产色婷婷电影| 国产又色又爽无遮挡免费看| 老司机福利观看| 男女午夜视频在线观看| 国产国语露脸激情在线看| 免费在线观看影片大全网站| 欧美人与性动交α欧美精品济南到| 亚洲熟女毛片儿| 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 99国产精品一区二区蜜桃av | 免费人成视频x8x8入口观看| 一级毛片女人18水好多| 好男人电影高清在线观看| 国产精品一区二区在线观看99| 法律面前人人平等表现在哪些方面| 999精品在线视频| 香蕉丝袜av| 国产日韩一区二区三区精品不卡| 精品亚洲成国产av| 首页视频小说图片口味搜索| 韩国精品一区二区三区| 极品人妻少妇av视频| 纯流量卡能插随身wifi吗| 国产欧美日韩一区二区精品| 悠悠久久av| 成人黄色视频免费在线看| 动漫黄色视频在线观看| 香蕉国产在线看| 热re99久久精品国产66热6| 超色免费av| 亚洲欧美日韩另类电影网站| 一边摸一边抽搐一进一小说 | 免费在线观看黄色视频的| 悠悠久久av| 后天国语完整版免费观看| 国产亚洲精品久久久久久毛片 | 日韩欧美国产一区二区入口| 性少妇av在线| 91成人精品电影| 久久久久久亚洲精品国产蜜桃av| 波多野结衣av一区二区av| 亚洲成a人片在线一区二区| 午夜免费成人在线视频| 欧美黄色淫秽网站| 亚洲国产精品sss在线观看 | 成人18禁高潮啪啪吃奶动态图| 极品少妇高潮喷水抽搐| 香蕉丝袜av| 天天添夜夜摸| 99riav亚洲国产免费| 午夜成年电影在线免费观看| 欧美日韩黄片免| 男女免费视频国产| netflix在线观看网站| 亚洲久久久国产精品| 夫妻午夜视频| 嫩草影视91久久| 欧美日韩乱码在线| 黄色丝袜av网址大全| 成人av一区二区三区在线看| 天天操日日干夜夜撸| 老司机深夜福利视频在线观看| 亚洲人成77777在线视频| 日日爽夜夜爽网站| 深夜精品福利| 黑人猛操日本美女一级片| 欧美在线黄色| 亚洲一区高清亚洲精品| 午夜视频精品福利| 亚洲黑人精品在线| 在线观看免费高清a一片| 日韩欧美三级三区| 欧洲精品卡2卡3卡4卡5卡区| 国产xxxxx性猛交| 黄色片一级片一级黄色片| 久久精品亚洲熟妇少妇任你| 在线免费观看的www视频| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 天天添夜夜摸| 男女高潮啪啪啪动态图| 午夜影院日韩av| 如日韩欧美国产精品一区二区三区| 一本综合久久免费| 深夜精品福利| 国产激情久久老熟女| 婷婷精品国产亚洲av在线 | 精品国产乱子伦一区二区三区| 免费少妇av软件| 成人黄色视频免费在线看| 欧美日韩av久久| 亚洲中文av在线| av一本久久久久| 亚洲欧美色中文字幕在线| 久久香蕉激情| 侵犯人妻中文字幕一二三四区| 免费观看人在逋| 国产成人免费无遮挡视频| 日韩三级视频一区二区三区| 91在线观看av| 久久久久久久精品吃奶| 久热爱精品视频在线9| 大香蕉久久成人网| 18禁观看日本| 久久久久久亚洲精品国产蜜桃av| 亚洲片人在线观看| 精品一区二区三区av网在线观看| 90打野战视频偷拍视频| 757午夜福利合集在线观看| a级毛片在线看网站| 建设人人有责人人尽责人人享有的| 久9热在线精品视频| 成人国产一区最新在线观看| 久久草成人影院| 日本欧美视频一区| 久久午夜综合久久蜜桃| 美女高潮到喷水免费观看| 欧美乱码精品一区二区三区| 夜夜爽天天搞| 999精品在线视频| 精品无人区乱码1区二区| 久久国产乱子伦精品免费另类| 亚洲精品国产一区二区精华液| 午夜免费观看网址| 久久 成人 亚洲| 香蕉丝袜av| 成年版毛片免费区| 欧美日韩视频精品一区| 精品亚洲成a人片在线观看| 国产亚洲欧美精品永久| 国产精品免费一区二区三区在线 | 在线永久观看黄色视频| 人人妻人人添人人爽欧美一区卜| 大陆偷拍与自拍| 免费日韩欧美在线观看| 99riav亚洲国产免费| 久久久久久人人人人人| 十分钟在线观看高清视频www| 一区福利在线观看| 亚洲欧美日韩高清在线视频| 久久中文字幕人妻熟女| 亚洲专区国产一区二区| 香蕉久久夜色| 国产成人精品无人区| 国产精品一区二区精品视频观看| 欧美人与性动交α欧美精品济南到| 精品人妻1区二区| 99久久综合精品五月天人人| 欧美老熟妇乱子伦牲交| 午夜日韩欧美国产| 国产成人精品久久二区二区91| 女同久久另类99精品国产91| 视频区图区小说| 巨乳人妻的诱惑在线观看| 国产高清国产精品国产三级| 国产一区二区激情短视频| 日韩免费高清中文字幕av| 亚洲国产欧美一区二区综合| 丝袜在线中文字幕| 国产单亲对白刺激| 亚洲,欧美精品.| 热re99久久国产66热| 亚洲av熟女| 激情视频va一区二区三区| 精品无人区乱码1区二区| 久久狼人影院| 精品国内亚洲2022精品成人 | 热re99久久国产66热| 男女免费视频国产| 亚洲精品国产区一区二| 国产精品亚洲av一区麻豆| 午夜福利一区二区在线看| 国产精品电影一区二区三区 | 激情视频va一区二区三区| 久久这里只有精品19| 一a级毛片在线观看| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 精品福利观看| 操出白浆在线播放| av在线播放免费不卡| 黄色怎么调成土黄色| 欧美激情高清一区二区三区| 日韩中文字幕欧美一区二区| 两人在一起打扑克的视频| 少妇被粗大的猛进出69影院| 一级毛片精品| 交换朋友夫妻互换小说| 午夜视频精品福利| 中出人妻视频一区二区| www.999成人在线观看| 亚洲专区字幕在线| 国产精华一区二区三区| 亚洲成人免费av在线播放| 国产精品一区二区在线不卡| 午夜福利在线免费观看网站| 国产精品免费视频内射| 国产麻豆69| 两个人免费观看高清视频| 国产1区2区3区精品| 亚洲av日韩在线播放| 不卡一级毛片| 成年女人毛片免费观看观看9 | 国产激情久久老熟女| 在线天堂中文资源库| 老司机深夜福利视频在线观看| 在线观看免费日韩欧美大片| 久久国产乱子伦精品免费另类| 最新美女视频免费是黄的| 亚洲性夜色夜夜综合| 国产亚洲av高清不卡| 日韩一卡2卡3卡4卡2021年| 久久精品国产亚洲av香蕉五月 | av网站免费在线观看视频| 日韩免费av在线播放| 免费久久久久久久精品成人欧美视频| 亚洲精品粉嫩美女一区| 亚洲色图av天堂| 日本一区二区免费在线视频| 亚洲片人在线观看| 热re99久久国产66热| 亚洲一区二区三区欧美精品| 极品少妇高潮喷水抽搐| xxx96com| 亚洲欧美一区二区三区黑人| 国产成人欧美| 美女扒开内裤让男人捅视频| 国产一区二区激情短视频| 精品福利观看| x7x7x7水蜜桃| 在线免费观看的www视频| 成人永久免费在线观看视频| 91国产中文字幕| 精品久久久久久,| 我的亚洲天堂| 亚洲aⅴ乱码一区二区在线播放 | 久久中文字幕一级| 国产无遮挡羞羞视频在线观看| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩精品亚洲av| 亚洲一区中文字幕在线| 人人妻人人澡人人看| 一进一出好大好爽视频| 午夜福利影视在线免费观看| 亚洲色图 男人天堂 中文字幕| 美女午夜性视频免费| 免费观看精品视频网站| 天天操日日干夜夜撸| www.自偷自拍.com| 正在播放国产对白刺激| 欧美激情高清一区二区三区| 美女视频免费永久观看网站| www.精华液| 黑丝袜美女国产一区| 在线观看www视频免费| 婷婷丁香在线五月| 女人久久www免费人成看片| 在线观看66精品国产| av片东京热男人的天堂| 99热只有精品国产| 老司机亚洲免费影院| 视频区图区小说| 日韩欧美一区二区三区在线观看 | 亚洲精品国产精品久久久不卡| 国产精品98久久久久久宅男小说| av片东京热男人的天堂| 一二三四社区在线视频社区8| 欧美不卡视频在线免费观看 | 欧美 亚洲 国产 日韩一| 一级作爱视频免费观看| 国产蜜桃级精品一区二区三区 | 99riav亚洲国产免费| 欧美av亚洲av综合av国产av| 丁香六月欧美| a在线观看视频网站| 在线视频色国产色| 久久青草综合色| 久久精品91无色码中文字幕| 午夜福利视频在线观看免费| 欧美最黄视频在线播放免费 | 国产精品秋霞免费鲁丝片| 欧美激情 高清一区二区三区| 在线十欧美十亚洲十日本专区| 人人妻人人澡人人看| 免费在线观看完整版高清| 欧美激情久久久久久爽电影 | 成人黄色视频免费在线看| 一区二区日韩欧美中文字幕| 国产亚洲精品一区二区www | 在线观看一区二区三区激情| 日日夜夜操网爽| 又大又爽又粗| 国产xxxxx性猛交| 丝袜美足系列| 久久性视频一级片| av在线播放免费不卡| 韩国av一区二区三区四区| 久9热在线精品视频| 国产欧美日韩一区二区三| 国产精品久久久久成人av|