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

    華北地區(qū)一次低槽冷鋒云系結(jié)構(gòu)和增雨條件研究

    2023-01-14 05:21:22史月琴趙俊杰孫玉穩(wěn)
    干旱氣象 2022年6期
    關(guān)鍵詞:云系鋒面冷鋒

    史月琴,趙俊杰,孫 晶,孫玉穩(wěn)

    (1.中國(guó)氣象局云降水物理與人工影響天氣重點(diǎn)開(kāi)放實(shí)驗(yàn)室,北京 100081;2.中國(guó)氣象局人工影響天氣中心,北京 100081;3.山西省氣象災(zāi)害防御技術(shù)中心,山西 太原 030032;4.河北省人工影響天氣辦公室,河北 石家莊 050021)

    引 言

    華北地區(qū)處于暖溫帶半潤(rùn)濕半干旱區(qū),年平均降水量400~800 mm,降水季節(jié)分配不均,年際變化大,旱災(zāi)發(fā)生頻繁。地處華北地區(qū)的河北、山西2017年春季發(fā)生嚴(yán)重干旱,對(duì)農(nóng)作物播種及生長(zhǎng)造成一定影響[1],抓住有利時(shí)機(jī)開(kāi)展人工增雨作業(yè)有利于緩解旱情、減輕干旱災(zāi)害影響。低槽冷鋒是河北省春秋季的主要降水天氣系統(tǒng),形成的大范圍層狀云系是當(dāng)?shù)厝斯ぴ鲇甑闹饕鳂I(yè)對(duì)象[2],但并不是云系的所有區(qū)域都適合開(kāi)展人工增雨作業(yè),如何確定合適的增雨時(shí)機(jī)和部位是科學(xué)實(shí)施增雨作業(yè)的關(guān)鍵。研究發(fā)現(xiàn)鋒前、鋒區(qū)、鋒后的云系結(jié)構(gòu)和降水機(jī)制差異明顯,位于地面鋒面上的窄冷鋒雨帶,其云系為冷暖混合云,自上而下分別是冰晶、雪、霰和過(guò)冷水、液態(tài)水,過(guò)冷水的分布局地性強(qiáng),而位于鋒線稍后的高空鋒區(qū)屬于寬冷鋒雨帶,常為高層冷云[3-8]。在云的發(fā)展、維持階段,上升氣流供應(yīng)的水汽可維持人工增雨的潛力[9],同時(shí)催化引起的水凝物相態(tài)轉(zhuǎn)化帶來(lái)潛熱釋放,進(jìn)而引起局部增溫和上升速度增加,也可促進(jìn)云和降水發(fā)展。春季出現(xiàn)的低槽冷鋒云系不同發(fā)展階段的云系結(jié)構(gòu)及增雨條件研究有助于提高人工增雨作業(yè)效率,進(jìn)而緩解當(dāng)?shù)睾登椤?/p>

    人工增雨不同于人工激發(fā)降水,考慮增雨作業(yè)的經(jīng)濟(jì)需求和實(shí)際作業(yè)效果,通常選擇雨強(qiáng)大于0.1 mm·h-1的雨區(qū)開(kāi)展[9]。自然云中冰晶濃度較低有利于增雨作業(yè),陶樹(shù)旺等[10]研究認(rèn)為機(jī)載2D-C探測(cè)的大粒子數(shù)濃度小于20 個(gè)·L-1時(shí),可確定為強(qiáng)可播區(qū),否則為可播區(qū)。

    經(jīng)典理論認(rèn)為,催化引入的人工冰晶可通過(guò)貝吉隆過(guò)程將過(guò)冷云水轉(zhuǎn)化為降水。胡志晉[9]指出人工冰晶還可使部分冰面過(guò)飽和水汽轉(zhuǎn)化為降水,這個(gè)觀點(diǎn)在碘化銀(AgI)催化數(shù)值模擬試驗(yàn)中得到了印證[11-13]。AgI 人工冰核有4 種核化機(jī)制[14],即接觸凍結(jié)、浸沒(méi)凍結(jié)核化消耗過(guò)冷水,凝華、凝結(jié)凍結(jié)核化消耗過(guò)飽和水汽。在沒(méi)有過(guò)冷水的區(qū)域,AgI 以凝華核化為主,發(fā)生條件是云中溫度低于-5 ℃且冰面飽和比大于1.04[15]。孫晶等[7]采用逐時(shí)雨強(qiáng)、上升運(yùn)動(dòng)區(qū)、過(guò)冷水含量、冰晶數(shù)濃度及溫度構(gòu)建增雨條件指標(biāo),發(fā)現(xiàn)西風(fēng)槽降水云系中槽線附近存在強(qiáng)可播區(qū),但該指標(biāo)體系沒(méi)有考慮云中冰面過(guò)飽和水汽對(duì)增雨的潛力。而在層狀云降水形成過(guò)程中,冰面過(guò)飽和水汽對(duì)凝華增長(zhǎng)有較大作用[4]。

    目前還鮮少有研究選擇一個(gè)特定區(qū)域,以時(shí)間演變的角度對(duì)經(jīng)過(guò)本地的低槽冷鋒云系展開(kāi)研究。本文以邢臺(tái)地區(qū)為例,研究不同降水階段云系結(jié)構(gòu)和增雨條件,以期為科學(xué)、精準(zhǔn)、安全實(shí)施人工增雨作業(yè)提供技術(shù)支撐。

    1 資料與方法

    所用的NCEP FNL 再分析資料(https://rda.ucar.edu/datasets/ds083.2/),水平分辨率為1°×1°,時(shí)間間隔6 h,用于分析天氣形勢(shì)并作為初始場(chǎng)和側(cè)邊界條件啟動(dòng)云降水顯式預(yù)報(bào)系統(tǒng)(Cloud and Precipitation Explicit Forecast System,CPEFS_v1.0)。FY-2F

    衛(wèi)星黑體亮溫?cái)?shù)據(jù)來(lái)源于風(fēng)云衛(wèi)星遙感數(shù)據(jù)服務(wù)網(wǎng)(http://data.nsmc.org.cn/portalsite/Data/Satellite.aspx)。中國(guó)地面自動(dòng)氣象站與CMORPH(Climate Prediction Center MORPHing technique,CMORPH)衛(wèi)星反演降水融合的逐小時(shí)降水量產(chǎn)品(V1.0)[16],水平分辨率為0.1°×0.1°,地面自動(dòng)站雨量資料均來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)。文中附圖涉及的地圖基于國(guó)家測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號(hào)為GS(2021)1328 號(hào)的標(biāo)準(zhǔn)地圖制作,底圖無(wú)修改。

    模擬采用CPEFS_v1.0 系統(tǒng),該系統(tǒng)基于WRF動(dòng)力框架耦合了中國(guó)氣象科學(xué)研究院(Chinese Academy of Meteorological Sciences,CAMS)云微物理方案[17-18]。模擬時(shí)段為5月22日08:00 至23日08:00,模式水平分辨率為3 km,中心點(diǎn)位于114.5°E、37.5°N,水平格點(diǎn)數(shù)400×350,垂直方向上分為不等間距32 層,模式頂層氣壓為50 hPa,時(shí)間步長(zhǎng)15 s,10 min 輸出一次模擬結(jié)果。采用MYJ 邊界層方案、RRTM 長(zhǎng)波輻射方案、Dudhia短波輻射方案,由于模式水平分辨率達(dá)到對(duì)流可分辨尺度,故沒(méi)有采用積云參數(shù)化方案。

    利用CPEFS_v1.0,結(jié)合實(shí)況資料,對(duì)2017年5月22日一次低槽冷鋒系統(tǒng)經(jīng)過(guò)邢臺(tái)站時(shí)的水汽、動(dòng)力、云宏微觀結(jié)構(gòu)特征進(jìn)行研究,在已往研究基礎(chǔ)上構(gòu)建增雨條件數(shù)值模式綜合識(shí)別指標(biāo),利用該指標(biāo)體系對(duì)不同降水階段與不同區(qū)域的增雨可播區(qū)進(jìn)行識(shí)別。

    2 天氣形勢(shì)與降水過(guò)程概況

    2017年5月22日08:00(北京時(shí),下同),我國(guó)西北地區(qū)東部到華北地區(qū)500 hPa 多短波擾動(dòng),由700 hPa 低渦中心延伸出的切變線經(jīng)過(guò)內(nèi)蒙古東南部—河北北部—山西北部,500 hPa槽線明顯落后于700 hPa切變線,為典型的后傾槽[圖1(a)]。850 hPa高度場(chǎng)上山西、河北處于倒槽影響下,河北北部—山西北、中部存在明顯的切變錢,山西與河北的中部、南部比濕超過(guò)9 g·kg-1;地面伴隨有冷鋒,地面鋒線位于遼寧中部—河北中部—山西中部—陜西中部一線,河北南部受偏南風(fēng)影響[圖1(b)]。大范圍低槽冷鋒云系呈帶狀分布,主要位于地面鋒線之后、500 hPa 高空槽前,08:00 邢臺(tái)上空還沒(méi)有云系覆蓋。

    低槽冷鋒在東移過(guò)程中發(fā)展加強(qiáng),22日14:00邢臺(tái)位于地面鋒線附近的云系前緣,20:00,700 hPa切變線壓在邢臺(tái)上空,邢臺(tái)處于地面鋒線后及500 hPa槽前的冷鋒云系中間部位[圖1(c)、圖1(d)]。23日02:00,500 hPa高空槽已移過(guò)邢臺(tái),邢臺(tái)地區(qū)處于冷鋒云系尾部。

    圖1 2017年5月22日08:00(a、b)與20:00(c、d)位勢(shì)高度、溫度、風(fēng)場(chǎng)、FY-2F衛(wèi)星黑體亮溫及比濕分布[(a、c)圖中,紅色虛等值線表示500 hPa溫度(單位:℃),藍(lán)色和黑色實(shí)等值線分別表示500 hPa和700 hPa位勢(shì)高度(單位:dagpm),風(fēng)矢(單位:m·s-1)表示700 hPa風(fēng)場(chǎng),填色表示FY-2F黑體亮溫(單位:K)。(b、d)圖中,紫色實(shí)等值線表示850 hPa位勢(shì)高度(單位:dagpm),填色表示比濕(單位:g·kg-1),風(fēng)矢(單位:m·s-1)表示10 m高度風(fēng)場(chǎng)]Fig.1 The distribution of geopotential height,temperature,wind field,black body temperature observed by satellite FY-2F and specific humidity at 08:00 BST (a,b) and 20:00 BST (c,d) May 22,2017(In fig.a and fig.c,the red dashed isolines are temperature (Unit:℃) at 500 hPa,the blue and black solid contours are geopotential height (Unit:dagpm) at 500 hPa and 700 hPa,respectively,wind vectors (Unit:m·s-1) are wind field at 700 hPa,the color shaded is black body temperature of FY-2F (Unit:K).In fig.b and fig.d,the purple solid contours are geopotential height (Unit:dagpm)at 850 hPa,the color shaded is specific humidity (Unit:g·kg-1) ,and wind vectors are wind field at 10 m height (Unit:m·s-1))

    受低槽冷鋒系統(tǒng)影響,5月22—23日華北地區(qū)雨帶呈東北—西南走向,與云帶走向一致,并隨系統(tǒng)自西北向東南移動(dòng)(圖2)。700 hPa 切變線與地面鋒線相距2 個(gè)緯距左右,致使靠近地面鋒線處的雨帶窄而雨強(qiáng)大[19],3 h 內(nèi)出現(xiàn)降水的區(qū)域約200~400 km寬,各地24 h降雨量普遍超過(guò)10 mm,其中山西南部、河南北部、河北西南部等地出現(xiàn)25~50 mm的大雨,河北西南部局地出現(xiàn)70~120 mm 的暴雨。此次過(guò)程是2017年入春以來(lái)最明顯的降水過(guò)程,河北及山西等地氣象部門抓住有利時(shí)機(jī),開(kāi)展了飛機(jī)和地面增雨作業(yè),助力各地緩解旱情。

    3 模擬結(jié)果驗(yàn)證

    邢臺(tái)站(114.36°E、37.18°N)位于太行山東側(cè),緊鄰太行山,海拔高度183 m。利用邢臺(tái)站10 min雨量資料、中國(guó)地面自動(dòng)氣象站與CMORPH 融合逐小時(shí)降水量產(chǎn)品(v1.0)作為實(shí)況對(duì)模擬結(jié)果進(jìn)行驗(yàn)證。

    從實(shí)況和模擬的逐3 h 雨量分布(圖2)可知,雨區(qū)呈東北—西南走向的狹長(zhǎng)帶狀分布,雨帶隨著鋒面系統(tǒng)逐漸東移南壓影響山西中部、南部及河北中部、河南、山東等地,降水大值區(qū)主要位于山西南部、河南北部及山西與河北交界處的太行山區(qū)。可以看出,CPEFS_v1.0 模擬出了靠近地面鋒線的窄冷鋒雨帶與高空鋒區(qū)造成的寬冷鋒雨帶,模擬的雨帶落區(qū)、范圍、形狀與實(shí)況基本一致,模擬的強(qiáng)降水中心累積雨量比實(shí)況略大。

    圖2 實(shí)況(左)和模擬(右)的5月22日3 h地面雨量分布(單位:mm)(黑點(diǎn)為邢臺(tái)站)Fig.2 The distribution of observed (the left) and simulated (the right) 3 h rainfall on May 22,2017 (Unit:mm)(The black dot is Xingtai station)

    邢臺(tái)站24 h 降雨量為42.6 mm,根據(jù)降水量等級(jí)標(biāo)準(zhǔn)[20]達(dá)到大雨量級(jí)。從邢臺(tái)站逐10 min 實(shí)況雨量演變(圖3)可知,降水從22日14:40 開(kāi)始,持續(xù)到23日02:40 結(jié)束。其中14:50—15:10 逐10 min雨量超過(guò)2.0 mm,15:20—17:00 逐10 min 雨量均大于或接近1.0 mm,此階段降水強(qiáng)度較大。

    圖3 2017年5月22日13:00至23日04:00邢臺(tái)站逐10 min降水量演變Fig.3 The evolution of observed and simulated 10-minute rainfall from 13:00 BST 22 to 04:00 BST 23 May,2017 at Xingtai station

    從CPEFS_v1.0 模擬結(jié)果看,模擬的逐10 min 雨量22日14:30 為0.4 mm,15:10 達(dá)到最大2.7 mm,比實(shí)況(2.9 mm)僅小0.2 mm;16:20—17:00 降水強(qiáng)度明顯減弱,17:10 后雨強(qiáng)變化平緩,23日04:00 降水結(jié)束。該系統(tǒng)較好地模擬了本次冷鋒降水的演變過(guò)程,模擬的降水開(kāi)始時(shí)間、峰值時(shí)段與實(shí)測(cè)基本一致。

    實(shí)況與模擬的逐10 min 雨量均體現(xiàn)出低槽冷鋒降水系統(tǒng)經(jīng)過(guò)邢臺(tái)本站的特征及變化,先出現(xiàn)的較強(qiáng)降水體現(xiàn)了靠近地面鋒線的窄冷鋒雨帶降水特征,之后出現(xiàn)降水強(qiáng)度小但較均勻、持續(xù)時(shí)間長(zhǎng)的寬冷鋒雨帶降水特征。

    4 不同降水時(shí)段云系特征

    按照邢臺(tái)站逐10 min 雨量變化,依據(jù)冷鋒云系特征,將22日13:30—17:00、22日17:10 至23日01:20、23日01:30—04:00劃分為鋒面抬升強(qiáng)雨期、深厚層云小雨期和降水消散期,分別對(duì)各階段的水汽、動(dòng)力條件及云系特征進(jìn)行分析。

    4.1 鋒面抬升強(qiáng)雨期

    從圖4(a)看出,邢臺(tái)站300 hPa 以下都有水汽輸送,22日14:20 降水開(kāi)始之前,850 hPa 以下低層比濕超過(guò)12 g·kg-1,17:00 之前700 hPa 以下比濕基本均大于9 g·kg-1。從14:00 開(kāi)始,近地面水汽通量最大值達(dá)24 g·hPa-1·cm-1·s-1,700 hPa 水汽通量達(dá)6 g·hPa-1·cm-1·s-1,表明由冷鋒造成的暖濕空氣強(qiáng)烈抬升階段水汽輸送旺盛,水汽含量充沛,水汽聚集時(shí)間早于降水出現(xiàn)時(shí)間,近地面水汽通量大于21 g·hPa-1·cm-1·s-1的大值存在時(shí)間與較強(qiáng)降水出現(xiàn)時(shí)間一致。

    圖4(b)顯示,22日14:00 后邢臺(tái)站近地層偏東風(fēng)轉(zhuǎn)為東北風(fēng),表明地面冷鋒移過(guò)邢臺(tái)站,降水出現(xiàn)在地面冷鋒過(guò)境之后。14:20—16:10,近地層?xùn)|北風(fēng)風(fēng)速達(dá)10 m·s-1,低層強(qiáng)勁的冷墊勢(shì)力有利于輻合抬升;700 hPa 為偏東風(fēng),風(fēng)速僅1~4 m·s-1,說(shuō)明700 hPa 切變線仍未過(guò)境;500 hPa 持續(xù)為西南氣流,風(fēng)速最大不超過(guò)6 m·s-1。700 hPa 與500 hPa 中高層引導(dǎo)氣流偏弱造成云系移動(dòng)緩慢。垂直運(yùn)動(dòng)方面,由于鋒面系統(tǒng)的輻合抬升,22日13:30之后上升運(yùn)動(dòng)明顯加強(qiáng),從近地層發(fā)展到200 hPa,大值中心出現(xiàn)在350 hPa附近(最大達(dá)1.0 m·s-1),伴隨降水出現(xiàn)產(chǎn)生下沉運(yùn)動(dòng)。之后高層仍存在上升運(yùn)動(dòng)區(qū)(高度可達(dá)100 hPa),中低層上升運(yùn)動(dòng)維持,但與高層上升運(yùn)動(dòng)區(qū)之間不連續(xù)。

    從圖4(c)可知,22日13:40 之前云系只是高層冰云,包含少量冰晶和雪晶。14:00—15:00,由于高層上升運(yùn)動(dòng)的存在,冰晶混合比快速增大,大值中心出現(xiàn)在250 hPa,最大達(dá)0.4 g·kg-1;在冰晶混合比大值中心下方區(qū)域雪晶混合比隨時(shí)間逐漸增加,其大值中心(最大值0.3 g·kg-1)出現(xiàn)在450~350 hPa。15:00 之后,300~150 hPa 高度范圍內(nèi)冰晶數(shù)濃度最高達(dá)5000個(gè)·L-1,但冰晶混合比減小,雪晶混合比變化不大;在600~500 hPa 高度范圍內(nèi)、溫度為-10~0 ℃時(shí),冰晶數(shù)濃度約20 個(gè)·L-1。冰晶混合比大值中心出現(xiàn)時(shí)間、區(qū)域與冷區(qū)上升運(yùn)動(dòng)大值中心出現(xiàn)時(shí)間、區(qū)域一致。霰粒子混合比最大達(dá)2.2 g·kg-1,大值中心出現(xiàn)在溫度-10~-5 ℃之間(高度600~500 hPa),并與過(guò)冷水混合比大值區(qū)重疊。

    圖4(d)顯示,22日13:30 后,暖區(qū)出現(xiàn)云水,云水混合比最大達(dá)0.7 g·kg-1,冷區(qū)出現(xiàn)豐富的過(guò)冷水(混合比最高達(dá)0.2 g·kg-1),云水混合比大值中心與上升運(yùn)動(dòng)大值區(qū)域基本一致,充沛的過(guò)冷水有利于霰的凇附增長(zhǎng)(collection of cloud droplets by graupel,ccg)[21]。如圖5所示,600~500 hPa高度區(qū)域ccg過(guò)程在此階段的平均轉(zhuǎn)化率超過(guò)0.0006 g·kg-1·s-1,是霰碰并雪晶過(guò)程(collection of snow by graupel,csg)增長(zhǎng)率的2 倍多。從14:10 開(kāi)始有微量雨水降落到地面,到14:40雨水混合比達(dá)最大約0.9 g·kg-1,雨滴逐漸增加,但落后于霰的增長(zhǎng),說(shuō)明霰的融化過(guò)程(melting of graupel to form rain,mgr)是降水形成的主要過(guò)程。低層充沛的云水含量,使雨滴在下落過(guò)程中可以通過(guò)碰并云滴增長(zhǎng)(collection of cloud droplets by rain,ccr),650 hPa 高度附近霰的融化過(guò)程轉(zhuǎn)化率是雨滴碰并云滴過(guò)程的2 倍多,同時(shí)低層800 hPa 附近雨滴仍可以收集云滴獲得增長(zhǎng)。伴隨著雨滴增長(zhǎng),14:50 在800~700 hPa 高度出現(xiàn)下沉運(yùn)動(dòng)中心,雨滴的拖曳加劇了下沉運(yùn)動(dòng)。

    圖4 2017年5月22日13:00—17:00邢臺(tái)站各物理量的時(shí)間-高度剖面(底部黑色粗柱狀為10 min降水量,單位:mm)(a)比濕(填色區(qū),單位:g·kg-1)及水汽通量(藍(lán)色等值線,單位:g·hPa-1·cm-1·s-1),(b)水平風(fēng)場(chǎng)(風(fēng)矢,單位:m·s-1)、垂直速度(紅色等值線,單位:m·s-1)(黃色填色區(qū)為雨水混合比大于0.01 g·kg-1區(qū)域),(c)冰晶(紅色等值線)、雪(填色區(qū))、霰(藍(lán)色等值線)混合比(單位:g·kg-1),(d)云水(填色區(qū))、雨水(紅色等值線)混合比(單位:g·kg-1)、冰晶數(shù)濃度(藍(lán)色等值線,單位:個(gè)·L-1)Fig.4 Time-height cross sections of physical quantities from 13:00 BST to 17:00 BST May 22,2017 at Xingtai station(The black thick bar at bottom is 10 min rainfall,Unit:mm)(a) specific humidity (the color shaded,Unit:g·kg-1) and water vapor flux (blue isolines,Unit:g·hPa-1·cm-1·s-1),(b) horizontal wind (wind vectors,Unit:m·s-1) and vertical velocity (red isolines,Unit:m·s-1) (the yellow shaded for the mixing ratio of rain water greater than 0.01 g·kg-1),(c) mixing ratio of ice crystal (red isolines),snow (the color shaded) and graupel(blue isolines) (Unit:g·kg-1),(d) mixing ratio of cloud water (the color shaded),rain water (red isolines)(Unit:g·kg-1) and number concentration of ice crystal (blue isolines,Unit:L-1)

    圖5 邢臺(tái)站微物理過(guò)程轉(zhuǎn)化率垂直廓線(黑色線為22日13:00—17:00平均,紅色線為22日17:10至23日01:20平均)Fig.5 Mean profiles of microphysical process conversion rates at Xingtai station(the black lines averaged from 13:00 BST to 17:00 BST May 22,the red lines averaged from 17:10 BST May 22 to 01:20 BST May 23)

    鋒面抬升強(qiáng)雨期,邢臺(tái)處于地面冷鋒之后、700 hPa 槽前,850 hPa 切變線過(guò)境邢臺(tái)前低層較強(qiáng)的輻合形成一致的上升運(yùn)動(dòng),云系發(fā)展在垂直方向上體現(xiàn)出直展特征。該階段500 hPa 以上由冰、雪晶組成,800~500 hPa中層為過(guò)冷水、雪、霰組成的冰水混合層,800 hPa 以下為液態(tài)水云,云系具有“催化-供給”結(jié)構(gòu)特征,出現(xiàn)重要的人工增雨云結(jié)構(gòu)條件[4],霰的融化過(guò)程和雨滴碰并云滴過(guò)程是形成降水的主要過(guò)程,冷云降水過(guò)程和暖云降水過(guò)程都重要。

    4.2 深厚層云小雨期

    圖6(a)顯示,22日17:00至23日01:00,600 hPa以上水汽含量及水汽輸送變化不大,但700 hPa 之下中低層水汽含量、水汽通量明顯小于因鋒面抬升形成的較強(qiáng)降雨階段。從圖6(b)可知,此階段850 hPa之下為東北氣流控制,17:50 前后700 hPa 切變線過(guò)境邢臺(tái),22:00 之前500 hPa 一直為西南風(fēng),風(fēng)速仍維持在4~6 m·s-1,22:00 開(kāi)始轉(zhuǎn)為偏西風(fēng)且持續(xù)到23日01:00,23日01:00 前后500 hPa 槽線過(guò)境邢臺(tái),之后轉(zhuǎn)為西北風(fēng)。小雨階段邢臺(tái)處在700 hPa槽后、500 hPa槽前。該階段上升運(yùn)動(dòng)主要出現(xiàn)在冷區(qū),上升速度明顯減小(大多數(shù)時(shí)間上升速度僅0.1 m·s-1),只在局部時(shí)間300 hPa 高度附近出現(xiàn)0.4 m·s-1的較強(qiáng)上升運(yùn)動(dòng)中心。

    圖6(c)和圖6(d)顯示,云頂高度維持在約12 km,云底高度抬升到約2 km。300~200 hPa 高層冰晶混合比及數(shù)濃度明顯減小,冰晶數(shù)濃度最多約100個(gè)·L-1。在400~300 hPa之間出現(xiàn)雪晶混合比大值中心,最大值約0.7 g·kg-1,雪晶混合比高于鋒面抬升造成的較強(qiáng)降雨階段,但霰粒子混合比明顯減小,最大僅0.3 g·kg-1,其大值中心出現(xiàn)在雪晶混合比大值中心的下方區(qū)域,霰碰并雪晶過(guò)程(csg)是霰增長(zhǎng)的主要過(guò)程。云水僅僅出現(xiàn)在800~500 hPa 之間,云水混合比隨時(shí)間推移逐漸減小,云水區(qū)域與上升運(yùn)動(dòng)區(qū)域密切相關(guān),22:00之前-10~0 ℃之間依然存在過(guò)冷水,但混合比減小,最大約0.1 g·kg-1。雨水隨高度下降而減少,表明低層不飽和,雨滴在下落過(guò)程中有蒸發(fā)現(xiàn)象,降水主要由霰的融化過(guò)程(mgr)產(chǎn)生,以冷云降水過(guò)程為主。

    圖6 2017年5月22日17:00至23日01:00邢臺(tái)站各物理量的時(shí)間-高度剖面(底部黑色粗柱狀為10 min降水量,單位:mm)(a)比濕(填色區(qū),單位:g·kg-1)及水汽通量(藍(lán)色等值線,單位:g·hPa-1·cm-1·s-1),(b)水平風(fēng)場(chǎng)(風(fēng)矢,單位:m·s-1)、垂直速度(紅色等值線,單位:m·s-1)(黃色填色為雨水混合比大于0.01 g·kg-1區(qū)域),(c)冰晶(紅色等值線)、雪(填色區(qū))、霰(藍(lán)色等值線)混合比(單位:g·kg-1),(d)云水(填色區(qū))、雨水(紅色等值線)混合比(單位:g·kg-1)及冰晶數(shù)濃度(藍(lán)色等值線,單位:個(gè)·L-1)Fig.6 Time-height cross sections of physical quantities from 17:00 BST May 22 to 01:00 BST May 23 at Xingtai station(The black thick bar at bottom is 10 min rainfall,Unit:mm)(a) specific humidity (the color shaded,Unit:g·kg-1) and water vapor flux (blue isolines,Unit:g·hPa-1·cm-1·s-1),(b) horizontal wind (wind vectors,Unit:m·s-1) and vertical velocity (red isolines,Unit:m·s-1) (the yellow shaded for the mixing ratio of rain water greater than 0.01 g·kg-1),(c) mixing ratio of ice crystal (red isolines),snow (the color shaded)and graupel (blue isolines) (Unit:g·kg-1),(d) mixing ratio of cloud water (the color shaded),rain water(red isolines) (Unit:g·kg-1) and number concentration of ice crystal (blue isolines,Unit:L-1)

    4.3 消散期

    當(dāng)500 hPa 槽線從邢臺(tái)上空過(guò)境(23日01:00)后,邢臺(tái)站逐漸由小雨轉(zhuǎn)為毛毛雨,降水過(guò)程至23日04:00停止。

    此階段邢臺(tái)站700 hPa以上水汽含量明顯降低。高層沒(méi)有冰相水凝物,云團(tuán)僅出現(xiàn)在700 hPa 以下暖區(qū),云頂高度小于3 km,云底高度約1 km,云水混合比高值區(qū)位于800 hPa 附近,最高達(dá)0.3 g·kg-1(圖略),推知消散階段的毛毛雨由暖云降水過(guò)程產(chǎn)生。

    5 增雨條件分析

    上述分析顯示,不同降水階段邢臺(tái)站上空云系結(jié)構(gòu)差異顯著,那么云系的增雨條件也必有差別。

    孫晶等[7]在研究西風(fēng)槽降水云系的增雨條件時(shí)總結(jié)了一套指標(biāo)體系,認(rèn)為當(dāng)同時(shí)滿足雨強(qiáng)大于0.1 mm·h-1、上升運(yùn)動(dòng)區(qū)(垂直速度大于0.01 m·s-1)、過(guò)冷水含量大于0.01 g·kg-1、冰晶數(shù)濃度大于20個(gè)·L-1、溫度為-20~0 ℃條件時(shí),云層具有可播性;若冰晶數(shù)濃度小于20 個(gè)·L-1時(shí),云層具有強(qiáng)可播性。但該指標(biāo)體系沒(méi)有考慮云中冰面過(guò)飽和水汽對(duì)增雨的潛力,在層狀云降水形成過(guò)程中,冰面過(guò)飽和水汽對(duì)凝華增長(zhǎng)有較大作用[4]。當(dāng)云中有過(guò)冷水時(shí)引入AgI 將會(huì)通過(guò)貝吉龍過(guò)程促進(jìn)冰晶增長(zhǎng),在沒(méi)有過(guò)冷水的區(qū)域,AgI 以凝華核化為主,冰面飽和比大于1.04是AgI凝華核化的必要條件[15]。

    在已有研究[7,9,11-13]基礎(chǔ)上,本文增加冰面飽和比作為增雨作業(yè)條件的指標(biāo)之一。選取雨強(qiáng)、上升運(yùn)動(dòng)區(qū)、溫度、過(guò)冷水混合比、冰面飽和比與冰晶數(shù)濃度6 個(gè)物理因子組成增雨作業(yè)條件指標(biāo),這組指標(biāo)在宏觀方面重點(diǎn)考慮了云系的發(fā)展階段、云中是否存在上升運(yùn)動(dòng)、是否產(chǎn)生降水,微觀方面主要考慮過(guò)冷水混合比、冰面飽和比、冰晶數(shù)濃度3 個(gè)AgI 核化增長(zhǎng)的關(guān)鍵物理因子。具體閾值設(shè)計(jì)如下:當(dāng)同時(shí)滿足雨強(qiáng)大于0.1 mm·h-1、上升速度大于0.01 m·s-1、溫度為-20~0 ℃、過(guò)冷水混合比大于0.01 g·kg-1或冰面飽和比大于1.04、冰晶數(shù)濃度大于20個(gè)·L-1條件時(shí),云系具有可播性;若冰晶數(shù)濃度小于20個(gè)·L-1時(shí),云系具有強(qiáng)可播性。

    利用構(gòu)建的增雨條件指標(biāo),對(duì)經(jīng)過(guò)邢臺(tái)站的此次低槽冷鋒云系的增雨時(shí)機(jī)和部位進(jìn)行識(shí)別。由冰晶、雪和霰混合比求和得到冰相水凝物的混合比,代表冰云分布。圖7(a)顯示,在鋒面抬升強(qiáng)雨期,增雨作業(yè)高度主要位于4.0~7.9 km 之間,與600 hPa高度以上過(guò)冷水出現(xiàn)的區(qū)域基本相同,其中4.0~5.5 km 高度層內(nèi)為強(qiáng)可播區(qū),結(jié)合圖4可知,該高度層內(nèi)過(guò)冷水混合比最高可達(dá)0.2 g·kg-1,冰晶數(shù)濃度不足20個(gè)·L-1,上升速度最高可達(dá)0.5 m·s-1;從增雨作業(yè)時(shí)機(jī)角度分析,地面降水出現(xiàn)之后的較強(qiáng)降雨期都有增雨作業(yè)時(shí)機(jī),而在22日14:00 前,云系以冰云為主,在-20~0 ℃云層中缺乏過(guò)冷水且冰面飽和比小于1.04,云系沒(méi)有增雨作業(yè)時(shí)機(jī)。深厚層云小雨期[圖7(b)],此階段可播區(qū)仍主要位于4.0~7.9 km 之間,強(qiáng)可播區(qū)出現(xiàn)在4.0~5.5 km 高度層內(nèi),與鋒面抬升強(qiáng)雨期基本一致;增雨時(shí)機(jī)出現(xiàn)在600 hPa 維持西南風(fēng)階段,當(dāng)600 hPa 轉(zhuǎn)為西北風(fēng)影響后,盡管雨水混合比大于0.01 g·kg-1,表明地面仍產(chǎn)生降水,但由于云中沒(méi)有過(guò)冷水且冰面飽和比小于1.04,云中沒(méi)有增雨作業(yè)條件。

    圖7 2017年5月22日13:00—17:00邢臺(tái)站鋒面抬升強(qiáng)雨期(a)及22日17:00至23日01:00深厚層云小雨期(b)物理量變化與可播區(qū)識(shí)別Fig.7 Variation of physical quantities and the identified glaciogenic seeding areas during the frontal uplifting heavy rain period from 13:00 BST to 17:00 BST 22 May 2017 (a) and the deep stratus cloud light rain period from 17:00 BST 22 to 01:00 BST 23 May 2017 (b) at Xingtai station

    結(jié)合圖4和圖6可知,當(dāng)云系處于鋒面抬升較強(qiáng)發(fā)展階段時(shí)增雨作業(yè)可主要在冰水混合區(qū)開(kāi)展,而當(dāng)云系處于深厚層云階段時(shí),由于云中上升速度減小且上升運(yùn)動(dòng)主要出現(xiàn)在冷區(qū),此階段過(guò)冷水含量減少且出現(xiàn)高度比鋒面抬升強(qiáng)雨期有所下降,增雨作業(yè)除了可在冰水混合區(qū)開(kāi)展之外,還可以在其上大約3 km 高度范圍內(nèi)冰面飽和比大于1.04 的區(qū)域開(kāi)展。

    總體看來(lái),邢臺(tái)站增雨作業(yè)時(shí)機(jī)主要出現(xiàn)在鋒面抬升形成的較強(qiáng)降雨期及深厚層云形成的小雨期,降水消散期云中沒(méi)有可播區(qū)。出現(xiàn)可播區(qū)時(shí)云系主要是處于發(fā)展及維持階段的冷暖混合云,云頂高度達(dá)12 km,云底高度維持在近地面到2 km 之間,云系冷層厚度約8 km,冰水混合層厚度約4 km,高層有冰晶和雪,中間層為含有過(guò)冷水的冰水混合區(qū),低層為液態(tài)水云,具有“催化-供給”云系結(jié)構(gòu)特征。

    6 結(jié) 論

    華北地區(qū)2017年春季出現(xiàn)較嚴(yán)重干旱,5月22—23日出現(xiàn)一次明顯降水過(guò)程,河北、山西等地氣象部門開(kāi)展飛機(jī)和地面增雨作業(yè)。本文利用CPEFS_v1.0系統(tǒng),針對(duì)此次低槽冷鋒降水過(guò)程,以河北邢臺(tái)站為代表研究了不同降水階段云結(jié)構(gòu)特征及增雨作業(yè)條件,得到如下結(jié)論:

    (1)低槽冷鋒云降水系統(tǒng)自西北向東南移動(dòng)經(jīng)過(guò)邢臺(tái)站時(shí),先出現(xiàn)10 min 雨量為1.0~3.0 mm 且持續(xù)約3 h的較強(qiáng)降水期,之后出現(xiàn)10 min雨量為0.2~0.7 mm 且持續(xù)約8 h 的小雨期,之后降水消散,模擬結(jié)果與實(shí)況一致。

    (2)因鋒面抬升形成的較強(qiáng)降雨階段,邢臺(tái)處于地面冷鋒之后、700 hPa 槽前深厚的西南氣流影響下,此階段與850 hPa 以下水汽通量大值(大于21 g·hPa-1·cm-1·s-1)存在時(shí)間一致,850 hPa 切變線過(guò)境邢臺(tái)前后云中出現(xiàn)上下貫穿的上升運(yùn)動(dòng)。500 hPa 以上高層冰云由冰晶和雪晶組成;800~500 hPa中層云為冰水混合云,過(guò)冷水與霰粒子混合比最高分別達(dá)0.2、2.2 g·kg-1;800 hPa 以下低層云為液態(tài)水云,云系具有“催化-供給”云結(jié)構(gòu)特征,霰的融化過(guò)程、雨滴碰并云滴過(guò)程是降水形成的主要過(guò)程。

    (3)在22日17:10 至23日01:20 深厚層云形成的小雨階段,邢臺(tái)處于700 hPa 槽后和500 hPa 槽前。700 hPa 之下中低層水汽含量及水汽通量明顯減小,上升運(yùn)動(dòng)主要出現(xiàn)在冷區(qū),上升速度最大僅0.4 m·s-1,云系依舊為冷暖混合云,冰晶、霰混合比明顯減小,雪晶混合比增加,-10~0 ℃之間過(guò)冷水混合比減小,最大約0.1 g·kg-1,低層云水混合比也減小,降水主要由霰的融化過(guò)程產(chǎn)生。

    (4)綜合利用逐時(shí)雨強(qiáng)、上升運(yùn)動(dòng)區(qū)、溫度、冰面飽和比、過(guò)冷水含量、冰晶數(shù)濃度6個(gè)物理量構(gòu)建增雨條件指標(biāo),發(fā)現(xiàn)本次邢臺(tái)地區(qū)低槽冷鋒云系適宜冷云催化增雨作業(yè)的時(shí)間主要出現(xiàn)在鋒面抬升強(qiáng)雨期與深厚層云小雨期,這時(shí)云系為處于發(fā)展及維持階段的冷暖混合云,可播區(qū)高度位于4.0~7.9 km,強(qiáng)可播區(qū)出現(xiàn)在4.0~5.5 km 高度過(guò)冷水充沛區(qū)域。

    猜你喜歡
    云系鋒面冷鋒
    四川盆地冷鋒特征及其對(duì)逆溫和大氣污染的影響
    熱聲耦合燃燒振蕩中火焰鋒面識(shí)別分析
    煤氣與熱力(2022年2期)2022-03-09 06:29:20
    2019年夏季長(zhǎng)江口及鄰近海域鋒面控制下葉綠素a的分布特征及其環(huán)境影響因素分析
    2020年江西汛期大暴雨衛(wèi)星云圖特征分析
    2019年5月26日朝陽(yáng)飛機(jī)人工增雨作業(yè)分析
    基于核心素養(yǎng)的高中地理“問(wèn)題式教學(xué)”——以“鋒面氣旋”為例
    哈爾濱地區(qū)冷鋒分析研究
    綠色科技(2018年22期)2019-01-15 05:36:08
    電影《戰(zhàn)狼2》中冷鋒人物形象及其現(xiàn)實(shí)意義
    廣西11—12月人工增雨天氣研究
    一次低槽冷鋒暴雨數(shù)值模擬和診斷分析
    亚洲av电影在线观看一区二区三区 | 亚洲高清免费不卡视频| 亚洲精品中文字幕在线视频 | 日本色播在线视频| av在线天堂中文字幕| 中文乱码字字幕精品一区二区三区 | 神马国产精品三级电影在线观看| 日韩av免费高清视频| 少妇猛男粗大的猛烈进出视频 | 内射极品少妇av片p| av在线蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 超碰av人人做人人爽久久| 欧美xxxx性猛交bbbb| 欧美一级a爱片免费观看看| 岛国毛片在线播放| 婷婷六月久久综合丁香| 国产免费视频播放在线视频 | 美女脱内裤让男人舔精品视频| 国产伦理片在线播放av一区| 亚洲欧美成人精品一区二区| 欧美不卡视频在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 国语对白做爰xxxⅹ性视频网站| 精品久久久精品久久久| 在线免费观看不下载黄p国产| 麻豆乱淫一区二区| 超碰97精品在线观看| 成年av动漫网址| 午夜久久久久精精品| 精品国内亚洲2022精品成人| a级一级毛片免费在线观看| 少妇熟女aⅴ在线视频| 三级男女做爰猛烈吃奶摸视频| 精品不卡国产一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产白丝娇喘喷水9色精品| 全区人妻精品视频| 国产爱豆传媒在线观看| 免费高清在线观看视频在线观看| 日韩精品青青久久久久久| 黄片无遮挡物在线观看| 又大又黄又爽视频免费| 久久久久久久久久黄片| 色综合站精品国产| 日本爱情动作片www.在线观看| 好男人视频免费观看在线| 高清欧美精品videossex| 亚洲成人久久爱视频| 超碰97精品在线观看| 亚洲av成人精品一区久久| 欧美激情国产日韩精品一区| 免费观看av网站的网址| 午夜福利在线在线| 国产高清有码在线观看视频| 尤物成人国产欧美一区二区三区| 成人毛片60女人毛片免费| 亚洲av福利一区| 午夜精品在线福利| 午夜精品一区二区三区免费看| 成年女人看的毛片在线观看| 日本猛色少妇xxxxx猛交久久| 麻豆成人午夜福利视频| 美女被艹到高潮喷水动态| 永久免费av网站大全| 亚洲av成人av| 韩国高清视频一区二区三区| 精品酒店卫生间| 欧美bdsm另类| 亚洲欧美清纯卡通| 99久久精品一区二区三区| 国产爱豆传媒在线观看| 最近中文字幕高清免费大全6| 欧美 日韩 精品 国产| 国产探花极品一区二区| 国产成人免费观看mmmm| 国产美女午夜福利| 亚洲在线自拍视频| 国产精品蜜桃在线观看| 午夜日本视频在线| 久久99蜜桃精品久久| 国产成人免费观看mmmm| 99久久中文字幕三级久久日本| 99久久精品热视频| 少妇裸体淫交视频免费看高清| 一区二区三区免费毛片| 国产在视频线精品| 亚洲国产色片| 晚上一个人看的免费电影| 精品一区二区三卡| 免费看光身美女| 大香蕉久久网| 国产精品99久久久久久久久| 亚洲成人av在线免费| 亚洲在久久综合| 大又大粗又爽又黄少妇毛片口| 日韩电影二区| 亚洲最大成人中文| 永久网站在线| 免费播放大片免费观看视频在线观看| 国产一区亚洲一区在线观看| 国产高清有码在线观看视频| 在线观看免费高清a一片| 男女国产视频网站| 亚洲欧美成人综合另类久久久| 午夜福利在线在线| 最近最新中文字幕免费大全7| 一级二级三级毛片免费看| 又爽又黄无遮挡网站| 黄色配什么色好看| 精品午夜福利在线看| 丰满少妇做爰视频| 最近的中文字幕免费完整| 水蜜桃什么品种好| 午夜激情欧美在线| 久久久久久久大尺度免费视频| 国产免费又黄又爽又色| 99视频精品全部免费 在线| 黑人高潮一二区| 好男人在线观看高清免费视频| 成人午夜精彩视频在线观看| 日韩精品有码人妻一区| 好男人在线观看高清免费视频| 看黄色毛片网站| 亚洲欧美精品专区久久| 亚洲av成人av| 亚洲精品色激情综合| av免费在线看不卡| 最近最新中文字幕免费大全7| 久久久久久久久久黄片| 黑人高潮一二区| 联通29元200g的流量卡| 国产精品99久久久久久久久| 国产成人精品久久久久久| 男的添女的下面高潮视频| 少妇被粗大猛烈的视频| 久久精品综合一区二区三区| 国产成人a∨麻豆精品| 亚洲精品成人av观看孕妇| 又黄又爽又刺激的免费视频.| 久久久久久久久久成人| 成人毛片60女人毛片免费| 成人国产麻豆网| 国产高清国产精品国产三级 | 国产精品一区二区在线观看99 | 亚洲国产成人一精品久久久| 永久网站在线| 一级二级三级毛片免费看| 日韩一区二区三区影片| 久久久久久久大尺度免费视频| 久久草成人影院| 天堂影院成人在线观看| 精品少妇黑人巨大在线播放| 我的老师免费观看完整版| 美女国产视频在线观看| 69人妻影院| 三级毛片av免费| 国产在视频线在精品| 伊人久久精品亚洲午夜| 日韩av不卡免费在线播放| 国产白丝娇喘喷水9色精品| 婷婷色麻豆天堂久久| 最近最新中文字幕免费大全7| 最近中文字幕高清免费大全6| 又爽又黄a免费视频| 亚洲精品aⅴ在线观看| av在线天堂中文字幕| 别揉我奶头 嗯啊视频| 人妻夜夜爽99麻豆av| 舔av片在线| 亚洲av电影不卡..在线观看| 91午夜精品亚洲一区二区三区| 久久精品久久精品一区二区三区| 亚洲丝袜综合中文字幕| 亚洲av成人精品一二三区| 男女视频在线观看网站免费| 久久国内精品自在自线图片| 中文字幕av成人在线电影| 少妇被粗大猛烈的视频| 日韩av免费高清视频| 亚洲aⅴ乱码一区二区在线播放| 熟女人妻精品中文字幕| 天堂影院成人在线观看| 国产伦理片在线播放av一区| 久久久久久九九精品二区国产| 日韩大片免费观看网站| 在线观看免费高清a一片| 国产一级毛片在线| 一级毛片aaaaaa免费看小| 大香蕉久久网| 国产免费福利视频在线观看| 少妇丰满av| 一夜夜www| 51国产日韩欧美| 99热这里只有是精品50| 国产女主播在线喷水免费视频网站 | av在线天堂中文字幕| 国产视频内射| 91精品一卡2卡3卡4卡| 2022亚洲国产成人精品| 亚洲av一区综合| 蜜桃亚洲精品一区二区三区| 九草在线视频观看| 日韩制服骚丝袜av| 久久久色成人| 成人午夜精彩视频在线观看| av播播在线观看一区| 成人午夜高清在线视频| 又爽又黄a免费视频| 纵有疾风起免费观看全集完整版 | 在线播放无遮挡| 中文天堂在线官网| 久久久久久久午夜电影| 亚洲精品第二区| 国产伦在线观看视频一区| 男女啪啪激烈高潮av片| 国产精品爽爽va在线观看网站| 成人二区视频| 亚洲在线自拍视频| 成人毛片a级毛片在线播放| 国国产精品蜜臀av免费| 欧美xxxx黑人xx丫x性爽| 伊人久久国产一区二区| 寂寞人妻少妇视频99o| 97超视频在线观看视频| 高清欧美精品videossex| 精品国内亚洲2022精品成人| 亚洲欧美精品自产自拍| 最近最新中文字幕免费大全7| 午夜老司机福利剧场| 成人鲁丝片一二三区免费| 九草在线视频观看| 日本午夜av视频| 精品久久久久久久末码| 观看免费一级毛片| 国产精品久久久久久久电影| 一级毛片 在线播放| 久久午夜福利片| 久久久久久久亚洲中文字幕| 久久精品人妻少妇| 少妇的逼好多水| 亚洲精品日本国产第一区| 80岁老熟妇乱子伦牲交| 91久久精品电影网| 精品一区在线观看国产| 亚洲一区高清亚洲精品| 日本熟妇午夜| 久久这里只有精品中国| 免费观看无遮挡的男女| 亚洲精品久久久久久婷婷小说| 精品久久久久久久久av| 午夜激情福利司机影院| 亚洲性久久影院| 亚洲欧美成人综合另类久久久| 久久久久九九精品影院| 亚洲天堂国产精品一区在线| 亚洲激情五月婷婷啪啪| 床上黄色一级片| 非洲黑人性xxxx精品又粗又长| 91狼人影院| 26uuu在线亚洲综合色| 18禁在线无遮挡免费观看视频| 亚洲精品456在线播放app| 黄色日韩在线| 亚洲成人中文字幕在线播放| 一级二级三级毛片免费看| 极品教师在线视频| 伦理电影大哥的女人| 久久这里只有精品中国| 2021少妇久久久久久久久久久| 国产色婷婷99| 男的添女的下面高潮视频| 国产国拍精品亚洲av在线观看| 日日摸夜夜添夜夜爱| 亚洲精品第二区| 久久久久免费精品人妻一区二区| 国产精品综合久久久久久久免费| 国产精品一区二区三区四区久久| 日本一二三区视频观看| 亚洲图色成人| 色综合色国产| 在线观看av片永久免费下载| 中文精品一卡2卡3卡4更新| av国产久精品久网站免费入址| 老女人水多毛片| www.av在线官网国产| 婷婷色麻豆天堂久久| 狂野欧美激情性xxxx在线观看| 在线免费观看不下载黄p国产| 欧美人与善性xxx| www.色视频.com| 大香蕉97超碰在线| 久久久久久久久久成人| 国产成人91sexporn| 成人av在线播放网站| 中文资源天堂在线| 国产精品一及| 亚洲精品第二区| 欧美+日韩+精品| av国产免费在线观看| 美女xxoo啪啪120秒动态图| 欧美高清性xxxxhd video| 国产一级毛片七仙女欲春2| 国产黄a三级三级三级人| 国产人妻一区二区三区在| 18+在线观看网站| 亚洲国产精品sss在线观看| 亚洲aⅴ乱码一区二区在线播放| 久久99热6这里只有精品| 国产精品不卡视频一区二区| 国产午夜精品久久久久久一区二区三区| 一级二级三级毛片免费看| 晚上一个人看的免费电影| 久久久久久久久久人人人人人人| 亚洲av.av天堂| 免费黄色在线免费观看| 日韩国内少妇激情av| 我的女老师完整版在线观看| 2021少妇久久久久久久久久久| 国产伦一二天堂av在线观看| 久久99精品国语久久久| 久久久精品94久久精品| 亚洲成人久久爱视频| 中文资源天堂在线| 久久久久精品久久久久真实原创| 欧美zozozo另类| 真实男女啪啪啪动态图| 亚洲激情五月婷婷啪啪| 在线观看免费高清a一片| 如何舔出高潮| 亚洲成人av在线免费| 九色成人免费人妻av| 亚洲精品成人av观看孕妇| 18+在线观看网站| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说 | 日韩国内少妇激情av| 日韩欧美三级三区| 精品人妻熟女av久视频| 欧美变态另类bdsm刘玥| 国产 一区精品| 99久久精品国产国产毛片| 亚洲最大成人中文| 国产精品久久久久久精品电影| 亚洲av免费高清在线观看| 亚洲在线观看片| 永久网站在线| 成人特级av手机在线观看| 亚洲欧美一区二区三区国产| 亚洲最大成人中文| 亚洲自偷自拍三级| 色视频www国产| 国产精品1区2区在线观看.| 国产 亚洲一区二区三区 | 亚洲精品影视一区二区三区av| 丝袜喷水一区| 国产熟女欧美一区二区| 神马国产精品三级电影在线观看| 少妇的逼好多水| 久久精品夜色国产| 日韩在线高清观看一区二区三区| 中文字幕av成人在线电影| 免费看美女性在线毛片视频| 亚洲av电影不卡..在线观看| 夫妻性生交免费视频一级片| 深夜a级毛片| 丰满乱子伦码专区| 国产精品国产三级国产专区5o| 特级一级黄色大片| 日日撸夜夜添| 人妻少妇偷人精品九色| 一个人看视频在线观看www免费| 免费播放大片免费观看视频在线观看| 亚洲国产色片| 天天一区二区日本电影三级| 国产欧美另类精品又又久久亚洲欧美| 成年免费大片在线观看| 成人亚洲欧美一区二区av| xxx大片免费视频| 51国产日韩欧美| 白带黄色成豆腐渣| 五月玫瑰六月丁香| 成人二区视频| 久久久精品欧美日韩精品| 亚洲综合精品二区| 伦精品一区二区三区| 国产精品国产三级国产av玫瑰| 国产成人aa在线观看| 日本熟妇午夜| 高清av免费在线| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| 国产免费视频播放在线视频 | 亚洲国产精品专区欧美| 最近的中文字幕免费完整| h日本视频在线播放| 男女国产视频网站| 晚上一个人看的免费电影| 亚洲在久久综合| 草草在线视频免费看| 亚洲精品成人av观看孕妇| 日本熟妇午夜| 色综合亚洲欧美另类图片| 亚洲四区av| 国产成人午夜福利电影在线观看| 91aial.com中文字幕在线观看| 80岁老熟妇乱子伦牲交| 欧美日韩在线观看h| 好男人视频免费观看在线| 国产色婷婷99| 国产三级在线视频| 午夜免费激情av| 国产欧美日韩精品一区二区| 男女边摸边吃奶| 一夜夜www| 亚洲国产精品国产精品| 日产精品乱码卡一卡2卡三| 简卡轻食公司| 欧美+日韩+精品| 精品久久国产蜜桃| 狠狠精品人妻久久久久久综合| 大陆偷拍与自拍| 五月玫瑰六月丁香| 秋霞在线观看毛片| 欧美不卡视频在线免费观看| 久久久久久久久久黄片| 日本黄色片子视频| 久久久久久久久久人人人人人人| 久久人人爽人人爽人人片va| 久久久久久久久久久丰满| 我的老师免费观看完整版| 日韩成人av中文字幕在线观看| 欧美一级a爱片免费观看看| 男人舔女人下体高潮全视频| 国产69精品久久久久777片| 男人舔奶头视频| 高清欧美精品videossex| 少妇高潮的动态图| 一级毛片电影观看| 真实男女啪啪啪动态图| eeuss影院久久| 最近的中文字幕免费完整| 国产av码专区亚洲av| 久久6这里有精品| 午夜免费男女啪啪视频观看| 日日撸夜夜添| 99久国产av精品| 国产爱豆传媒在线观看| 97热精品久久久久久| 日韩欧美一区视频在线观看 | 草草在线视频免费看| 国产精品蜜桃在线观看| 久久99蜜桃精品久久| 午夜激情欧美在线| av.在线天堂| 91精品伊人久久大香线蕉| 99re6热这里在线精品视频| 久久人人爽人人爽人人片va| av线在线观看网站| 国产 一区精品| 99热6这里只有精品| 麻豆久久精品国产亚洲av| 成人无遮挡网站| 亚洲最大成人中文| 国产 一区精品| 插阴视频在线观看视频| 高清视频免费观看一区二区 | 久久久久久久国产电影| 久久久久久九九精品二区国产| 一个人免费在线观看电影| 国产亚洲av片在线观看秒播厂 | 校园人妻丝袜中文字幕| 久久精品熟女亚洲av麻豆精品 | 欧美3d第一页| 国产极品天堂在线| 国产美女午夜福利| 纵有疾风起免费观看全集完整版 | 在线观看人妻少妇| av国产免费在线观看| 精品久久久久久成人av| 日韩精品有码人妻一区| 亚洲图色成人| 久久精品久久久久久噜噜老黄| 麻豆av噜噜一区二区三区| 亚洲激情五月婷婷啪啪| 国产乱来视频区| 亚洲真实伦在线观看| 久久精品国产亚洲网站| 久久精品人妻少妇| 成人av在线播放网站| 插逼视频在线观看| 男女视频在线观看网站免费| 亚洲最大成人中文| 亚洲精品久久午夜乱码| 午夜免费男女啪啪视频观看| av又黄又爽大尺度在线免费看| 内地一区二区视频在线| 亚洲自拍偷在线| 国产精品国产三级国产专区5o| 中文字幕av成人在线电影| 女人久久www免费人成看片| 偷拍熟女少妇极品色| 国产激情偷乱视频一区二区| 亚洲国产欧美在线一区| 欧美潮喷喷水| 男女国产视频网站| 成人一区二区视频在线观看| 亚洲欧洲国产日韩| 免费黄网站久久成人精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产成人一精品久久久| 少妇人妻一区二区三区视频| 久久久久久久国产电影| 哪个播放器可以免费观看大片| 最近的中文字幕免费完整| 日本免费在线观看一区| 亚洲高清免费不卡视频| 老司机影院毛片| 毛片女人毛片| 一区二区三区高清视频在线| 国产高清三级在线| 国产av在哪里看| 人妻少妇偷人精品九色| 国产免费福利视频在线观看| 99久国产av精品国产电影| 在线天堂最新版资源| 久久国内精品自在自线图片| 国产精品久久久久久久电影| 国产综合精华液| 人妻一区二区av| 一级毛片久久久久久久久女| 日韩电影二区| 美女国产视频在线观看| 在线观看人妻少妇| 日本猛色少妇xxxxx猛交久久| 波多野结衣巨乳人妻| 国产黄a三级三级三级人| 人人妻人人看人人澡| 亚洲精品色激情综合| 精品国内亚洲2022精品成人| 乱人视频在线观看| 晚上一个人看的免费电影| eeuss影院久久| 精品亚洲乱码少妇综合久久| 国产一区亚洲一区在线观看| 天堂av国产一区二区熟女人妻| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 2018国产大陆天天弄谢| a级一级毛片免费在线观看| 色综合站精品国产| 成人美女网站在线观看视频| 亚洲国产欧美在线一区| 国产日韩欧美在线精品| 久久人人爽人人爽人人片va| 赤兔流量卡办理| 超碰av人人做人人爽久久| 搞女人的毛片| 嫩草影院新地址| av福利片在线观看| 亚洲av免费高清在线观看| 中文字幕av成人在线电影| 麻豆av噜噜一区二区三区| 久久久久精品久久久久真实原创| 久久久久久久久久久丰满| av免费在线看不卡| 老师上课跳d突然被开到最大视频| 亚洲国产av新网站| 亚洲精品一区蜜桃| 国内揄拍国产精品人妻在线| 国产69精品久久久久777片| 黄片wwwwww| 国产在线一区二区三区精| 色综合亚洲欧美另类图片| 青春草视频在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美激情在线99| 天堂中文最新版在线下载 | 亚洲av免费在线观看| 成人综合一区亚洲| 男女边吃奶边做爰视频| 中国美白少妇内射xxxbb| 国产高清有码在线观看视频| 国产av码专区亚洲av| 欧美成人一区二区免费高清观看| 成人亚洲欧美一区二区av| 亚洲四区av| 久久99精品国语久久久| 99热这里只有精品一区| 免费不卡的大黄色大毛片视频在线观看 | 欧美三级亚洲精品| 亚洲最大成人手机在线| 国产黄频视频在线观看| 成年av动漫网址| 777米奇影视久久| 人妻一区二区av| 午夜福利视频精品| 日韩欧美 国产精品| 亚洲精品456在线播放app| 少妇被粗大猛烈的视频| 午夜精品国产一区二区电影 | 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 日本wwww免费看| 成人一区二区视频在线观看| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 久久久久网色| 久久久午夜欧美精品| 国产精品蜜桃在线观看| 国产精品人妻久久久影院| 国产淫语在线视频| 免费大片18禁| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片|