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

    Aircraft Measurements of the Microphysical Properties of Stratiform Clouds with Embedded Convection

    2021-06-04 08:46:50TuanjieHOUHengchiLEIYoujiangHEJiefanYANGZhenZHAOandZhaoxiaHU
    Advances in Atmospheric Sciences 2021年6期

    Tuanjie HOU, Hengchi LEI, Youjiang HE, Jiefan YANG, Zhen ZHAO, and Zhaoxia HU

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

    2Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science and Technology, Nanjing 210044, China

    3Institute of Atmospheric Environment, Chinese Research Academy of Environmental Sciences, Beijing 100012, China

    ABSTRACT

    Key words:aircraft measurements, liquid water content, ice particle habit, particle size distributions

    1.Introduction

    Stratiform clouds associated with frontal systems under large-scale synoptic conditions frequently occur over northern China. These clouds are commonly multilayered mixedphase clouds that persist from several hours to two days.The precipitating stratiform clouds in He-Tao cyclones(Wang et al., 1982) generally consist of upper altostratus and lower stratocumulus, with cloud top temperatures varying from ?8°C to ?24°C. The winter frontal clouds (You and Liu, 1995) in the northwestern Xinjiang Uygur autonomous region often consist of below-front stratocumulus,above-front altostratus, and cloud bands in the upper atmosphere. Some stratiform clouds associated with cold frontal systems begin with the formation of upper-level cirrus above 6 km, and postfrontal precipitation is produced from altostratus and stratocumulus clouds (Hou et al., 2016). An evident dry layer, with a thickness of 1 km between the upper altostratus and lower stratocumulus clouds, was observed in the dissipating stage of one stratiform cloud case in Jilin Province (Hou et al., 2010).

    To reveal the microphysical characteristics of stratiform clouds, in situ aircraft measurements have been conducted in recent years in China (Ma et al., 2007; Hou et al.,2014; Guo et al., 2015; Yang and Lei, 2016). High ice particle concentrations with a magnitude of 10Lwere observed at temperatures between ?3°C and ?6°C in nimbostratus clouds (Zhao and Lei, 2014; Yang et al., 2017),which could be explained by the Hallett-Mossop ice multiplication process (Hallett and Mossop, 1974). Yang et al.(2014) showed that high ice particle concentrations occur due to breakup dendrite crystal-crystal collisions in the region of ?10°C to ?15°C in stratiform clouds. The active riming and aggregation processes at temperatures warmer than?9.5°C (Zhu et al., 2015) were responsible for the broader particle size distributions within embedded convective regions as opposed to stratiform regions. These previous studies have provided insights into ice particle size distributions and growth modes, but much remains to be learned. The physical mechanisms that control supercooled droplet size distributions and secondary ice particle processes are not well understood. The distribution of cloud droplets influences the development of cloud microphysics via conversion to ice particles. Knowledge of supercooled water is therefore important for evaluating the optimum temperatures for riming occurrence and the potential of clouds for precipitation enhancement.

    One important feature of stratiform clouds over northern China is the existence of embedded convective cells with enhanced reflectivity. The radar returns within the convective region are approximately 10 to 20 dBZ larger than those observed in the surrounding stratiform region (Zhu et al., 2015; Yang et al., 2017), indicating either high ice particle concentrations or a heavy degree of riming.However, few observations are available variations in the particle size distributions in the vertical direction and within different regions of clouds over northern China. Field(1999) studied the evolution of ice particle size distributions based on observations of Lagrangian spiral descents within an altostratus cloud and found that aggregation dominated particle growth at temperatures between ?10°C and?20°C. Several previous studies (e.g., Herzegh and Hobbs,1980; Stith et al., 2002; Lawson and Zuidema, 2009; Crosier et al., 2011) have supported enhanced ice particle growth in the embedded convective region. Herzegh and Hobbs (1980) found that riming was important in the convective cells of warm frontal clouds where the updrafts reached 60 cm s. The embedded convection in the frontal Arctic cloud system (Lawson and Zuidema, 2009) produced more graupel particles and larger snowflakes. The ice concentrations in the convective region of the stratus (Crosier et al.,2011) were more than one order of magnitude greater than the concentration of primary ice in the overlaying stratus.The embedded convection significantly affected the ice microphysical properties and precipitation development.However, the differences in the supercooled water content and particle size distributions between the convective and stratiform regions are not well understood.

    The Taihang Mountains Precipitation Enhancement and Hail Suppression Experiment, which is an observational field program, was conducted to the east of the Taihang Mountains during 2017 and 2019. The experiment focused on observations of stratiform clouds with embedded convection and hailstorms to provide insights into precipitation formation processes and to contribute to water resource management. For stratiform clouds, aircraft and ground observations were performed to reveal the microphysical characteristics and subsequently assess the potential for weather modification.

    In this paper, we present aircraft measurements of a precipitating stratiform cloud case over northern China on 22 May 2017. The King-350 aircraft conducted spiral ascents and descents at different regions of the cloud system. The aircraft flew to a height of 7.26 km, which is higher than the previous measurement height of 6 km from the Y-12E and Xiayan 3A aircraft. Microphysical parameters including supercooled liquid water content, particle concentrations, habits,and combined size distributions, were used to demonstrate the differences between the convective and stratiform cloud regions.

    The primary goal of this study is to learn more about how supercooled water and particle size distributions vary from upper to lower levels and from convective to stratiform regions. The cloud droplet probe (CDP) and the King probe were used together in this study for quantitative liquid water measurements and cloud phase identification.Through a comparison of particle imagery and associated concentrations, we aim to investigate the temperature region where riming and the associated rime splintering occur.Rime splintering, known as the Hallett-Mossop process (Hallett and Mossop, 1974), has been used to explain the high ice particle concentrations in some clouds. The Hallett-Mossop effect occurs between approximately ?3°C and ?8°C in the presence of supercooled droplets smaller than 13 μm and larger than 25 μm in diameter. We present small and large particle concentration variations in the vertical combined with particle habits to discuss the broader particle size distributions observed in convective regions compared to those observed in stratiform regions.

    The remainder of this paper is organized as follows. In section 2, the instrumentation and data processing are presented. An overview of the synoptic features, radar information, and flight paths are presented in section 3. The results,including the distribution of the supercooled liquid water content (LWC), particle size distributions, ice particle habits,and growth processes are analyzed in sections 4 and 5.Finally, the study is summarized and conclusions are presented in section 6.(King et al., 1978), incomplete response to the large end of the cloud droplet spectrum (Feind et al., 2000), and is subject to contamination by ice particles in mixed-phase conditions. The presence of ice crystals generally biases the hotwire LWC by less than 0.025 g m, but the difference can be as much as 0.1 g m(Cober et al., 1995).

    Considering the aforementioned uncertainties in King LWCs, LWC measurements from the CDP were also presented in this study. The CDP provided particle concentration measurements for sizes smaller than 50 μm and could be used to derive liquid water content. Since the CDP has difficulty distinguishing small ice particles from droplets during measurements, the droplet spectra measured by the CDP were inevitably contaminated by ice particles for mixedphase conditions. It also had uncertainty in sizing droplets;therefore, these biases should be considered when analyzing the CDP estimated liquid water content. Comparisons of LWCs between the CDP and a Nevzorov hotwire probe(Faber et al., 2018) suggested that the CDP-estimated LWC exceeded that measured by the Nevzorov by approximately 20%. The CDP and hotwire have been used together to estimate LWC in previous studies (Heymsfield and Willis, 2014;Wang et al., 2015; Giangrande et al., 2016).

    Figure 1 shows a comparison of King LWCs and CDPestimated LWCs. The CDP-estimated LWC was broadly consistent with that measured by the King probe. The best linear fit was LWC= 0.02 + 0.82 LWCwith a correlation coefficient of 0.87.

    The CDP and King probe were used together in this study to identify the existence of supercooled liquid water for discriminating between mixed-phase and ice clouds. We only analyzed clouds colder than 0°C in the study, not referring to liquid-only cloud regions. Ice particles were observed in all the flight legs listed in the study; therefore,the criteria to identify the mixed-phase and ice clouds did not include particle concentrations from the 2D-S. The criteria are listed in Table 2. The minimum CDP concentration of 10 cmand the King LWC of 0.1 g mwere chosen to identify mixed-phase clouds following Lance et al. (2010). When one of the above two criteria was met, we referred to the cloud as a mixed-phase cloud. When the CDP concentration was between 1 and 10 cmand the King-LWC was less than 0.1 g m, we could not state with confidence whether supercooled liquid water existed in the cloud. Although the CDP was equipped with anti-shattering tips, the probe could still be contaminated by shattering artifacts in the presence of ice. Ice particles may cause artifi-

    2.Instrumentation and data processing

    This study focused on a typical stratiform rainfall event over northern China on 22 May 2017. The airborne data for this case were collected using a King-350 aircraft (Yang et al., 2019) from the Weather Modification Office of Hebei Province, China. The instruments used in the study (Table 1)included a King hotwire probe (LWC-300), a CDP, a 2D-S optical array probe, a high-volume precipitation spectrometer (HVPS), and a cloud particle imager (CPI) probe.

    The CDP, which is manufactured by Droplet Measurement Technologies, Inc. (DMT), measures particles with a maximum dimension (D) between 2 and 50 μm. It has a bin width of 1 μm for particles smaller than 14 μm and a bin width of 2 μm for particles larger than 14 μm. The CDP features an open-path laser that is passed between two arms with anti-shattering tips; therefore, ice crystal shattering does not substantially affect the CDP size distribution measurements (Lance et al., 2010).

    The 2D-S has a size range of 10 μm to 3 mm, and the HVPS has a size range of 150 μm to 2 cm. The 2D-S consists of two independent 128 photodiode arrays that are occulted when a particle passes through the sample volume. The averaged data from the two channels were used in this study. To reduce the effects of large ice particle shattering,2D-S probe tips are designed to reduce the number of shattered particles reaching the probe sample volume and apply postprocessing techniques that use an interarrival time algorithm to remove shattered particles (Lawson, 2011;Lawson et al., 2015).

    The CPI, which has a resolution of 2.3 μm, provided high-resolution cloud particle images for particles in the size range from 10 to 2000 μm. The probe was used to assess ice particle habits and to identify the occurrence of riming.

    2.1.LWC measurements and phase determination

    The King probe is a widely used probe for measuring liquid water. The LWC in this study refers to that measured by the King probe unless otherwise specified. However, the King hotwire also has limitations in the LWC measurements, including uncertainty in collecting small droplets < 5 μmcial counts of up to 1 cmor higher in pure ice clouds. For strict exclusion of supercooled liquid water from ice clouds,clouds were only identified as an ice cloud region when the CDP concentration was less than 1 cmand the King LWC was less than 0.1 g m. In this way, we attempted to eliminate the contamination of supercooled liquid water when calculating the ice water content (IWC) from the observed ice particle concentrations.

    Table 1. Overview of the airborne instruments.

    Fig. 1. Comparison of the King LWCs and CDP-estimated LWCs. Each data point represents a 5-s time average. The red line is the linear best fit line with a correlation coefficient of 0.87. .

    Table 2. The criteria to identify mixed-phase and ice clouds.

    2.2.IWC calculation

    For IWC measurements, the aircraft in Hebei was equipped with a Nevzorov hotwire TWC/LWC probe.However, due to the problem of a baseline drift for the Nevzorov TWC, its accuracy still needs testing. Therefore,the Nevzorov probe was not used in this study. Instead, we calculated the IWC from the particle size distribution measurements.

    The particle mass is usually expressed in the form m =a D, where D is taken to be a representation of the particle size and a and b are empirically derived terms. Different sets of parameters have been derived from previous studies(Locatelli and Hobbs, 1974; Mitchell et al., 1990; Brown and Francis, 1995; Baker and Lawson, 2006; Heymsfield et al., 2010; Wang et al., 2015). The relationship from Brown and Francis (1995) has been widely used in IWC calculations. However, their IWC estimates were based on observations primarily in the temperature range from ?20°C to?30°C, with particle sizes between 200 and 800 μm. The relationship from Baker and Lawson (2006) was developed based on the dataset from Mitchell et al. (1990). This dataset was collected at the surface and was not necessarily representative of the various ice particles inside the clouds. In comparison, the relationship from Heymsfield et al. (2010) was from the more recent observations of clouds with broader temperatures ranging from 0°C to ?60°C, and the particle sizes were between 100 and 2000 μm. Therefore, we used the mass-dimensional relationship of Heymsfield et al. (2010)to calculate the ice water content. The relationship is of the form m(D) = 0.00528 D, where m represents the mass of an ice particle in grams and D is the maximum particle dimension in centimeters, which is valid for both stratiform and convective clouds over a wide range of temperatures and particle sizes. Note that localized high ice particle concentrations may lead to deviations, and the coefficient a is higher by approximately a factor of 2 for convectively generated versus stratiform ice clouds.

    The particle size distributions were derived from a combination of the CDP, 2D-S, and HVPS measurements. The combined particle size distribution was then used to derive the ice water content. Figure 2 shows example size distributions from the three instruments at three different temperatures. According to the figure, the overlap between 10 and 30 μm in the examples allows for the combination of the CDP and 2D-S. For smoothing of the combined size distribution, the breakpoint between the CDP and 2D-S was chosen as 25 μm. and the breakpoint between the 2D-S and HVPS was 1105 μm. The combined particle size distribution included particles with maximum sizes between 2 and 24 μm,25 and 1105 μm and 1125 and 2 cm, from the CDP, 2D-S,and HVPS, respectively. The concentration data from the CDP, 2D-S, and HVPS were averaged over 5-s intervals to reduce spurious variability.

    3.Synoptic overview and flight paths

    On 22 May 2017, a cold frontal system approached northern China. The 500 hPa analysis at 0000 UTC showed a trough to the northwest of Hebei Province (Fig. 3). The study region, located in southern Hebei Province, was characterized by winds primarily from the west and southwest at 500 hPa. As the trough moved eastward, it deepened and produced wind shifts from northerly to southwesterly at the middle and lower levels. With the passage of the strengthening frontal system with cold-air advection, large-scale precipitation occurred over central and northern China.

    Figure 4 presents the observed 4-h accumulated precipitation from 0300 to 1500 UTC (valid at 0700, 1100, and 1500 UTC 22 May 2017) from the rain gauges. Figure 4a shows a southwest-northeast oriented band of accumulated rainfall over the central and northern provinces, with the precipitation center to the southwest of Hebei Province. For the study region, light rainfall with a peak value of 52.7 mm was found in the northwestern part, indicating the initiation of precipitation. In comparison, more rainfall was observed during the process of southward movement between 0700 and 1100 UTC (Fig. 4b). Rainfall covered most of the study region, with a higher maximum value of 81.4 mm. By 1500 UTC (Fig. 4c), the precipitation center had already moved out of Hebei Province. The 4-h accumulated precipitation in the study region decreased to generally less than 25 mm.Therefore, a relatively higher precipitation center with 4-h rainfall values larger than 25 mm was primarily located in the southwestern part of the study region from 0700 to 1100 UTC.

    Fig. 2. Sample size distributions showing the overlap among the CDP, 2D-S, and HVPS at (a) 0742:21 UTC, (b)0746:46 UTC, and (c) 0808:56 UTC. .

    Fig. 3. Geopotential height (units: m) and winds (units: m s?1)for 500 hPa at 0000 UTC 22 May 2017 from the National Centers for Environmental Prediction reanalysis data. The red square indicates the study region. .

    Considering the synoptic setting and precipitation formation, aircraft measurements were conducted in the study region from 0703 to 1032 UTC. The flying height did not exceed 7.26 km due to the restrictions of the aircraft. By the time of aircraft measurement, the cloud-top height was approximately 8 km, and the 0°C layer was at approximately 3.8 to 4.0 km. The aircraft departed from Zhengding airport, which was very close to the Shijiazhuang radar station(indicated with a red dot in Fig. 4, longitude: 114.7°E; latitude: 38.3°N; 134.8 m above MSL). Then, the aircraft flew to the southwest of the radar center, where evident precipitation was observed. Figure 5 shows the aircraft track between 0703 and 1022 UTC. The aircraft flight pattern first consisted of a gradual ascent from 2 km to 7.26 km between 0737 and 0753 UTC. After that, the aircraft continued to fly southwest toward the embedded convective cells at 7.26 km. The aircraft then descended to 2.1 km between 0808 and 0823 UTC and ascended from 2.1 to 7.26 km between 0831 and 0848 UTC. The aircraft then made another descent, ascent, and horizontal penetration at 7.26 km in the process of flying back to the airport. As the study focused on ice particle characteristics, the data at temperatures colder than 0°C (at approximately 4 km) were divided into six flight legs with colored traces, as shown in Fig. 5.The flight leg features, including the beginning and end times, altitudes, and cloud temperatures are listed in Table 3.The flight paths in the present study included three upward spirals, two downward spirals, and one constant-altitude penetration at 7.26 km.

    To clearly illustrate the cloud development and flight track, Figure 6 presents the reflectivity fields from the Shijiazhuang radar observations with an overlay of the flight track at the ending times of the four legs. At 0754 UTC (Fig. 6a), a southwest-northeast oriented cloud region with a general reflectivity of 20-30 dBZ was observed. The scattered region to the southwest of the radar center had higher values of 35-40 dBZ, indicating the existence of convection within the stratiform cloud system. The first spiral ascent (leg 1) from 4.2 to 7.2 km between 0741 and 0753 UTC was at the echo region of approximately 20-30 dBZ.According to the rain gauge observations, precipitation started at 0700 UTC below the aircraft measurement region.The 1-h accumulated precipitation during 0700 and 0800 UTC below leg 1 was 7.6 mm.

    Fig. 4. Observed 4-h accumulated precipitation (units: mm) from (a) 0300 to 0700 UTC, (b) 0700 to 1100 UTC and(c) 1100 to 1500 UTC 22 May 2017. The red square indicates the study region, and the red dot represents the location of the radar site at Shijiazhuang. .

    Fig. 5. Aircraft track (black trace) starting from S between 0703 and 1022 UTC. The six flight legs shown in Table 2 are highlighted in colored traces. .

    Table 3. Flight leg features, including beginning and end times, altitudes, and cloud temperatures.

    The cloud system then moved south. The most evident change was the gradual organization of the convective region in the southwestern part, from a few scattered small regions to one large echo region at 0848 (Fig. 6b) with reflectivity values between 30 and 40 dBZ. Therefore, the first spiral descent (leg 2) between 0808 and 0818 UTC was at one of the small convective regions, while the second spiral ascent (leg 3) between 0836 and 0848 UTC was conducted at the edge of the large convective region. The 1-h accumulated precipitation during 0800 and 0900 UTC, below legs 2 and 3, was 10.2 mm and 6.8 mm, respectively.

    By 0930 UTC (Fig. 6c), the southwestern convective region had moved south and the region with values between 35 and 40 dBZ was much smaller, indicating a weakening trend. The horizontal penetration at 7.36 km between 0848 and 0908 UTC was in the region with values of approximately 30 dBZ.

    Fig. 6. PPI radar display of the radar reflectivity factor (Units: dBZ) from the Shijiazhuang radar 0.7° elevation with an overlay of the flight track (red traces starting from S) between (a) 0703 and 0753 UTC, (b) 0703 and 0934 UTC,(c) 0703 and 0848 UTC, and (d) 0703 and 1004 UTC.

    As the aircraft flew back to the airport, the second spiral descent (leg 5) between 0927 and 0934 was in the more homogeneous stratiform region with reflectivity values of approximately 25 dBZ. Similarly, the third spiral ascent (leg 6) was also in the stratiform region with values of 20 dBZ (Fig. 6d). The 1-h accumulated precipitation during 0900 and 1000 UTC, below the legs 5 and 6, was 6.1 mm and 0.6 mm, respectively.

    Therefore, the five spirals including legs 1, 5, and 6 were in the more homogeneous stratiform regions, while the spirals on legs 2 and 3 were in one embedded convective cell and at the edge of the convective region, respectively.The hourly precipitation rates below the convective region were higher, with values up to 10 mm h. In contrast, relatively lower precipitation rates were observed below the other cloud regions. The precipitation rates below leg 6 were the lowest, at less than 1 mm h. Thus, the cloud system was in an active state during the measurements.

    4.Supercooled liquid water

    The distribution of supercooled liquid water content within mixed-phase clouds can be complex. The evolution of the supercooled liquid water was closely related to the cloud developmental stages. The presence of supercooled liquid water provided conditions for the growth of nucleated particles, and it was an important parameter to determine whether clouds were suitable for seeding.

    Figure 7 shows the vertical distribution of the CDP-measured particle concentrations, the CDP LWCs, and the King LWCs on legs 1, 2, and 3. On flight leg 1 (Fig. 7a), the droplets were primarily measured at temperatures warmer than ?11°C, with concentrations peaking at 116 cm. The droplets at this stage were relatively small since the concentrations of droplets larger than 24 μm were less than 0.5 cm.Therefore, the derived CDP LWCs (Fig. 7b) were extremely low, at no more than 0.04 g m. Similarly, the King probe(Fig. 7c) also measured a very low LWC with a peak of 0.03 g m. Hence, the ice region was determined to be in the temperature regions colder than ?11.4°C, where the CDP concentrations were less than 0.25 cmand both the CDP and King LWCs were less than 0.01 g m.

    Fig. 7. Vertical distribution of the (a) CDP concentrations, (b) CDP LWC and (c) King LWC on leg 1; (d) CDP concentrations, (e) CDP LWC and (f) King LWC on leg 2; (g) CDP concentrations, (h) CDP LWC and (i) King LWC on leg 3. The red solid horizontal lines represent the transition to ice clouds. The regions between the dashed and solid lines represent the unidentified clouds. .

    Notably, at temperatures between ?6.6°C and ?8.5°C,the CDP-measured cloud particle concentrations were extremely low, at less than 0.5 cm. No supercooled liquid water was detected from either the CDP or the King probe.According to our criteria for measuring LWCs, that region also belonged to the ice region. The 2D-S observed ice crystal concentrations in that region varied from 8 to 20 Lwith mean volume diameters from 476 to 1100 μm, and the large particles were primarily dendritic and broad-branched crystals. That indicated vigorous growth of ice particles.The well-distinguished crystal branches suggested that vapor deposition was the dominant growth process. Generally, dendrites and sector plates grow under water-saturated conditions. Therefore, it should have been a mixed-phase region with extremely low liquid water contents at an earlier stage and the ice particle growth depleted cloud droplets.

    On leg 2 (Fig. 7d), the CDP-measured total droplet concentration at temperatures warmer than ?6°C was typically greater than 40 cm, peaking at 231 cm. Large droplets with diameters larger than 24 μm had concentrations that varied from 1.4 to 14 cmand were found within the temperatures ranging from ?4.4°C to ?4.9°C. A comparison of LWCs (Figs. 7e and 7f) from the CDP and King probes suggested that both probes detected evident supercooled cloud water at temperatures between ?4°C and ?5°C, with peaks of 0.45 and 0.57 g mat the same temperature of ?4°C.The two probes were consistent at temperatures warmer than ?5°C. However, the CDP detected virtually no liquid water at temperatures colder than ?6°C since small particle concentrations were less than 0.4 cm. In contrast, the King probe still measured LWCs varying from 0.05 to 0.2 g m,which was probably a baseline drift due to contamination from ice crystals. Previous research suggested that it is difficult for the King probe to accurately measure small LWCs in mixed-phase conditions. To minimize the errors caused by the baseline drift of the King probe, Cober et al. (1995)artificially zeroed the results from the King probe when the FSSP concentrations were less than 6 cm.

    CPI images suggested the existence of rimed ice particles at temperatures between ?10°C and ?1°C (Fig. 8).Unfortunately, most of the CPI images taken between 0808:31 and 0812:20 UTC (?15.7°C to ?9.9°C) on leg 2 were out of focus. According to Fig. 8. moderately rimed particles were observed at ?9.7°C and ?6°C. Since broadbranched crystals were predominant, aggregates of the broad-branched crystals and plates or capped columns(Figs. 8a and 8b) were observed. The rimed particle observed at ?5.9°C (0814:08 UTC, Fig. 8c) was an assemblage of a few plates.

    Large radiating plates without riming were found in both the embedded convective and stratiform regions. Since plates predominated at temperatures from 0°C to ?4°C and?8°C to ?22°C, combined with the vertical motions, it was reasonable to see plate-like crystals and plate aggregates at temperatures between ?4°C and ?6°C. The complex platelike assemblages (e.g., Fig. 8d) that originated at temperatures colder than ?8°C experienced riming and aggregation during descent.

    Fig. 8. Example CPI images of ice particles in the embedded convective region. .

    Another feature between ?4°C and ?5°C was the evident increase in needles and needle combinations. The needles coexisted with partially rimed particles during this stage. Riming on needles was not easy, as the onset of riming on columnar particles was dependent on the length along the minor axis (Ono, 1969). At widths between 50 and 90 μm, the onset of riming was critical, and riming was unlikely to occur on columnar crystals with widths of less than 50 μm. A few lightly rimed columns were also observed from the CPI. The column in Fig. 7i had a width of 87 μm, indicating a transition from pristine to rimed columns. Moderately rimed columns (Figs. 8j and 8k) were also observed.

    Heavily rimed particles with roundness values of 0.71 and 0.56 are shown in Figs. 8e and 8g. At temperatures between ?3°C and ?6°C, the ice particles were generally moderately to heavily rimed. This is understandable since the LWC increased to over 0.1 g mat temperatures between 0°C and ?5°C.

    Therefore, the rimed particles from CPI images suggested the existence of supercooled liquid water in the convective region from ?6°C to ?10°C. For the strict exclusion of mixed-phase regions, the ice region on leg 2 was determined to be associated temperatures colder than ?15°C,where the CDP concentrations were less than 0.4 cmand the King LWC was less than 0.1 g m.

    On leg 3 (Fig. 7g), the cloud droplet concentrations peaked at 322 cmat ?2.5°C but decreased rapidly to less than 1 cmat temperatures colder than ?6°C. The supercooled water was primarily distributed at temperatures warmer than ?6°C. A comparison of LWCs from the CDP and King probe suggested consistent results at levels lower than ?6°C, with peak values of 0.67 and 0.56 cm, respectively. For the LWC at temperatures colder than ?6°C, the King probe once again showed systematically higher values than the CDP due to its response to ice crystals. Therefore, temperatures between ?6°C and ?14.3°C were again classified as unidentified regions. According to the criteria of the CDP concentrations less than 1 cmand King LWCs less than 0.1 g m, the ice region occurred at temperatures colder than ?14.3°C.

    Figure 9 shows the vertical distribution of the CDP-measured particle concentrations, the CDP LWCs, and the King LWCs on legs 5 and 6. During flight legs 5 and 6, the concentrations of the cloud particles decreased dramatically, with low values of 94 cmat ?2.9°C and 40 cmat ?1.5°C,respectively. Since most cloud particles were small at temperatures warmer than ?4°C and the concentrations of particles between 24 and 50 μm were predominantly less than 0.5 cm, the measured LWCs were less than 0.1 g m. Accordingly, the top of the supercooled liquid water for legs 5 and 6 was identified at lower levels of ?5°C and ?3.6°C.

    For the stratiform region, radiating plates and rimed plate aggregations were commonly observed at temperatures colder than ?4°C, where ice particles grew predominantly via vapor deposition. By comparison, only a few moderately rimed particles were found at temperatures warmer than ?3.5°C, and riming did not contribute significantly to particle growth. Aggregates of plates could also be seen in lower regions of the cloud.

    Therefore, for the cloud region observed on leg 1, the top of the supercooled liquid water extended to ?11°C. The occurrence of dendritic precipitation particles in the ice region at temperatures between ?6.6°C and ?8.5°C on leg 1 suggested that the growth of ice particles might have depleted the supercooled liquid water before the measurement. The supercooled liquid water in the embedded convective cells primarily occurred at temperatures warmer than?10°C with maxima of approximately 0.57 g m. For flight leg 1, the precipitation started during the period between 0700 and 0800 UTC, which coincided with the same time period as the aircraft measurements. In contrast, for legs 2 and 3, the aircraft measurements started 1-h later than the precipitation onset time. The particle growth associated with the precipitation falling at an earlier time might have removed cloud droplets at temperatures between ?10°C and?15°C in the convective region.

    5.Particle size distributions

    5.1.Vertical variation in particle size distributions

    Figure 10 shows the combined particle size distributions at approximately ?3°C, ?5°C, ?10°C, and ?15°C across the five spirals. The particle size distributions were from the 5-s averaged values. The particle size distributions for leg 1 (Fig. 10a) became generally broader from ?15°C to ?3.2°C, with the maximum dimensions increasing from almost 2000 μm to 6000 μm. Large particles with maximum dimensions greater than 1000 μm increased with decreasing height, suggesting ice crystal growth from small to large sizes within the cloud. The concentrations of 300 μm particles first peaked in the upper cloud at ?15°C and then in the lower cloud at ?3.2°C, while the concentrations were relatively lower in the middle cloud layer between?5°C and ?10°C. The first peak could be explained by the initial ice nucleation which produced more small ice crystals in the upper part than in the lower part. The second peak at?3.2°C corresponded with a significant increase in particle(25 to 505 μm) concentrations from 1 Lat ?5°C to 70 Lat ?3.2°C. An overview of 2D-S images suggested predominantly columns and needles at approximately ?3°C. The high particle concentrations were assumed to be associated with secondary ice production processes. However, the cloud conditions did not fit into the Hallett-Mossop criteria due to the lack of large supercooled droplets. Knight (2012) suggested that another ice multiplication mechanism that does not involve riming occurs at temperatures of approximately?5°C. The thin needles could shatter on impact with other crystals and provide splinters similar to the Hallett-Mossop process. We will further evaluate the production mechanism of the needles within the mixed-phase stratiform clouds using more in situ observations.

    Fig. 9. Vertical distribution of the (a) CDP concentrations, (b) CDP LWC and (c) King LWC on leg 5; (d) CDP concentrations, (e) CDP LWC, and (f) King LWC on leg 6. The red horizontal lines represent the discrimination between the mixed-phase and ice clouds. .

    For the convective region (leg 2, Fig. 10b), the increase in the concentration of particles larger than 500 μm could be seen at ?10°C and ?15°C, and the particle size distributions at all the levels became even broader than those on leg 1.The significant increases in larger particles were associated with dendrites and broad-branched crystals. In comparison,for the stratiform regions (Figs. 10d and 10e), the sizes of the largest particles and concentrations of particles in the size range 200 to 1000 μm decreased dramatically. Therefore, the convective region had broader particle size distributions than the stratiform region.

    Table 4 lists the measurements for the selected temperature regions from legs 2 and 5. A comparison of the precipitation particle mean volume diameters from the HVPS is shown in the table. The particle mean volume diameter is primarily used in microphysics parametizations and related calculations. It can be expressed as

    where ρ is the density of air, Qis the mixing ratio of hydrometer species, ρis the density of the hydrometer species,and Nis the number concentration of the hydrometer species.

    It suggested that precipitation particle sizes in the convective region (leg 2) were generally larger for each temperature than those in the stratiform region (leg 5). This difference corresponded with larger maximum particle dimensions and wider particle size distributions in the convective region. The most significant increases in the particle maximum sizes on leg 2 were from ?15°C to ?10°C, with values increasing from 5925 to 10 950 μm. These differences were associated with particle habit variations from broadbranched crystals at ?15°C to large dendrites and their aggregates at ?10°C. In comparison, the maximum precipitation particle dimensions from ?15°C to ?10°C on leg 6 also increased, but lower values were observed from 2400 to 5025 μm since the proportion of smaller plate-like crystals increased in the stratiform region. Therefore, a significant shift in the particle size distributions toward larger particles occurred in the convective region from ?15°C to ?10°C.

    Fig. 10. Combined particle size distributions at various temperatures on (a) leg 1, (b) leg 2, (c) leg 3, (d) leg 5, and(e) leg 6. .

    Another difference between the convective and stratiform regions was the particle concentrations with a size range of 25-505 μm at temperatures between ?3°C and?5°C. High particle concentrations of 357 and 108 Lat?3°C and ?5°C, respectively, were observed on leg 2, while the particle concentrations on leg 5 were only 1.6 and 4.3 Lat ?3°C and ?5°C, respectively. This finding indicated secondary ice particle processes in the convective cells.Many previous studies (Rangno and Hobbs, 2001; Heymsfield and Willis, 2014; Field et al., 2017) have discussed the occurrence of columns and needles since they were thought to be associated with secondary ice production. In situ observations of stratiform clouds suggested that needles occurred most frequently at temperatures between 0°C and ?5°C(Korolev et al., 2000). Woods et al. (2008) found that needle-like crystals increased suddenly at ?7°C and that the rime-splintering process was likely responsible for new particle formation at temperatures warmer than ?10°C. We observed a large number of columns and needles across leg 2 at temperatures between ?3°C and ?6°C. The combined presence of large droplets exceeding 24 μm in diameter and the evident increase in needles within embedded convective cells suggested that the Hallett-Mossop process was responsible for the presence of columns and needles at approximately ?5°C. In contrast, needles were rarely observed on leg5, where the cloud conditions did not fit into the Hallett-Mossop criteria due to the lack of large supercooled droplets.

    Table 4. List of measurements for selected temperature regions from legs 2 and 5.

    A comparison of column and needle occurrence across the five spirals suggested that evident columns and needles occurred on legs 1, 2, and 6, but not on legs 3 and 5, which indicated inhomogeneous cloud conditions. It would be interesting to study the transportation of cloud particles with the help of horizontal penetration at temperatures warmer than?10°C. Heymsfield and Willis (2014) found that needles were observed primarily in regions of low liquid water content (< 0.1 g m) and weak vertical motions (?1 to 1 m s).This result may not apply to the embedded convective regions since the supercooled liquid water content on leg 2 was higher, peaking at 0.57 g m.

    5.2.Variation in particle size distributions at -16°C

    From 0848 to 0908 UTC, the aircraft performed five in situ circles at 7.26 km. Fig. 10 shows the observed radar reflectivity at 0848 UTC and the flight track. Radar observations at 0848 UTC (Fig. 11a) suggested that the clouds had a general reflectivity of 20-30 dBZ. Embedded convective cells with maximum values of 37 dBZ were observed in the southern part. The sample region was at the edge of the embedded convective region. Each in situ cycle lasted for approximately four minutes (Fig. 11b). Measured reflectivity values in the sample region varied from 26 dBZ at 0848 UTC, 30 dBZ at 0854 UTC, 36 dBZ at 0900 UTC to 33 dBZ at 0906 UTC. The 1-h accumulated precipitation below the sample region increased from 6.8 mm from 0800 to 0900 UTC to a peak value of 8.8 mm from 0900 to 1000 UTC.Therefore, aircraft measurements during this period corresponded to steady and moderate surface precipitation.

    The CDP-measured droplet concentrations at 7.26 km were mostly less than 0.4 cm, which was almost below the detectable limits. The King LWC at the same height was less than 0.1 g m. Hence, the sampled region at 7.26 km from both the spirals and the horizontal penetration belonged to the ice cloud region.

    Figure 12 shows the mean combined particle size distributions at 7.26 km from the five circles. A general decrease in concentrations for ice particles in the size range of 25 to 1105 μm could be seen. An overview of 2D-S concentrations suggested that concentrations of small particles in the size range of 25 to 505 μm decreased from 2.91 Lon circle 1 to 0.45 Lon circle 5, and concentrations of larger particles in the size range of 505 to 1105 μm also decreased,with values from 9.42 to 1.92 L. At the same time, the HVPS mean volume diameters generally increased from 941 μm on circle 1 to 1022 μm on circle 5. The growth of small ice crystals likely contributed to the increase in the sizes of larger particles. The 2D-S images suggested a mixture of particle habits varying from plates and broadbranched crystals to dendrites and irregulars. The well-distinguished branches suggested that most particles were pristine crystals growing predominantly via vapor deposition. Since plates and broad-branched crystals were the common crystal types at ?16°C, aggregates of such particles could be seen.

    Fig. 11. (a) PPI radar display of the radar reflectivity factor (in dBZ) at 0848 UTC and (b) the flight track during 0848 and 0908 UTC. The flight track is also marked by red circles in (a). .

    Fig. 12. Mean combined particle size distributions at 7.26 km from the horizontal penetration. .

    Figure 13 shows the time series of IWC during 0848 and 0908 UTC. According to Fig. 13, the estimated IWC fluctuated from the peak value of 0.8 g mto relatively lower values of less than 0.2 g m. Relatively higher IWCs, including 0.8 g mat 0856 UTC, 0.54 g mat 0900 UTC, and 0.55 g mat 0904 UTC, occurred in regions where there were many large radiating crystals (Fig. 14). The HVPS mean volume diameters in the three regions were over 1400 μm, which were much higher than the mean values of approximately 1100 μm. The large particle sizes indicated that the ice particles grew vigorously at that level. With the falling out of these large crystals, it was reasonable to see the general decrease in IWC from 0848 to 0908 UTC.

    Fig. 13. Time series of IWC during 0848 and 0908 UTC.

    Fig. 14. Example 2D-S images observed on flight leg 4 at 7.26 km.

    6.Summary and conclusions

    In this study, we presented the microphysical features of one mixed-phase stratiform cloud system from the Taihang Mountains Precipitation Enhancement and Hail Suppression Experiment. Based on the examination of the characteristics of supercooled liquid water and particle size distributions in stratiform precipitating regions and embedded convective cells, we have discussed riming and associated rime splintering as well as their impact on ice particle habits and concentrations.

    Supercooled liquid water was found predominantly at temperatures warmer than ?6°C and decreased substantially at colder temperatures. The top of the supercooled liquid water cloud region was observed to be at temperatures colder than ?11°C and warmer than ?15°C. Significant amounts of supercooled liquid water with maximum values of up to 0.57 g mwere observed from ?3°C to ?6°C in the embedded convective cells, while the peak LWC was generally less than 0.1 g min the stratiform regions.

    For the convective region, moderately to heavily rimed ice particles were observed at temperatures warmer than -10°C. The rimed particles grew predominantly from plates and broad-branched crystals, with few from columns. For the stratiform region, ice particles grow predominantly through vapor deposition, and only a few moderately rimed particles are found at temperatures warmer than ?3.5°C.Thus, riming does not contribute significantly to particle growth in the stratiform regions.

    There are two distinct differences in particle size distributions between convective and stratiform regions. One difference is that the convective region has generally broader particle size distributions than the stratiform region, especially at temperatures colder than ?10°C. For the case in this study, a significant shift in particle size distributions toward larger particles occurs in the convective region from?15°C to ?10°C, with the maximum particle dimensions increasing from less than 6000 μm to over 1 cm. This change is due to habit variations from broad-branched crystals to large dendrites and their aggregates.

    The other difference is that large concentrations of small particles (25-205 μm) are observed at temperatures between ?3°C and ?5°C in the convective region, where rimed ice particles and needles coexist. The cloud conditions within the convective region fit into the Hallett-Mossop criteria, and the associated rime splintering is assumed to be responsible for the observed high concentrations at those temperature ranges. In comparison, the cloud conditions within the stratiform regions do not fit into the Hallett-Mossop criteria due to the lack of large supercooled droplets, and the small particle concentrations are considerably lower.

    The updraft data are useful to demonstrate the intensity of the embedded convection and the substantial changes in crystal habits due to vertical motions. As the accuracy of vertical wind measurements has not been well assessed, the updraft data were not analyzed in this study. Further aircraft measurements are required to analyze the spatial variability of cloud microstructure in stratiform clouds with embedded convection. Modeling work is required to test the impact of the riming process on particle size distributions and particle fall speeds.

    Acknowledgements. This work was supported by the National Key Research and Development Program of China (Grant Nos. 2019YFC1510300 and 2018YFC1507900) and the National Natural Science Foundation of China (Grant Nos. 41575131).

    97精品久久久久久久久久精品| 又大又爽又粗| 亚洲精品一区蜜桃| 午夜福利免费观看在线| 成年人黄色毛片网站| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜添小说| 韩国高清视频一区二区三区| 久久久精品免费免费高清| 国产亚洲欧美在线一区二区| 亚洲精品日韩在线中文字幕| 可以免费在线观看a视频的电影网站| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 亚洲熟女毛片儿| 免费在线观看影片大全网站| 亚洲欧美色中文字幕在线| 日韩免费高清中文字幕av| 在线亚洲精品国产二区图片欧美| 国产xxxxx性猛交| 亚洲视频免费观看视频| 日韩视频一区二区在线观看| 中亚洲国语对白在线视频| 后天国语完整版免费观看| 国产麻豆69| 亚洲avbb在线观看| 丰满人妻熟妇乱又伦精品不卡| 黄色怎么调成土黄色| 91九色精品人成在线观看| 免费日韩欧美在线观看| 老司机靠b影院| 国产精品影院久久| kizo精华| 中文字幕人妻丝袜制服| 国产野战对白在线观看| 日本精品一区二区三区蜜桃| 久久人妻熟女aⅴ| 中国美女看黄片| 制服诱惑二区| 国产成人啪精品午夜网站| 免费人妻精品一区二区三区视频| 中文字幕人妻熟女乱码| 丰满人妻熟妇乱又伦精品不卡| 久久ye,这里只有精品| 国产日韩一区二区三区精品不卡| 美国免费a级毛片| 手机成人av网站| 亚洲国产欧美日韩在线播放| 欧美黑人精品巨大| 成人亚洲精品一区在线观看| 久久精品成人免费网站| 少妇人妻久久综合中文| 淫妇啪啪啪对白视频 | 国产亚洲欧美在线一区二区| 亚洲精品久久久久久婷婷小说| 操出白浆在线播放| av国产精品久久久久影院| 热99国产精品久久久久久7| 人人妻,人人澡人人爽秒播| avwww免费| 满18在线观看网站| 亚洲天堂av无毛| 日韩三级视频一区二区三区| 亚洲av电影在线进入| av超薄肉色丝袜交足视频| 日本av免费视频播放| 久久综合国产亚洲精品| 色婷婷久久久亚洲欧美| av线在线观看网站| 久久久精品免费免费高清| 一级片免费观看大全| 人人妻人人澡人人看| 国产成人啪精品午夜网站| 9热在线视频观看99| www.自偷自拍.com| 久久久久久久国产电影| 亚洲人成77777在线视频| 亚洲av欧美aⅴ国产| 亚洲精品粉嫩美女一区| 久久天躁狠狠躁夜夜2o2o| 精品人妻在线不人妻| 国产精品久久久久久精品古装| 亚洲精品一卡2卡三卡4卡5卡 | 女性被躁到高潮视频| 最新的欧美精品一区二区| 欧美日本中文国产一区发布| 久久精品国产a三级三级三级| 黄网站色视频无遮挡免费观看| 美国免费a级毛片| www.熟女人妻精品国产| 国产精品秋霞免费鲁丝片| 80岁老熟妇乱子伦牲交| 国产av又大| 国产无遮挡羞羞视频在线观看| 9热在线视频观看99| 人妻人人澡人人爽人人| 女人爽到高潮嗷嗷叫在线视频| 精品国产乱子伦一区二区三区 | 精品免费久久久久久久清纯 | 亚洲精品粉嫩美女一区| 亚洲av成人一区二区三| xxxhd国产人妻xxx| 欧美国产精品一级二级三级| 成人国产av品久久久| 韩国精品一区二区三区| 2018国产大陆天天弄谢| 高清黄色对白视频在线免费看| 青春草视频在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 一边摸一边做爽爽视频免费| 在线十欧美十亚洲十日本专区| av电影中文网址| 色播在线永久视频| 麻豆乱淫一区二区| 午夜影院在线不卡| 色94色欧美一区二区| 精品国产乱码久久久久久男人| 69av精品久久久久久 | 免费观看av网站的网址| 日韩一卡2卡3卡4卡2021年| 午夜免费鲁丝| 黄色视频在线播放观看不卡| 自线自在国产av| 在线看a的网站| 视频在线观看一区二区三区| 欧美少妇被猛烈插入视频| 国产在视频线精品| 狠狠精品人妻久久久久久综合| 午夜福利一区二区在线看| 久久人妻福利社区极品人妻图片| bbb黄色大片| 国产片内射在线| 亚洲av片天天在线观看| 日韩制服丝袜自拍偷拍| 美女高潮到喷水免费观看| 免费观看a级毛片全部| 搡老乐熟女国产| 女性被躁到高潮视频| 亚洲少妇的诱惑av| 亚洲国产精品999| 麻豆国产av国片精品| 精品一区在线观看国产| a 毛片基地| 在线观看人妻少妇| 51午夜福利影视在线观看| 女人精品久久久久毛片| 中文资源天堂在线| 一区二区三区高清视频在线| 色综合亚洲欧美另类图片| 青草久久国产| 国产人伦9x9x在线观看| 国产精品亚洲美女久久久| 亚洲欧美日韩高清专用| 日本在线视频免费播放| 黄色a级毛片大全视频| 亚洲中文字幕日韩| 亚洲欧美日韩无卡精品| 成人国产一区最新在线观看| 成人三级黄色视频| 国产成人精品无人区| 国产精品亚洲av一区麻豆| 欧美中文日本在线观看视频| 精品国产乱子伦一区二区三区| 又紧又爽又黄一区二区| 免费高清视频大片| 我要搜黄色片| 久久久久性生活片| 日本黄色视频三级网站网址| 午夜a级毛片| 俄罗斯特黄特色一大片| 蜜桃久久精品国产亚洲av| 看免费av毛片| 色综合婷婷激情| 99国产精品一区二区蜜桃av| 国产亚洲精品久久久久5区| 国内揄拍国产精品人妻在线| 成年版毛片免费区| 观看免费一级毛片| 99久久精品国产亚洲精品| 一本一本综合久久| 美女 人体艺术 gogo| 亚洲欧美日韩高清专用| 男插女下体视频免费在线播放| 一本一本综合久久| 欧美色视频一区免费| 亚洲精华国产精华精| 视频区欧美日本亚洲| 久久中文字幕人妻熟女| 国产精品精品国产色婷婷| 最近最新中文字幕大全免费视频| 91av网站免费观看| 亚洲精品美女久久久久99蜜臀| 午夜福利免费观看在线| 色播亚洲综合网| 十八禁人妻一区二区| 精品熟女少妇八av免费久了| 99热这里只有是精品50| 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 无遮挡黄片免费观看| 久久中文字幕人妻熟女| 精品欧美一区二区三区在线| 亚洲成av人片在线播放无| 成人欧美大片| 男人舔奶头视频| 九九热线精品视视频播放| 免费在线观看成人毛片| 欧美久久黑人一区二区| 日本黄色视频三级网站网址| 免费高清视频大片| 三级毛片av免费| 日韩欧美在线乱码| 在线视频色国产色| 久久精品国产亚洲av高清一级| 成熟少妇高潮喷水视频| 国产精品国产高清国产av| 色噜噜av男人的天堂激情| 欧美在线一区亚洲| 日本在线视频免费播放| 久久香蕉精品热| 免费在线观看完整版高清| a级毛片在线看网站| 久久伊人香网站| 亚洲专区国产一区二区| 国内揄拍国产精品人妻在线| 色精品久久人妻99蜜桃| 国产在线精品亚洲第一网站| 老汉色∧v一级毛片| 久久香蕉国产精品| 97超级碰碰碰精品色视频在线观看| 一本综合久久免费| 精品久久久久久,| 精品国产乱码久久久久久男人| e午夜精品久久久久久久| 亚洲成av人片在线播放无| 日本成人三级电影网站| 日本精品一区二区三区蜜桃| 国产高清有码在线观看视频 | 老司机午夜福利在线观看视频| 国产又黄又爽又无遮挡在线| 亚洲天堂国产精品一区在线| 哪里可以看免费的av片| 亚洲成人国产一区在线观看| 国产真人三级小视频在线观看| 91在线观看av| 亚洲欧美日韩东京热| 午夜精品一区二区三区免费看| 99久久99久久久精品蜜桃| 成年女人毛片免费观看观看9| 久久欧美精品欧美久久欧美| 日本成人三级电影网站| 小说图片视频综合网站| 伊人久久大香线蕉亚洲五| 国产三级中文精品| 女同久久另类99精品国产91| 99国产极品粉嫩在线观看| 好男人电影高清在线观看| 国产免费av片在线观看野外av| 丰满的人妻完整版| 成人国产一区最新在线观看| 久久这里只有精品中国| or卡值多少钱| 欧美日韩乱码在线| 国产不卡一卡二| 黄色a级毛片大全视频| 亚洲中文日韩欧美视频| 国产精品乱码一区二三区的特点| 欧美三级亚洲精品| 巨乳人妻的诱惑在线观看| 日韩三级视频一区二区三区| 午夜亚洲福利在线播放| 亚洲成人久久爱视频| 中文字幕最新亚洲高清| 9191精品国产免费久久| 小说图片视频综合网站| 男人舔奶头视频| 深夜精品福利| 中文字幕熟女人妻在线| 久久精品夜夜夜夜夜久久蜜豆 | 欧美另类亚洲清纯唯美| 777久久人妻少妇嫩草av网站| 欧美日本亚洲视频在线播放| 免费搜索国产男女视频| 男女那种视频在线观看| 亚洲 国产 在线| 男女之事视频高清在线观看| 色综合欧美亚洲国产小说| 国产av麻豆久久久久久久| 他把我摸到了高潮在线观看| 亚洲精品久久国产高清桃花| 国产片内射在线| 国产精品久久久av美女十八| 桃色一区二区三区在线观看| 午夜两性在线视频| 长腿黑丝高跟| 美女 人体艺术 gogo| 国产久久久一区二区三区| 岛国在线观看网站| 久久人人精品亚洲av| 99在线人妻在线中文字幕| 国产69精品久久久久777片 | 亚洲人成伊人成综合网2020| 99国产精品一区二区三区| 国产又色又爽无遮挡免费看| 日韩大码丰满熟妇| 国产成人精品久久二区二区免费| 亚洲精品久久成人aⅴ小说| 丝袜人妻中文字幕| 一级毛片高清免费大全| av欧美777| 这个男人来自地球电影免费观看| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线| 成人亚洲精品av一区二区| 国产一区二区三区视频了| 叶爱在线成人免费视频播放| 久久热在线av| 国产在线观看jvid| 在线观看美女被高潮喷水网站 | 国产高清videossex| www国产在线视频色| 亚洲国产精品合色在线| 亚洲av美国av| 欧美黄色淫秽网站| 欧美日韩瑟瑟在线播放| 波多野结衣高清作品| 麻豆成人午夜福利视频| 在线观看免费视频日本深夜| www.自偷自拍.com| 精品免费久久久久久久清纯| 99国产精品99久久久久| 国产又色又爽无遮挡免费看| 麻豆国产av国片精品| 老司机深夜福利视频在线观看| 99久久精品热视频| 国产97色在线日韩免费| 国产精品乱码一区二三区的特点| 精品国产美女av久久久久小说| 无人区码免费观看不卡| 成年版毛片免费区| 国产精品日韩av在线免费观看| 国产真实乱freesex| 久久久久国产精品人妻aⅴ院| 黄片小视频在线播放| 麻豆国产av国片精品| 久久精品91蜜桃| 这个男人来自地球电影免费观看| 一本综合久久免费| 99热只有精品国产| 日韩欧美国产在线观看| 无遮挡黄片免费观看| 免费人成视频x8x8入口观看| 国产精品乱码一区二三区的特点| 亚洲熟妇熟女久久| 国产精品久久电影中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 一个人免费在线观看的高清视频| 高清在线国产一区| 久久久精品国产亚洲av高清涩受| 男人舔奶头视频| 18禁黄网站禁片午夜丰满| 黄色 视频免费看| 啪啪无遮挡十八禁网站| 免费一级毛片在线播放高清视频| 久久久水蜜桃国产精品网| 性欧美人与动物交配| 久久香蕉精品热| 深夜精品福利| 精品日产1卡2卡| 两个人的视频大全免费| 51午夜福利影视在线观看| 蜜桃久久精品国产亚洲av| 亚洲国产看品久久| av超薄肉色丝袜交足视频| 国产97色在线日韩免费| 不卡一级毛片| 一二三四社区在线视频社区8| 1024手机看黄色片| av福利片在线观看| 国产精品自产拍在线观看55亚洲| 精品一区二区三区av网在线观看| 精品久久久久久,| 丝袜人妻中文字幕| 国产激情偷乱视频一区二区| 亚洲精品在线观看二区| 午夜免费激情av| 欧美性长视频在线观看| 免费高清视频大片| 亚洲第一欧美日韩一区二区三区| 国内少妇人妻偷人精品xxx网站 | 国产午夜精品久久久久久| 最新在线观看一区二区三区| 身体一侧抽搐| 国产高清有码在线观看视频 | 国产高清有码在线观看视频 | 国产1区2区3区精品| 最好的美女福利视频网| www.www免费av| 男女做爰动态图高潮gif福利片| 日本 欧美在线| 中出人妻视频一区二区| 日韩中文字幕欧美一区二区| 久久精品aⅴ一区二区三区四区| 国产熟女午夜一区二区三区| 嫩草影视91久久| 狂野欧美白嫩少妇大欣赏| 又粗又爽又猛毛片免费看| 午夜激情福利司机影院| 午夜免费成人在线视频| 最新在线观看一区二区三区| 日韩 欧美 亚洲 中文字幕| www.精华液| 在线永久观看黄色视频| 精品一区二区三区av网在线观看| 国产一区二区激情短视频| 午夜日韩欧美国产| 天天添夜夜摸| 老司机深夜福利视频在线观看| 久久精品国产亚洲av高清一级| 国产精品一区二区三区四区免费观看 | 88av欧美| 亚洲精品国产一区二区精华液| 亚洲全国av大片| 一夜夜www| 午夜日韩欧美国产| 国产精品一区二区三区四区久久| а√天堂www在线а√下载| 中文字幕久久专区| 午夜免费观看网址| 中文字幕高清在线视频| 国产精品亚洲一级av第二区| 狠狠狠狠99中文字幕| 免费一级毛片在线播放高清视频| 亚洲专区字幕在线| 久久香蕉国产精品| 高潮久久久久久久久久久不卡| 色噜噜av男人的天堂激情| 免费在线观看黄色视频的| 免费在线观看视频国产中文字幕亚洲| 亚洲中文av在线| 精品欧美国产一区二区三| 91九色精品人成在线观看| 午夜免费成人在线视频| 成人亚洲精品av一区二区| 久久精品国产亚洲av香蕉五月| 十八禁人妻一区二区| 亚洲av五月六月丁香网| 久久人妻福利社区极品人妻图片| 国产精品亚洲av一区麻豆| 悠悠久久av| 精品一区二区三区四区五区乱码| 老熟妇乱子伦视频在线观看| 午夜精品久久久久久毛片777| 亚洲欧洲精品一区二区精品久久久| 每晚都被弄得嗷嗷叫到高潮| 亚洲七黄色美女视频| 熟女少妇亚洲综合色aaa.| 不卡av一区二区三区| 免费在线观看完整版高清| 国产区一区二久久| 久久久久亚洲av毛片大全| 啪啪无遮挡十八禁网站| 亚洲国产精品合色在线| 天堂√8在线中文| 国产视频内射| 久久精品亚洲精品国产色婷小说| 男女下面进入的视频免费午夜| 亚洲成a人片在线一区二区| 国产精品久久电影中文字幕| 精品久久久久久,| 国模一区二区三区四区视频 | 丰满人妻一区二区三区视频av | 欧美日韩一级在线毛片| 中文字幕熟女人妻在线| 色老头精品视频在线观看| 国产成人影院久久av| 免费在线观看影片大全网站| 亚洲免费av在线视频| 国产高清激情床上av| 久久久久国产一级毛片高清牌| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人不卡在线观看播放网| 日韩中文字幕欧美一区二区| 亚洲精品久久国产高清桃花| 久久香蕉激情| 国产伦在线观看视频一区| 午夜福利在线在线| 制服诱惑二区| 999久久久精品免费观看国产| 午夜福利在线观看吧| 亚洲色图av天堂| av在线天堂中文字幕| 日韩欧美在线乱码| 又黄又粗又硬又大视频| 九九热线精品视视频播放| 欧美最黄视频在线播放免费| 老汉色av国产亚洲站长工具| 免费在线观看成人毛片| 午夜视频精品福利| 色综合站精品国产| 久久久久免费精品人妻一区二区| or卡值多少钱| 亚洲av中文字字幕乱码综合| 久久中文字幕人妻熟女| 亚洲一区二区三区不卡视频| 亚洲av日韩精品久久久久久密| 成人亚洲精品av一区二区| avwww免费| 亚洲成a人片在线一区二区| 黄色 视频免费看| 午夜福利在线在线| 色综合亚洲欧美另类图片| 国产成年人精品一区二区| 亚洲精品粉嫩美女一区| 午夜福利视频1000在线观看| 亚洲在线自拍视频| 美女大奶头视频| 日韩大码丰满熟妇| 中文字幕av在线有码专区| 国产av在哪里看| 99热只有精品国产| 亚洲人成电影免费在线| 日韩大尺度精品在线看网址| 国语自产精品视频在线第100页| 九九热线精品视视频播放| 淫秽高清视频在线观看| 亚洲人与动物交配视频| 曰老女人黄片| 亚洲国产欧美人成| 欧美国产日韩亚洲一区| 女同久久另类99精品国产91| 黄频高清免费视频| 国产区一区二久久| 曰老女人黄片| 日韩 欧美 亚洲 中文字幕| 一级作爱视频免费观看| 国产1区2区3区精品| 999久久久精品免费观看国产| 俄罗斯特黄特色一大片| 特大巨黑吊av在线直播| 午夜两性在线视频| 国产午夜福利久久久久久| 久久这里只有精品中国| 欧美人与性动交α欧美精品济南到| 看片在线看免费视频| 亚洲成人精品中文字幕电影| 这个男人来自地球电影免费观看| 男女那种视频在线观看| 可以免费在线观看a视频的电影网站| 一区二区三区激情视频| 深夜精品福利| 国产亚洲av高清不卡| 波多野结衣高清作品| 色尼玛亚洲综合影院| 午夜久久久久精精品| 男人舔奶头视频| 母亲3免费完整高清在线观看| 日日摸夜夜添夜夜添小说| 日韩精品免费视频一区二区三区| 一进一出抽搐动态| 在线观看免费视频日本深夜| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| 亚洲在线自拍视频| 国产不卡一卡二| 免费在线观看成人毛片| 欧美日韩黄片免| 久久亚洲真实| 无限看片的www在线观看| 亚洲熟妇中文字幕五十中出| 动漫黄色视频在线观看| 最好的美女福利视频网| 制服诱惑二区| 少妇熟女aⅴ在线视频| 床上黄色一级片| 国产精品永久免费网站| 村上凉子中文字幕在线| 一级作爱视频免费观看| 国产精品99久久99久久久不卡| 成在线人永久免费视频| 亚洲人成网站高清观看| 999久久久精品免费观看国产| 亚洲黑人精品在线| 亚洲国产欧美网| 国产欧美日韩精品亚洲av| 99久久综合精品五月天人人| 一进一出抽搐gif免费好疼| 日本 欧美在线| 脱女人内裤的视频| 国产免费av片在线观看野外av| 国产单亲对白刺激| 欧美+亚洲+日韩+国产| 午夜精品在线福利| 国产亚洲精品久久久久久毛片| 午夜亚洲福利在线播放| 国产亚洲精品综合一区在线观看 | 午夜福利免费观看在线| 最新美女视频免费是黄的| 国产高清激情床上av| 免费在线观看亚洲国产| ponron亚洲| 给我免费播放毛片高清在线观看| 亚洲一区中文字幕在线| 长腿黑丝高跟| 精品一区二区三区av网在线观看| 巨乳人妻的诱惑在线观看| 首页视频小说图片口味搜索| 男人舔奶头视频| 日本一区二区免费在线视频| 999精品在线视频| 毛片女人毛片| 国产精品久久久久久久电影 | 老司机福利观看| 精品久久久久久,| 黄色 视频免费看|