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

    遼寧開原龍卷強(qiáng)對(duì)流過程的擾動(dòng)天氣環(huán)境

    2021-05-07 13:06:24錢維宏艾陽陳笑晨
    地球物理學(xué)報(bào) 2021年5期
    關(guān)鍵詞:龍卷強(qiáng)對(duì)流氣流

    錢維宏, 艾陽, 陳笑晨

    1 北京大學(xué)大氣與海洋科學(xué)系, 北京 100871 2 中國氣象局廣州熱帶海洋氣象研究所/廣東省區(qū)域數(shù)值天氣預(yù)報(bào)重點(diǎn)實(shí)驗(yàn)室, 廣州 510641 3 廈門市海峽氣象開放重點(diǎn)實(shí)驗(yàn)室, 廈門 361013 4 福建省氣候中心, 福州 350001

    0 引言

    夏季,我國東北常有高空冷渦活動(dòng)(孫力等,1994),是有可能形成局地強(qiáng)對(duì)流的最大母體天氣系統(tǒng).在常規(guī)天氣圖上,上部對(duì)流層冷渦西側(cè)的橫槽,會(huì)形成冷空氣的堆積,一旦橫槽轉(zhuǎn)豎,冷空氣南下與低層大氣中的暖濕氣團(tuán)交匯形成層結(jié)不穩(wěn)定,其不穩(wěn)定能量的釋放會(huì)形成強(qiáng)對(duì)流天氣.在東北冷渦控制的背景下,冷渦南側(cè)的次天氣尺度短波槽有可能誘發(fā)龍卷(王寧等,2014;王婷婷等,2017).同時(shí),有研究表明700~500 hPa之間較大的溫差和較強(qiáng)的低空(850~925 hPa)急流都有利于龍卷的發(fā)生(王秀明等,2015).

    我國夏季龍卷風(fēng)發(fā)生頻次最高的區(qū)域是江蘇的北部(范雯杰和俞小鼎,2015).蘇北里下河水網(wǎng)地與其北側(cè)黃河故道沙臺(tái)高地,和其東側(cè)黃海沖積臺(tái)地之間存在著下墊面的屬性差異.前者形成的氣流輻合線沿蘇北灌溉渠呈東西走向,而后者形成的海陸風(fēng)輻合線平行于海岸線呈東南至西北走向,蘇北的龍卷風(fēng)多數(shù)發(fā)生在上述兩條近地面氣流輻合帶上(錢維宏等,2017).這樣的氣候環(huán)境條件一旦遇上中尺度渦旋,龍卷和冰雹等強(qiáng)對(duì)流天氣就有可能發(fā)生.我國其他地區(qū)發(fā)生的龍卷并不具備蘇北這樣有利的氣候輻合條件,所以龍卷的發(fā)生地理位置相對(duì)分散.開原位于我國東北平原的沈陽與長春兩個(gè)省會(huì)城市之間,東西兩側(cè)是山地,西南暖濕氣流容易從低層進(jìn)入東北平原谷地,遇上對(duì)流層中上部冷渦活動(dòng)時(shí),對(duì)流天氣極易發(fā)生.

    在世界上,龍卷發(fā)生年頻次最高的地區(qū)是美國,2010年前的30年中平均每年發(fā)生1100多個(gè)龍卷.目前,龍卷的預(yù)警仍主要依靠雷達(dá)探測(cè).但由于單個(gè)龍卷的時(shí)空尺度較小,天氣雷達(dá)難以直接觀測(cè)到孕育龍卷發(fā)生的小尺度渦旋特征和更大的渦旋環(huán)境場,如中氣旋的結(jié)構(gòu)特征.與中國的每次單個(gè)龍卷有所不同,美國龍卷多數(shù)成群爆發(fā),即會(huì)有多個(gè)龍卷風(fēng)形成在一個(gè)天氣尺度或近地面次天氣(百千米)尺度的低壓系統(tǒng)中(Knupp et al., 2014).在美國南部,隨著低壓系統(tǒng)的長距離移動(dòng),其中可能同時(shí)或先后發(fā)生多個(gè)龍卷,即龍卷群發(fā)(Corfidi et al., 2010;Mercer et al., 2012).

    龍卷屬于小概率極端天氣.極端天氣是大氣中擾動(dòng)能量緩慢積累到一個(gè)臨界值時(shí)的快速釋放.觀測(cè)的大氣變量中包含隨太陽高度角與下墊面達(dá)到熱力-動(dòng)力平衡時(shí)的瞬變氣候態(tài)和瞬時(shí)擾動(dòng)分量,兩個(gè)分量都具有真實(shí)的物理含義.所以,極端天氣的發(fā)生與當(dāng)前大氣中的瞬時(shí)擾動(dòng)能量的積累有關(guān),表現(xiàn)以溫度擾動(dòng)為代表的內(nèi)能擾動(dòng)、以風(fēng)擾動(dòng)為代表的動(dòng)能擾動(dòng)等.怎樣從觀測(cè)的和數(shù)值模式預(yù)報(bào)產(chǎn)品中分離掉瞬變氣候參考態(tài),得到瞬時(shí)擾動(dòng)量分布的異常特征呢?2012年以來使用的大氣變量擾動(dòng)分解法提供了這種可能(丑紀(jì)范,2012).利用模式輸出產(chǎn)品,分解法已經(jīng)在國內(nèi)外的多種極端天氣的分析和預(yù)報(bào)中得到了應(yīng)用,如對(duì)北京暴雨(Jiang et al., 2016)、中國東部熱浪(Chen et al., 2017)、東北夏季低溫(Qian and Jiang,2014)和南方雨雪冰凍(錢維宏和張宗婕,2012)等方面的分析研究.因此,對(duì)龍卷風(fēng)這樣的小尺度罕見極端天氣,我們希望在它的發(fā)生之前有可捕捉的中尺度或大尺度擾動(dòng)環(huán)流和擾動(dòng)變量分布的環(huán)境形勢(shì).遼寧開原當(dāng)?shù)赜袣v史觀測(cè)記載的強(qiáng)龍卷例子并不多.對(duì)2019年7月3日下午發(fā)生在遼寧開原的龍卷例子,我們?cè)噲D從實(shí)時(shí)分析的大氣變量和歐洲中期天氣預(yù)報(bào)中心(ECMWF)模式產(chǎn)品中分離出擾動(dòng)分量,并結(jié)合地面觀測(cè)、雷達(dá)圖像和衛(wèi)星云圖演變的特征,尋找有利于龍卷發(fā)生的擾動(dòng)天氣環(huán)境.

    1 龍卷特征與資料和方法

    根據(jù)遼寧省氣象局在2019年7月4日的新聞通報(bào),受東北冷渦及蒙古氣旋前側(cè)低壓帶的共同影響,7月3日下午鐵嶺地區(qū)出現(xiàn)了龍卷、冰雹、短時(shí)雷暴降水等強(qiáng)對(duì)流天氣.天氣雷達(dá)顯示,強(qiáng)回波從昌圖縣進(jìn)入開原市西北部的金鉤子鎮(zhèn),向東南方向移動(dòng)到開原市區(qū)附近.開原市國家氣象觀測(cè)站在7月3日17時(shí)25分至37分觀測(cè)到開原市區(qū)出現(xiàn)龍卷,氣象站測(cè)得的最大瞬時(shí)風(fēng)速為23 m·s-1,同時(shí)出現(xiàn)了1小時(shí)25 mm的強(qiáng)降水,并且在開原市區(qū)出現(xiàn)了冰雹和雷電.在沈陽雷達(dá)站的7月3日17時(shí)22分0.5°仰角回波圖像上,開原市北部出現(xiàn)鉤狀回波并持續(xù)到17時(shí)44分,回波上呈現(xiàn)出中氣旋結(jié)構(gòu).遼寧省氣象局對(duì)事后的天氣影響評(píng)估中更具體地指出,龍卷始發(fā)于開原市金鉤子鎮(zhèn)金英村北約0.5 km,消失于中固鎮(zhèn)清水溝子村北約1.5 km,全程直線距離約15 km,強(qiáng)度為EF4級(jí).龍卷、鉤狀回波和中氣旋式樣的回波,它們?nèi)叩某叨炔煌?,但出現(xiàn)的時(shí)間和地點(diǎn)一致.此次強(qiáng)對(duì)流天氣過程導(dǎo)致了人員傷亡和房屋倒塌.

    本文使用了多套資料.其一是美國國家環(huán)境預(yù)報(bào)中心(NCEP)提供的全球預(yù)報(bào)系統(tǒng)(GFS)的實(shí)時(shí)分析資料(https:∥www.emc.ncep.noaa.gov/GFS/.php),時(shí)間間隔是3 h,水平分辨率為1.0°×1.0°經(jīng)緯格點(diǎn),垂直分31個(gè)氣壓層,對(duì)風(fēng)場資料內(nèi)插到了0.5 h分辨率上.其二是ECMWF模式產(chǎn)品提供的每天0000和1200世界時(shí)(UTC)分析場和對(duì)未來10天時(shí)間間隔為6 h(0000、0600、1200和1800UTC)的預(yù)報(bào)產(chǎn)品,水平分辨率為0.75°×0.75°經(jīng)緯格點(diǎn),垂直分8層(http:∥tigge.ecmwf.int/).其三是中國國家氣象中心天氣預(yù)報(bào)技術(shù)研發(fā)室對(duì)沈陽氣象雷達(dá)觀測(cè)的圖像分析.其四是日本葵花8氣象衛(wèi)星紅外云圖,時(shí)間間隔是10 min(http:∥www.jma.go.jp/en/gms/).

    擾動(dòng)分解方法已經(jīng)在多種類極端天氣的分析中有了介紹(錢維宏,2012a,b;錢維宏和于甜甜,2013;錢維宏等,2016,2017).這一方法的思想來源于中華祖先對(duì)氣候年循環(huán)中24節(jié)氣和72候的時(shí)間段劃分.年循環(huán)的氣候等分成24節(jié)氣和72候,可以稱為古代氣候鐘,它是古代人對(duì)中原地區(qū)地面氣候年循環(huán)的描述.依此,現(xiàn)代人可以利用當(dāng)今全球高時(shí)空分辨率多年觀測(cè)資料的積累,建立一個(gè)現(xiàn)代氣候鐘.它可以把年循環(huán)的氣候態(tài)細(xì)分成365天/年×24個(gè)時(shí)次/天=8760個(gè)時(shí)次的瞬變氣候態(tài),它包含年循環(huán)過程中的逐日變化和日循環(huán)過程中的逐時(shí)變化.用實(shí)況大氣變量或模式產(chǎn)品變量減去瞬變氣候變量就得到瞬時(shí)擾動(dòng).瞬變氣候和24節(jié)氣氣候都是不需要預(yù)報(bào)的,要預(yù)報(bào)的是相對(duì)瞬變氣候的瞬時(shí)擾動(dòng).雷達(dá)回波和云圖上的對(duì)流實(shí)體可以看作為極端天氣或擾動(dòng)天氣的實(shí)況.

    雷達(dá)回波中的反射率因子反映了大氣中的云和降水中的各種水汽凝結(jié)物對(duì)電磁波的后向散射強(qiáng)度,其回波反映了氣象目標(biāo)內(nèi)部降水粒子的尺度和密度分布.圖1a是國家氣象中心天氣預(yù)報(bào)技術(shù)研發(fā)室對(duì)2019年7月3日17時(shí)44分32秒沈陽雷達(dá)0.5°仰角觀測(cè)的回波強(qiáng)度分析.雷達(dá)回波呈東西帶狀,最強(qiáng)的回波在最西部分的尖緣.圖1a中的回波總體分布與850 hPa上的擾動(dòng)風(fēng)切變(黃箭頭)均為可分辨的中尺度系統(tǒng).為了細(xì)化回波的特征,圖1b是對(duì)雷達(dá)回波西端部分的放大,鉤狀回波特征明顯.鉤狀回波也對(duì)應(yīng)著一個(gè)擾動(dòng)風(fēng)下的氣旋式擾動(dòng)環(huán)流.所以,鉤狀回波(圖1b)和帶狀回波(圖1a)的形成有著共同的擾動(dòng)大背景.雖然氣象觀測(cè)站在17時(shí)25分至37分觀測(cè)到開原市區(qū)出現(xiàn)龍卷,但在其后的8 min內(nèi)龍卷的鉤狀回波并沒有消失.蘇北阜寧龍卷也有類似的鉤狀回波(錢維宏等,2017),但開原龍卷在17時(shí)44分32秒沿最強(qiáng)位置線(圖1b中白直線)的回波高度和強(qiáng)度都比阜寧的高和強(qiáng)(回波垂直剖面圖略).圖1b中的鉤狀回波結(jié)構(gòu)和圖1a中的總體回波分布與常規(guī)天氣圖和擾動(dòng)天氣圖上天氣系統(tǒng)之間的聯(lián)系,是需要進(jìn)一步考察的內(nèi)容.

    圖1 國家氣象中心天氣預(yù)報(bào)技術(shù)研發(fā)室對(duì)2019年7月3日17時(shí)44分32秒沈陽雷達(dá)0.5°仰角觀測(cè)回波的分析:(a)總體回波和(b)局部放大的回波(陰影,10 DBZ間隔). (a)中白色圓圈指示所在的龍卷風(fēng)暴,而(b)中的框區(qū)為龍卷風(fēng)暴主體部分. 黃色箭頭示意17∶ 30 BT 的850 hPa擾動(dòng)風(fēng)矢量Fig.1 0.5° elevated radar echo over Shenyang at 17∶ 44∶ 32 BT 3 July 2019 from (a) regional and (b) local scopes made by the Research and Development Office of Weather Forecast Technology, National Meteorological Center of China. The white solid circle in (a) denotes the location of the tornado and the white rectangle in (b) is the main body of the tornado. The yellow arrows denote the anomalous winds at 17∶ 30 BT at 850 hPa

    2 常規(guī)天氣圖上的形勢(shì)分析

    來自實(shí)時(shí)分析資料和模式產(chǎn)品繪制的常規(guī)天氣圖仍然是國內(nèi)外氣象業(yè)務(wù)部門預(yù)報(bào)和預(yù)警重要天氣的首選工具.在每天兩個(gè)時(shí)次(北京08和20時(shí))的常規(guī)天氣圖上,我們尋找當(dāng)日龍卷發(fā)生前08時(shí)500和850 hPa以及存在于高空天氣圖和地面形勢(shì)場中可能的強(qiáng)對(duì)流天氣潛在信息.圖2a中,500 hPa層的東北冷渦渦旋中心(C1和L1)位于黑龍江省的中部,冷的溫度中心處于低壓中心的西側(cè),而蒙古國在高壓系統(tǒng)(H)的控制下,開原(三角)的對(duì)流層中部處于位勢(shì)高度場的西高東低之間的強(qiáng)西北氣流中.這種溫壓場配置很容易在渦旋中心的西側(cè)形成橫槽,并有冷空氣的堆積.在2019年7月3日08時(shí)的500 hPa常規(guī)天氣圖上,三條階梯槽(圖2a中的粗虛線)由北向南分別出現(xiàn)在東北、渤海到江淮地區(qū),有利于冷空氣南下.

    低層850 hPa層上(圖2b),開原所在緯度帶的天氣形勢(shì)是,蒙古國的西部有一個(gè)氣旋(L2),蒙古國的東部受高壓脊(H)的影響,而我國東北有一個(gè)冷渦(C1和L1),開原的西側(cè)有一條低壓槽(低壓帶)或西北風(fēng)與偏西風(fēng)的切變線(粗虛線).這就是氣象部門在新聞通報(bào)中提到的天氣形勢(shì):開原的龍卷發(fā)生受到了東北冷渦及蒙古氣旋前側(cè)低壓帶的共同影響.東北三省的中東部地區(qū)低層大氣都在西南暖平流的影響下.中低層500 hPa和850 hPa之間的溫壓場與風(fēng)場的配置關(guān)系是有利于大氣層結(jié)不穩(wěn)定度加強(qiáng)的.類似的中低層大氣溫壓場特征在以往的東北冷渦背景下的強(qiáng)對(duì)流天氣形勢(shì)中也有反映(王寧等,2014;王婷婷等,2017).類似的水平階梯槽和垂直不穩(wěn)定形勢(shì)都會(huì)反映在圖2所示的廣大范圍的常規(guī)天氣圖上.但是,預(yù)報(bào)員難以確定圖1中與帶狀回波和鉤狀回波有緊密聯(lián)系的擾動(dòng)系統(tǒng)的具體位置和時(shí)間.

    圖2 2019年7月3日08時(shí)ECMWF模式產(chǎn)品中的實(shí)時(shí)分析場(a)500 hPa和(b)850 hPa位勢(shì)高度(黑色等值線,40gpm間隔)、溫度(灰色等值線,4℃間隔)和風(fēng)羽(單位:m·s-1).粗虛線為高度槽或風(fēng)的切變線,字母W和C表示暖和冷的位置,H和L指示高壓脊和低渦. 黑色三角為龍卷發(fā)生地開原Fig.2 Total geopotential height (black contours, 40 gpm interval), total temperature (grey contours, 4 ℃ interval) and wind barbs (unit: m·s-1) at 08∶ 00 BT 3 July 2019 based on the ECMWF analysis products at (a) 500 hPa and (b) 850 hPa. Thick dashed line denotes the height trough or wind shear. The letters W/C and H/L denote the centers of warm/cold masses and the locations of height ridge and low vortex. The black triangle denotes the location of tornado

    上述常規(guī)天氣圖形勢(shì)的分析表明,該時(shí)間段內(nèi)東北地區(qū)是有利于強(qiáng)對(duì)流天氣孕育的,但仍然難以推斷極端天氣,如強(qiáng)對(duì)流風(fēng)暴可能發(fā)生的具體位置和時(shí)間.所以,即使這樣的天氣形勢(shì)再次出現(xiàn)在常規(guī)天氣圖上,預(yù)報(bào)員也難以提前做出龍卷可能在哪個(gè)省的哪個(gè)地方發(fā)生的預(yù)報(bào)和預(yù)警.由于極端天氣的發(fā)生反映了大氣中擾動(dòng)能量積累到一定程度時(shí)的快速釋放,我們希望通過擾動(dòng)天氣圖的形勢(shì)分析,尋找可能發(fā)生強(qiáng)風(fēng)暴的局地環(huán)境潛勢(shì).

    3 大氣擾動(dòng)變量垂直分布結(jié)構(gòu)

    為了節(jié)省篇幅,我們?cè)谝韵碌姆治鲋胁辉倮L制大氣實(shí)時(shí)變量和模式產(chǎn)品的總量(原始場)垂直剖面圖,而是直接給出垂直的和水平的擾動(dòng)變量分析圖,在時(shí)間描述上使用北京時(shí).在現(xiàn)代條件下,如能對(duì)龍卷發(fā)生地有提前幾十個(gè)小時(shí)的預(yù)示征兆和臨近一個(gè)小時(shí)的預(yù)警信號(hào),對(duì)防災(zāi)減災(zāi)是很有意義的.從ECMWF模式預(yù)報(bào)時(shí)刻的實(shí)時(shí)分析到每6小時(shí)間隔的模式預(yù)報(bào)產(chǎn)品,我們逐一做了擾動(dòng)變量的分解分析.該模式提前6、12、18、24、30、36 h都預(yù)報(bào)出了可與實(shí)況比較一致的擾動(dòng)形勢(shì).于是,我們繼續(xù)分析ECMWF模式提前42 h和提前48 h預(yù)報(bào)的擾動(dòng)天氣形勢(shì).圖3是經(jīng)過龍卷發(fā)生地附近,由ECMWF模式提前42 h和提前48 h預(yù)報(bào)的高度擾動(dòng)、溫度擾動(dòng)、比濕擾動(dòng)和西風(fēng)擾動(dòng),繪制的東西方向和南北方向的垂直剖面分析.從垂直剖面圖3a和3b上看,東北地區(qū)上空的冷渦非常深厚,從700 hPa向上伸展到了整個(gè)對(duì)流層,溫度距平的負(fù)值中心C2在250 hPa,只在對(duì)流層的低層存在弱的擾動(dòng)高壓H1和在近地面水平方向上有冷暖氣團(tuán)的對(duì)比(W1與C1).圖3中,龍卷發(fā)生在模式預(yù)報(bào)的近地面冷暖擾動(dòng)氣團(tuán)的交界處,東北側(cè)是冷氣團(tuán),西南側(cè)是暖氣團(tuán),從對(duì)流層延伸下來的擾動(dòng)低壓槽(粗棕虛線)正好分隔了近地面層的冷暖擾動(dòng)氣團(tuán).中國東部地區(qū)夏季的暴雨帶多位于高度擾動(dòng)槽線和擾動(dòng)風(fēng)的輻合線上,但在常規(guī)天氣圖上暴雨帶位于總高度槽和總風(fēng)輻合線的右側(cè),低空急流的左側(cè)(Qian et al., 2013;Ai and Qian, 2020).它們之間的比較表明,擾動(dòng)天氣圖上的擾動(dòng)變量分布能夠確定地指示暴雨發(fā)生的位置.龍卷發(fā)生位置通常在冷暖氣團(tuán)交界處的偏暖氣團(tuán)一側(cè),這是在分析美國多次龍卷和蘇北多次龍卷的擾動(dòng)天氣圖時(shí)得到的認(rèn)識(shí)(錢維宏等,2017;Qian et al., 2019).所以,圖3中龍卷發(fā)生前的溫壓場空間配置關(guān)系再次證實(shí)了之前中國龍卷和美國龍卷的分析結(jié)果.同樣,中國東北龍卷與蘇北龍卷的擾動(dòng)形勢(shì)也相同,但它們?cè)诮孛鎸拥睦渑瘹鈭F(tuán)對(duì)比都比美國的弱小,反映的是擾動(dòng)內(nèi)能不如美國的大,所以只發(fā)生了單個(gè)龍卷,而美國的擾動(dòng)氣團(tuán)可以同時(shí)出現(xiàn)多個(gè)龍卷(Qian et al., 2019).

    2019年7月3日20時(shí),龍卷對(duì)地面的影響已經(jīng)結(jié)束了兩個(gè)半小時(shí),但龍卷在空中的背景擾動(dòng)渦旋還存在.ECMWF模式預(yù)報(bào)的第48小時(shí)(3日20時(shí)),對(duì)流層中高層的擾動(dòng)冷渦結(jié)構(gòu)在龍卷出現(xiàn)前后的幾個(gè)小時(shí)中沒有大的變化,近地面冷暖氣團(tuán)的交界線和高層延伸下來的擾動(dòng)高度槽線正好就在龍卷發(fā)生地(圖3c和3d).可見,龍卷出現(xiàn)時(shí),高層高度擾動(dòng)延伸下來的槽線和近地面冷暖氣團(tuán)分界線剛好定位在開原附近.這一結(jié)果表明,當(dāng)前數(shù)值模式可以在短期內(nèi)提前幾十個(gè)小時(shí)預(yù)示有可能孕育強(qiáng)對(duì)流天氣發(fā)生的時(shí)間和位置.ECMWF模式能夠在短期內(nèi)預(yù)報(bào)出可與實(shí)況較為一致的天氣環(huán)流形勢(shì),但是預(yù)報(bào)員難以直接從常規(guī)天氣圖上看到可指示有利于龍卷孕育的環(huán)境形勢(shì).龍卷的發(fā)生是局地的,但它的孕育需要一個(gè)大尺度的擾動(dòng)環(huán)境背景.這個(gè)擾動(dòng)環(huán)流系統(tǒng)從較大尺度上提供了局地強(qiáng)對(duì)流天氣可能發(fā)生的時(shí)間和地點(diǎn)的預(yù)示,它也為雷達(dá)和衛(wèi)星高時(shí)空分辨率圖像的檢測(cè)與預(yù)警提供了信息.這種預(yù)示信息在常規(guī)天氣圖上和模式預(yù)報(bào)產(chǎn)品中的表現(xiàn)沒有在擾動(dòng)天氣圖中的清晰.

    在實(shí)際的對(duì)流天氣和強(qiáng)降水分析中,人們同時(shí)考慮由總場基本量計(jì)算出的多種物理指標(biāo)和環(huán)流形勢(shì):如動(dòng)力學(xué)的渦度、散度和上升速度;熱力學(xué)的多種垂直穩(wěn)定度指數(shù)(K指數(shù)和對(duì)流有效位能(CAPE)等);水汽濕度條件中的溫度露點(diǎn)差、水汽通量、假相當(dāng)位溫、強(qiáng)水汽輻合等要素;中低層大氣中的次天氣或中尺度天氣系統(tǒng)的天氣流型(槽、切變線、急流位置、逆溫層、冷暖鋒面、低層濕舌等);最后把這些進(jìn)行搭配組合,做出對(duì)流天氣的判斷(唐曉文等,2010;張小玲等,2012;劉國忠等,2013).考慮到這些物理量和流型不是相互獨(dú)立的,擾動(dòng)法直接使用基本量(溫、壓、濕、風(fēng))中的擾動(dòng)分量空間分布組合特征,它們能夠在業(yè)務(wù)中直觀、快捷地指出有利于對(duì)流天氣發(fā)生的位置和強(qiáng)度隨時(shí)間的變化.

    圖3e是ECMWF模式從7月1日20時(shí)起報(bào)的第42小時(shí),即3日14時(shí)的比濕擾動(dòng)在南北剖面上的分布.龍卷發(fā)生日(9個(gè)小時(shí)前)的當(dāng)天上午,近地面層只有弱的正比濕中心(W1).龍卷發(fā)生后的3日20時(shí),比濕擾動(dòng)的正中心接近地面(圖3f中W1).與以前的分析結(jié)果一致,東北龍卷發(fā)生的低層水汽含量和濕層厚度比江淮和華南龍卷的低(王秀明等,2015;錢維宏等,2017).低層大氣擾動(dòng)風(fēng)上,龍卷發(fā)生日和發(fā)生地的南側(cè)為南風(fēng)擾動(dòng),北側(cè)為北風(fēng)擾動(dòng).龍卷發(fā)生前,開原上空的對(duì)流層存在異常的西風(fēng)(圖3e中P),而在北側(cè)有更強(qiáng)的東風(fēng)擾動(dòng)(圖3e中N),對(duì)流層大氣中存在很強(qiáng)的水平風(fēng)速切變(圖3e中P與N之間的密集風(fēng)速等值線).在850 hPa上,環(huán)繞龍卷發(fā)生地的東西風(fēng)擾動(dòng)切變?nèi)匀磺宄?在龍卷消亡后,這樣的風(fēng)切變也減弱了(圖3f中P與N之間的稀疏風(fēng)速等值線).龍卷是大氣中發(fā)生的擾動(dòng)動(dòng)能高度集中的小尺度環(huán)流系統(tǒng).這種擾動(dòng)動(dòng)能的高度集中必須有一個(gè)更大的擾動(dòng)動(dòng)能環(huán)境,龍卷發(fā)生前的大氣中(圖3e)就具有比龍卷發(fā)生后的大氣中(圖3f)更有利的大尺度擾動(dòng)動(dòng)能環(huán)境.

    4 大氣擾動(dòng)變量水平分布

    擾動(dòng)天氣圖上的大氣擾動(dòng)系統(tǒng)具有可識(shí)別的空間三維結(jié)構(gòu)特征.我們進(jìn)一步考察龍卷發(fā)生前ECMWF模式提前24 h內(nèi)預(yù)報(bào)的擾動(dòng)形勢(shì)在水平擾動(dòng)天氣圖上的分布.圖4是ECMWF模式預(yù)報(bào)的850 hPa風(fēng)擾動(dòng)和比濕擾動(dòng).3日02時(shí)(圖4a),東北擾動(dòng)渦旋的低層、吉林省的西北部和渤海三地各有一個(gè)擾動(dòng)比濕中心.一條擾動(dòng)風(fēng)切變線(紅實(shí)線)從東北擾動(dòng)冷渦向西南延伸到開原的北側(cè).3日08時(shí)至14時(shí)(圖4b和圖4c),擾動(dòng)風(fēng)切變線與其北側(cè)的最強(qiáng)擾動(dòng)風(fēng)(紅箭頭)一起,配合比濕中心到達(dá)開原.3日20時(shí)(圖4d),擾動(dòng)風(fēng)切變線東移到了吉林省和遼寧省的東部邊境附近.7月3日的四個(gè)時(shí)次,擾動(dòng)風(fēng)切變線的尾端配合最強(qiáng)的擾動(dòng)氣流(紅箭頭)由西北向東南掃過了開原地區(qū).切變線尾端的強(qiáng)偏北擾動(dòng)氣流與切變線南側(cè)的偏西氣流之間具有較大的正渦度擾動(dòng),它是有利于強(qiáng)對(duì)流天氣發(fā)生的潛勢(shì)區(qū)域.

    ECMWF模式2019年7月1日20時(shí)起報(bào)的3日14時(shí)和20時(shí)的風(fēng)擾動(dòng)和比濕擾動(dòng)也很好地指示了上述850 hPa擾動(dòng)風(fēng)切變線隨時(shí)間掃過了開原地區(qū).3日14時(shí)(圖4e),ECMWF模式提前42 h預(yù)報(bào)出了擾動(dòng)風(fēng)切變線和一個(gè)弱的濕擾動(dòng)中心位于開原地區(qū).3日20時(shí)(圖4f),ECMWF模式提前48 h預(yù)報(bào)出了擾動(dòng)風(fēng)切變線向東移動(dòng)到了吉林省和遼寧省的東部邊緣.從單一擾動(dòng)變量看,如擾動(dòng)風(fēng)切變線掃過的范圍遠(yuǎn)比發(fā)生龍卷的地方大,所以單一擾動(dòng)變量難以指示龍卷的發(fā)生地.龍卷的發(fā)生不但需要有氣候環(huán)境條件,還需要多擾動(dòng)變量的綜合環(huán)境條件.擾動(dòng)變量分布上,龍卷發(fā)生在最強(qiáng)擾動(dòng)北風(fēng)與擾動(dòng)西風(fēng)具有最大擾動(dòng)渦度的地方,又配合在一個(gè)正的擾動(dòng)比濕條件的地方.這些中尺度的擾動(dòng)環(huán)境為小尺度龍卷的可能發(fā)生準(zhǔn)備了條件.在很多的類似環(huán)境條件下,發(fā)生的可能不是龍卷,而是冰雹或雷暴等極端天氣.

    圖4 ECMWF模式2019年7月2日20時(shí)起報(bào)的3日(a)02時(shí),(b)08時(shí),(c)14時(shí)和(d)20時(shí)的850 hPa風(fēng)擾動(dòng)(單位:m·s-1,紅線為切變線,紅色箭頭指示偏北大風(fēng))和比濕擾動(dòng)(陰影,單位:10-3 kg·kg-1,W為濕中心); (e)和(f)為ECMWF模式2019年7月1日20時(shí)起報(bào)的3日14時(shí)和20時(shí)的風(fēng)擾動(dòng)和比濕擾動(dòng)Fig.4 Anomalous wind vector (units: m·s-1) and anomalous specific humidity (shading, units: 10-3 kg·kg-1) at 850 hPa based on the ECMWF model prediction at (a) 02∶ 00 BT 3, (b) 08∶ 00 BT 3, (c) 14∶ 00 BT 3, (d) 20∶ 00 BT on 3 July 2019 initiated from 20∶ 00 BT on 2 July 2019; (e) and (f) same as (a) except at 14∶ 00 BT on 3 and 20∶ 00 BT on 3 July 2019 initiated from 20∶ 00 BT on 1 July 2019. The red lines are the shear lines of winds at 850 hPa, and the red arrows indicate the strong northerlies. The letter W denotes the wet center of specific humidity anomalies

    5 云系演變與擾動(dòng)氣流分布

    上述分析的基本大氣擾動(dòng)變量,包括擾動(dòng)風(fēng)分布和擾動(dòng)比濕只是反映有可能出現(xiàn)強(qiáng)對(duì)流天氣的環(huán)境條件,尚難以判斷強(qiáng)對(duì)流的性質(zhì),如龍卷與一般性的對(duì)流,也難以確定強(qiáng)對(duì)流的更具體位置.在可預(yù)見的強(qiáng)對(duì)流天氣環(huán)境條件下,我們?cè)噲D從雷達(dá)回波和衛(wèi)星云圖上加密監(jiān)測(cè)可能出現(xiàn)的龍卷或其他強(qiáng)對(duì)流天氣的發(fā)生位置.作為在預(yù)報(bào)實(shí)踐中的可能應(yīng)用,圖5給出了2019年7月3日11時(shí)至21時(shí)日本葵花8氣象衛(wèi)星紅外云圖演變與美國GFS實(shí)時(shí)分析的850 hPa擾動(dòng)氣流的對(duì)應(yīng)關(guān)系.國家氣象站觀測(cè)到的龍卷出現(xiàn)時(shí)間在17時(shí)25分到37分,天氣雷達(dá)觀測(cè)到的強(qiáng)鉤狀回波在17時(shí)44分32秒.圖5a是11∶ 00時(shí)的云圖與擾動(dòng)氣流,最強(qiáng)的云帶沿?cái)_動(dòng)氣流輻合線,呈現(xiàn)東西向分布,出現(xiàn)在開原的東北方.12∶ 00時(shí)(圖略),云帶增強(qiáng),云帶西端尖緣部分清晰,云帶北側(cè)為東北擾動(dòng)氣流,南側(cè)為偏西擾動(dòng)氣流.在這樣的氣旋式擾動(dòng)氣流下,13∶ 00時(shí)云帶的西端尖緣部分向西南方向發(fā)展(圖5b).從15∶ 00時(shí)(圖5c),云帶進(jìn)一步增強(qiáng)向西南方向延伸,逐漸向開原靠近,云帶北側(cè)為較強(qiáng)的西北擾動(dòng)氣流,而南側(cè)為較弱的西南擾動(dòng)氣流.17∶ 00時(shí)是氣象站觀測(cè)到龍卷前的25 min,云帶的西端更為強(qiáng)盛,西端尖緣部分到達(dá)開原,西端云帶北側(cè)的強(qiáng)偏北擾動(dòng)氣流與云帶南側(cè)的弱偏西擾動(dòng)氣流之間形成氣旋式的擾動(dòng)渦度(圖5d).17∶ 30時(shí)是氣象站觀測(cè)到龍卷和雷達(dá)上有鉤狀回波的時(shí)刻,最強(qiáng)的云帶西端尖緣部分與擾動(dòng)氣流輻合線一起在開原.從18∶ 00時(shí)(圖5e)至19∶ 00時(shí),隨著擾動(dòng)氣流輻合線北側(cè)擾動(dòng)氣流的減弱,云帶變寬.云帶在20∶ 00時(shí)(圖5f)至21∶ 00時(shí)變得更松散,云帶西端位置也向南移動(dòng),離開開原.

    在2019年7月3日11∶ 00時(shí)至21∶ 00時(shí)的10個(gè)小時(shí)云圖強(qiáng)度和位置演變與850 hPa擾動(dòng)風(fēng)場的關(guān)系中,我們能夠快速地判斷它們之間的密切聯(lián)系.

    圖5 2019年7月3日(a) 11∶ 00時(shí),(b) 13∶ 00時(shí),(c) 15∶ 00時(shí),(d) 17∶ 00時(shí),(e) 18∶ 00時(shí)和(f) 20∶ 00時(shí)的葵花8云圖與GFS分析產(chǎn)品的850 hPa擾動(dòng)風(fēng)隨時(shí)間的變化. 紅三角指示開原位置. 紅線為偏北風(fēng)與偏西風(fēng)之間的擾動(dòng)氣流輻合線Fig.5 Himawari-8 image and anomalous wind vector at 850 hPa from GFS analysis product at (a) 11∶ 00 BT, (b) 13∶ 00 BT, (c) 15∶ 00 BT, (d) 17∶ 00 BT, (e) 18∶ 00 BT, and (f) 20∶ 00 BT on 3 July 2019. The red triangles denote the location of Kaiyuan and the red lines are convergence line of anomalous winds between northerlies and westerlies

    在雷達(dá)回波上,中氣旋的形成(圖1a)和強(qiáng)度變化也反映與這條云帶的演變對(duì)應(yīng).云帶亮度和位置變化也能夠反映對(duì)流的強(qiáng)度和移動(dòng),而真正背后的動(dòng)力學(xué)機(jī)制是擾動(dòng)氣流的輻合(圖1中的黃箭頭和圖5中的藍(lán)箭頭),即直接相關(guān)于擾動(dòng)氣流渦度(圖略).這樣的變化關(guān)系能夠反映出對(duì)流天氣強(qiáng)度和位置隨時(shí)間變化的動(dòng)力學(xué)原因,但困難仍然是,難以確認(rèn)是否會(huì)出現(xiàn)龍卷,還是一般的雷暴天氣.對(duì)照?qǐng)D1a與圖5e和圖5f中各時(shí)刻的東西向雷達(dá)回波與云帶西端的位置對(duì)應(yīng)關(guān)系,它們都形成在較強(qiáng)的北側(cè)偏北擾動(dòng)氣流與較弱的南側(cè)偏西擾動(dòng)氣流之間的正擾動(dòng)渦度區(qū)中.在龍卷發(fā)生的前后近10個(gè)小時(shí)中,最亮云帶和最強(qiáng)雷達(dá)回波的位置都能用850 hPa上的擾動(dòng)氣流輻合線和動(dòng)力學(xué)指標(biāo)量(擾動(dòng)渦度)來描述.在理論和應(yīng)用上,強(qiáng)對(duì)流性極端天氣(冰雹、龍卷和雷暴)的發(fā)生是潮濕氣團(tuán)受到強(qiáng)烈抬升的結(jié)果.它們是正的擾動(dòng)比濕與正的擾動(dòng)渦度重疊的地方,可用一個(gè)濕渦度擾動(dòng)(MVA)的指標(biāo)量表征(Qian et al., 2015),它也在本文的圖3至圖4中得到了再次驗(yàn)證.

    6 結(jié)論與討論

    高影響的極端天氣預(yù)報(bào)是世界性難題,其中龍卷預(yù)報(bào)更是難題中的難題,因?yàn)樗某叨刃?,但引發(fā)的災(zāi)害極大.本文對(duì)遼寧省開原地區(qū)發(fā)生的這例龍卷強(qiáng)對(duì)流過程的分析,從有效預(yù)報(bào)時(shí)間和擾動(dòng)天氣形勢(shì)分析的角度關(guān)注來自數(shù)值模式預(yù)報(bào)產(chǎn)品中的擾動(dòng)變量空間結(jié)構(gòu),從天氣過程演變和局地天氣特點(diǎn)的角度關(guān)注衛(wèi)星云圖和雷達(dá)回波隨實(shí)況擾動(dòng)氣流的變化趨勢(shì),可為相關(guān)部門提供一條預(yù)報(bào)強(qiáng)對(duì)流天氣發(fā)生的參考思路.

    (1)綜合比較當(dāng)前常規(guī)與擾動(dòng)天氣分析結(jié)果.我們首先可以從常規(guī)天氣圖上,或中-短期數(shù)值模式產(chǎn)品中,判斷一般性的對(duì)流不穩(wěn)定形勢(shì).從具體的做法出發(fā),預(yù)報(bào)員可以從圖2中分析500 hPa東北冷渦中心西側(cè)的橫槽與冷氣團(tuán)中心的配置關(guān)系.一旦橫槽帶著冷空氣轉(zhuǎn)豎,同時(shí)低層850 hPa有暖氣團(tuán)的存在或發(fā)展,即表明大氣層結(jié)不穩(wěn)定度在增加.如果又采用擾動(dòng)天氣分析方法,那么擾動(dòng)風(fēng)的強(qiáng)中心、擾動(dòng)風(fēng)切變線(或輻合線)、擾動(dòng)位勢(shì)高度槽、擾動(dòng)比濕的正中心和冷暖擾動(dòng)氣團(tuán)的邊界位置,它們的綜合分析會(huì)大大提高對(duì)強(qiáng)對(duì)流天氣發(fā)生時(shí)間和地點(diǎn)的判斷.

    (2)基于模式產(chǎn)品的擾動(dòng)變量分析.根據(jù)中國氣象局下發(fā)給各省市的ECMWF模式預(yù)報(bào)產(chǎn)品,地市縣預(yù)報(bào)員可以查看模式預(yù)報(bào)產(chǎn)品的當(dāng)前時(shí)刻和未來短期時(shí)刻的擾動(dòng)變量垂直分布配置(圖3)和水平分布特征(圖4),確定本地未來短期內(nèi)是否會(huì)受到多變量擾動(dòng)系統(tǒng)的影響.當(dāng)前ECMWF提供的未來10天模式預(yù)報(bào)產(chǎn)品中,對(duì)未來6天內(nèi)夏季中國的擾動(dòng)系統(tǒng)對(duì)降水預(yù)報(bào)的指示能力,要大大好于降水產(chǎn)品的直接預(yù)報(bào)(Qian et al., 2013).隨著將來極端天氣預(yù)報(bào)中人工智能技術(shù)的運(yùn)用,建立如本文類似的歷史擾動(dòng)變量空間分布結(jié)構(gòu)與地面極端天氣事件庫,有望在未來對(duì)不同種類的極端天氣環(huán)境形勢(shì)的識(shí)別和預(yù)報(bào)有所幫助.

    (3)衛(wèi)星云圖與擾動(dòng)風(fēng)定位極端天氣.衛(wèi)星云圖和雷達(dá)回波上看到的實(shí)體圖像,物理上反映的是大氣中云和降水中的各種水汽凝結(jié)物的存在.但它們也反映了大氣中擾動(dòng)能量緩慢積累到一個(gè)臨界值時(shí)的爆發(fā)或釋放.雷達(dá)回波上看到的實(shí)體已經(jīng)是強(qiáng)天氣爆發(fā)的實(shí)況了(圖1).云圖與擾動(dòng)風(fēng)的綜合分析可以得到一個(gè)隨時(shí)間演變的過程,也就是擾動(dòng)能量集中與釋放的過程(圖5).云圖中,強(qiáng)對(duì)流處云頂溫度低和邊緣清晰,表明強(qiáng)對(duì)流天氣發(fā)展強(qiáng)烈,是龍卷、冰雹和雷雨大風(fēng)可以定位的地方.實(shí)時(shí)監(jiān)測(cè)云圖和回波與擾動(dòng)風(fēng)場的變化過程,如擾動(dòng)氣流輻合線的變化或濕渦度擾動(dòng)的位置,是臨近龍卷等強(qiáng)對(duì)流天氣預(yù)警的最佳線索.當(dāng)然,用高時(shí)空分辨率的區(qū)域數(shù)值模式預(yù)報(bào)短期的擾動(dòng)風(fēng)分布和計(jì)算擾動(dòng)風(fēng)物理量(如濕渦度擾動(dòng)),對(duì)預(yù)報(bào)強(qiáng)對(duì)流天氣的發(fā)生時(shí)間和位置是有意義的.

    (4)擾動(dòng)天氣圖揭示的極端天氣本質(zhì).極端天氣,特別是龍卷的發(fā)生,是大氣中擾動(dòng)能量緩慢積累到一個(gè)臨界值時(shí)的快速釋放.擾動(dòng)能量就體現(xiàn)在溫度擾動(dòng)、位勢(shì)高度擾動(dòng)、風(fēng)擾動(dòng)和比濕擾動(dòng)的分布中.這些基本擾動(dòng)變量對(duì)應(yīng)龍卷發(fā)生時(shí)的空間分布特征能夠很好地反映在擾動(dòng)天氣圖上,也存在于常規(guī)天氣圖中,但難以顯現(xiàn)在常規(guī)天氣圖上.如果模式能在短期內(nèi)正確預(yù)報(bào)出總場變量,那么就可以分解并繪制出擾動(dòng)天氣圖.預(yù)報(bào)員就可以迅速地從擾動(dòng)天氣圖上識(shí)別和預(yù)報(bào)極端天氣的種類和強(qiáng)度.中國和美國的龍卷和暴雨例子中(錢維宏等,2017;Qian et al., 2019),強(qiáng)對(duì)流天氣也都發(fā)生在高層高度擾動(dòng)延伸下來的槽線和近地面冷暖氣團(tuán)分界線的位置,表明擾動(dòng)變量分布具有動(dòng)力學(xué)意義.溫度擾動(dòng)和高度擾動(dòng)之間滿足靜力平衡關(guān)系,即依據(jù)靜力平衡關(guān)系可以用高度擾動(dòng)計(jì)算出溫度擾動(dòng),也可以依據(jù)地轉(zhuǎn)關(guān)系用高度擾動(dòng)計(jì)算出風(fēng)擾動(dòng).溫度擾動(dòng)中的冷暖氣團(tuán)強(qiáng)度對(duì)比就是大氣中的不穩(wěn)定能量,它釋放時(shí)產(chǎn)生的對(duì)流運(yùn)動(dòng)與正的水汽異常配合就會(huì)發(fā)生龍卷、冰雹和強(qiáng)降水等極端天氣.

    致謝國家氣象中心天氣預(yù)報(bào)技術(shù)研發(fā)室朱文劍給予的雷達(dá)圖像支持,遼寧省氣象臺(tái)和科技處給予的開原龍卷相關(guān)實(shí)況信息.

    猜你喜歡
    龍卷強(qiáng)對(duì)流氣流
    氣流的威力
    哈爾濱2020年一次局地強(qiáng)對(duì)流天氣分析
    黑龍江氣象(2021年2期)2021-11-05 07:06:54
    海南龍卷統(tǒng)計(jì)特征
    突發(fā)性強(qiáng)對(duì)流天氣快速識(shí)別預(yù)警改進(jìn)方法
    我國龍卷預(yù)警試驗(yàn)的前沿陣地
    ——訪佛山市龍卷風(fēng)研究中心炎利軍主任
    龍卷的雷達(dá)探測(cè)研究進(jìn)展
    青藏高原東北部地區(qū)一次強(qiáng)對(duì)流天氣特征分析
    龍卷過程中雷達(dá)合成切變產(chǎn)品的分析
    固體運(yùn)載火箭變軌發(fā)動(dòng)機(jī)噴管氣流分離研究
    飛片下的空氣形成的“超強(qiáng)高速氣流刀”
    免费看日本二区| 在线观看66精品国产| 久久精品人妻少妇| 丰满乱子伦码专区| 少妇人妻精品综合一区二区 | 99国产极品粉嫩在线观看| 国产男人的电影天堂91| 一级毛片我不卡| 夜夜夜夜夜久久久久| 国产激情偷乱视频一区二区| 1000部很黄的大片| av又黄又爽大尺度在线免费看 | 国产真实乱freesex| 久久久久久久久久久丰满| 亚洲欧美日韩东京热| 国内少妇人妻偷人精品xxx网站| 久久久久久久久久黄片| 99在线视频只有这里精品首页| 精品人妻一区二区三区麻豆| 欧美三级亚洲精品| 亚洲国产欧洲综合997久久,| 春色校园在线视频观看| 联通29元200g的流量卡| 中文字幕av成人在线电影| 国产精品免费一区二区三区在线| 99国产精品一区二区蜜桃av| 国内久久婷婷六月综合欲色啪| 成人高潮视频无遮挡免费网站| 国产亚洲5aaaaa淫片| 成人特级av手机在线观看| 精品一区二区三区视频在线| 麻豆一二三区av精品| 人妻夜夜爽99麻豆av| 亚洲18禁久久av| 免费观看的影片在线观看| 日韩欧美国产在线观看| 日本成人三级电影网站| 男女做爰动态图高潮gif福利片| 深夜精品福利| 国产精品,欧美在线| 尾随美女入室| 国产在线男女| 日韩欧美精品免费久久| 国产亚洲5aaaaa淫片| 久久欧美精品欧美久久欧美| h日本视频在线播放| 色哟哟哟哟哟哟| 2021天堂中文幕一二区在线观| 国产精品av视频在线免费观看| 欧美精品国产亚洲| 少妇裸体淫交视频免费看高清| 午夜老司机福利剧场| 三级毛片av免费| 亚洲美女搞黄在线观看| 国产激情偷乱视频一区二区| 亚洲成人av在线免费| 黄片无遮挡物在线观看| 男女做爰动态图高潮gif福利片| 日本成人三级电影网站| 熟女人妻精品中文字幕| 99久久精品国产国产毛片| 秋霞在线观看毛片| 少妇熟女欧美另类| 久久精品久久久久久久性| 久久精品国产鲁丝片午夜精品| 小说图片视频综合网站| 久久精品91蜜桃| 亚洲成人中文字幕在线播放| 麻豆国产97在线/欧美| 天美传媒精品一区二区| 亚洲国产精品成人综合色| 久久久久网色| 如何舔出高潮| 国产伦精品一区二区三区四那| 精品少妇黑人巨大在线播放 | 久久6这里有精品| 亚洲熟妇中文字幕五十中出| 精品一区二区免费观看| 日韩欧美一区二区三区在线观看| 变态另类丝袜制服| 青春草亚洲视频在线观看| 久久久精品欧美日韩精品| 一级毛片我不卡| 又爽又黄无遮挡网站| 亚洲人成网站高清观看| 国产av一区在线观看免费| 美女高潮的动态| 国产老妇伦熟女老妇高清| 天天躁日日操中文字幕| 色视频www国产| 午夜精品一区二区三区免费看| 看黄色毛片网站| 夜夜看夜夜爽夜夜摸| 亚洲国产精品sss在线观看| 亚洲国产欧美在线一区| 免费观看人在逋| 麻豆国产av国片精品| 久久草成人影院| 欧美极品一区二区三区四区| 亚洲天堂国产精品一区在线| 欧美极品一区二区三区四区| 桃色一区二区三区在线观看| 麻豆国产av国片精品| 能在线免费看毛片的网站| 99热这里只有是精品在线观看| 色尼玛亚洲综合影院| 国产探花极品一区二区| 国产精品久久久久久久电影| 九九在线视频观看精品| 欧美另类亚洲清纯唯美| 乱人视频在线观看| 国产精品综合久久久久久久免费| 亚洲天堂国产精品一区在线| 日日摸夜夜添夜夜爱| 2021天堂中文幕一二区在线观| 亚洲av成人精品一区久久| 99久久无色码亚洲精品果冻| 91aial.com中文字幕在线观看| 国产精品野战在线观看| 亚洲欧美清纯卡通| 国产高清激情床上av| 人妻制服诱惑在线中文字幕| 美女 人体艺术 gogo| 亚洲av免费在线观看| 国产精品一区二区在线观看99 | 日本爱情动作片www.在线观看| 在线观看av片永久免费下载| 久久久a久久爽久久v久久| 天天躁日日操中文字幕| 免费观看a级毛片全部| 国产成人一区二区在线| 久久精品夜色国产| 天堂√8在线中文| 亚洲中文字幕日韩| 国产精品久久电影中文字幕| 精品免费久久久久久久清纯| 亚洲精品日韩在线中文字幕 | 久久欧美精品欧美久久欧美| 一级毛片aaaaaa免费看小| 久久九九热精品免费| 日日撸夜夜添| 狠狠狠狠99中文字幕| 亚洲av中文av极速乱| 免费人成在线观看视频色| 99久久精品国产国产毛片| 丰满人妻一区二区三区视频av| 看十八女毛片水多多多| av天堂中文字幕网| 女的被弄到高潮叫床怎么办| 亚洲av第一区精品v没综合| 最好的美女福利视频网| 亚洲精华国产精华液的使用体验 | 亚洲欧美日韩东京热| 99久国产av精品| 一夜夜www| 成人三级黄色视频| 内射极品少妇av片p| 色视频www国产| 少妇裸体淫交视频免费看高清| 久久人人精品亚洲av| 伦精品一区二区三区| 欧美3d第一页| 欧美日韩乱码在线| 亚洲欧美精品专区久久| 久久久久免费精品人妻一区二区| 少妇人妻精品综合一区二区 | 一个人观看的视频www高清免费观看| 亚洲成人中文字幕在线播放| 成人高潮视频无遮挡免费网站| 欧美成人精品欧美一级黄| 99国产精品一区二区蜜桃av| 精品国产三级普通话版| 精品不卡国产一区二区三区| 又爽又黄a免费视频| 能在线免费看毛片的网站| 99热只有精品国产| 国产一级毛片七仙女欲春2| 精品久久久久久久久久免费视频| 久久九九热精品免费| 精品久久久久久久久av| 欧美又色又爽又黄视频| 亚洲,欧美,日韩| 欧美成人免费av一区二区三区| 免费看美女性在线毛片视频| 亚洲自拍偷在线| 看十八女毛片水多多多| 国产又黄又爽又无遮挡在线| 99在线视频只有这里精品首页| 色综合站精品国产| av天堂在线播放| 国产免费男女视频| 亚洲国产精品久久男人天堂| 久久久久免费精品人妻一区二区| 大型黄色视频在线免费观看| 男人舔女人下体高潮全视频| 久久精品国产清高在天天线| 免费不卡的大黄色大毛片视频在线观看 | 一本一本综合久久| 超碰av人人做人人爽久久| 男女做爰动态图高潮gif福利片| 国内精品美女久久久久久| 日韩在线高清观看一区二区三区| 国产单亲对白刺激| 亚洲va在线va天堂va国产| 成年女人永久免费观看视频| 哪个播放器可以免费观看大片| 亚洲av二区三区四区| 国产在线精品亚洲第一网站| 亚洲人与动物交配视频| 非洲黑人性xxxx精品又粗又长| 大又大粗又爽又黄少妇毛片口| 免费一级毛片在线播放高清视频| av在线亚洲专区| 网址你懂的国产日韩在线| 九九在线视频观看精品| 国产爱豆传媒在线观看| 蜜臀久久99精品久久宅男| a级毛片免费高清观看在线播放| 能在线免费观看的黄片| 国产高清有码在线观看视频| 亚洲国产欧洲综合997久久,| 国产 一区精品| 欧美精品一区二区大全| 大香蕉久久网| 男女无遮挡免费网站观看| 如何舔出高潮| 中文乱码字字幕精品一区二区三区| 国产成人av激情在线播放 | 成人影院久久| 99热网站在线观看| 久久国产精品男人的天堂亚洲 | 美女xxoo啪啪120秒动态图| 亚洲综合色惰| 久久99精品国语久久久| 国产熟女欧美一区二区| 久久亚洲国产成人精品v| 午夜日本视频在线| 看十八女毛片水多多多| 80岁老熟妇乱子伦牲交| 久久精品国产鲁丝片午夜精品| 午夜激情福利司机影院| 美女福利国产在线| 国产片内射在线| 99久久中文字幕三级久久日本| 成人手机av| 大话2 男鬼变身卡| 亚洲中文av在线| 赤兔流量卡办理| av黄色大香蕉| 欧美3d第一页| 欧美日韩av久久| 热re99久久国产66热| 成人18禁高潮啪啪吃奶动态图 | 青春草亚洲视频在线观看| 最近最新中文字幕免费大全7| 只有这里有精品99| 人妻夜夜爽99麻豆av| 大香蕉久久网| 久久亚洲国产成人精品v| 亚洲av成人精品一区久久| 中文字幕亚洲精品专区| 久久久久国产精品人妻一区二区| 亚洲综合精品二区| av在线播放精品| 国产亚洲最大av| 中国三级夫妇交换| 国产欧美日韩综合在线一区二区| 亚洲精品亚洲一区二区| 大陆偷拍与自拍| 久久鲁丝午夜福利片| 日韩欧美一区视频在线观看| 在线观看人妻少妇| 久久久久久人妻| 超碰97精品在线观看| 久久久久网色| 国产在视频线精品| 一二三四中文在线观看免费高清| 亚洲av电影在线观看一区二区三区| 亚洲av.av天堂| 美女主播在线视频| 蜜臀久久99精品久久宅男| 成人毛片60女人毛片免费| 大陆偷拍与自拍| 精品酒店卫生间| 人人妻人人爽人人添夜夜欢视频| av专区在线播放| 国产一区有黄有色的免费视频| 九九爱精品视频在线观看| 久久99热6这里只有精品| 狂野欧美激情性bbbbbb| 国产精品嫩草影院av在线观看| 亚洲婷婷狠狠爱综合网| 久久久久久久精品精品| 国产视频首页在线观看| 在线观看三级黄色| 国产精品国产三级国产av玫瑰| 亚洲欧美成人综合另类久久久| 欧美日韩视频精品一区| 国产黄频视频在线观看| 亚洲一区二区三区欧美精品| 中文天堂在线官网| 久久人妻熟女aⅴ| 在线天堂最新版资源| 亚洲国产欧美在线一区| 国产精品一区二区在线不卡| 一级a做视频免费观看| 中文精品一卡2卡3卡4更新| 精品国产一区二区久久| 日韩精品有码人妻一区| av福利片在线| 美女主播在线视频| 婷婷色综合大香蕉| 观看av在线不卡| 亚洲色图 男人天堂 中文字幕 | 三级国产精品片| 一边摸一边做爽爽视频免费| 纯流量卡能插随身wifi吗| 99精国产麻豆久久婷婷| 久久久久人妻精品一区果冻| 中文乱码字字幕精品一区二区三区| 精品亚洲乱码少妇综合久久| 十八禁网站网址无遮挡| 人妻系列 视频| 国产毛片在线视频| 母亲3免费完整高清在线观看 | 亚洲精品色激情综合| 桃花免费在线播放| 黑人高潮一二区| 国产成人精品婷婷| 欧美另类一区| 青青草视频在线视频观看| 啦啦啦中文免费视频观看日本| 亚洲精品av麻豆狂野| 国产免费一区二区三区四区乱码| 嫩草影院入口| 一本久久精品| a级毛片黄视频| 一级黄片播放器| 久久亚洲国产成人精品v| 视频在线观看一区二区三区| av又黄又爽大尺度在线免费看| 欧美日韩国产mv在线观看视频| 精品99又大又爽又粗少妇毛片| 欧美人与善性xxx| 国产精品国产三级专区第一集| 久久久精品94久久精品| 99热这里只有精品一区| 亚洲久久久国产精品| 免费高清在线观看日韩| 高清黄色对白视频在线免费看| 亚洲精品乱码久久久久久按摩| av播播在线观看一区| 国产伦精品一区二区三区视频9| av免费在线看不卡| 2018国产大陆天天弄谢| 狂野欧美激情性xxxx在线观看| 日韩成人av中文字幕在线观看| 美女内射精品一级片tv| h视频一区二区三区| 日韩伦理黄色片| 亚洲四区av| 欧美人与善性xxx| 久久国产精品大桥未久av| 天天操日日干夜夜撸| 国产午夜精品久久久久久一区二区三区| 婷婷色综合大香蕉| 日韩中文字幕视频在线看片| 搡老乐熟女国产| 热re99久久国产66热| 亚洲国产精品成人久久小说| 国产69精品久久久久777片| 精品99又大又爽又粗少妇毛片| 免费少妇av软件| 欧美激情 高清一区二区三区| 久久久久久久久久久免费av| 两个人的视频大全免费| 国产精品一国产av| 亚洲av男天堂| a级毛片免费高清观看在线播放| 一二三四中文在线观看免费高清| 一本一本综合久久| 欧美3d第一页| 亚洲国产毛片av蜜桃av| 日日摸夜夜添夜夜添av毛片| 久久国产精品大桥未久av| 久久久久久久久久久久大奶| 精品久久久久久久久亚洲| 一本色道久久久久久精品综合| 久久久精品94久久精品| 国产精品久久久久久av不卡| 亚洲精品国产av成人精品| 日韩一区二区三区影片| 中文字幕人妻熟人妻熟丝袜美| 80岁老熟妇乱子伦牲交| 亚洲av二区三区四区| 亚洲av中文av极速乱| 99热网站在线观看| 亚洲av电影在线观看一区二区三区| videossex国产| 久久人人爽人人片av| 国产日韩欧美亚洲二区| 亚洲,欧美,日韩| 日本黄大片高清| 男人操女人黄网站| 一级片'在线观看视频| 中文字幕av电影在线播放| 精品人妻在线不人妻| 亚洲精华国产精华液的使用体验| 99视频精品全部免费 在线| 国产精品成人在线| 国产男人的电影天堂91| 亚洲,一卡二卡三卡| 最近2019中文字幕mv第一页| av线在线观看网站| 麻豆成人av视频| 纵有疾风起免费观看全集完整版| 国产极品粉嫩免费观看在线 | 黄色毛片三级朝国网站| 久久精品夜色国产| 免费黄频网站在线观看国产| 国产一区二区三区综合在线观看 | 久久精品国产亚洲av涩爱| av网站免费在线观看视频| 国产一区二区在线观看av| 国产综合精华液| 51国产日韩欧美| 老女人水多毛片| 麻豆成人av视频| 日韩一本色道免费dvd| 97超视频在线观看视频| 国产免费视频播放在线视频| 麻豆乱淫一区二区| 午夜福利视频在线观看免费| www.av在线官网国产| 国产极品天堂在线| 99久久中文字幕三级久久日本| 最近中文字幕高清免费大全6| 九色亚洲精品在线播放| 少妇精品久久久久久久| 精品酒店卫生间| 欧美日韩在线观看h| 欧美另类一区| a级毛片黄视频| 国产高清三级在线| 国产精品欧美亚洲77777| 97在线视频观看| 久久久久久人妻| 亚洲精品美女久久av网站| 九九久久精品国产亚洲av麻豆| 大香蕉97超碰在线| 蜜桃在线观看..| 人人妻人人澡人人看| 亚洲国产精品成人久久小说| 免费高清在线观看日韩| 精品久久久久久久久亚洲| 国产精品秋霞免费鲁丝片| 久久精品久久精品一区二区三区| 91久久精品国产一区二区成人| 少妇熟女欧美另类| 久久久欧美国产精品| 22中文网久久字幕| 欧美精品国产亚洲| 五月天丁香电影| 中文字幕久久专区| 夜夜看夜夜爽夜夜摸| 日本vs欧美在线观看视频| 91精品国产国语对白视频| 嘟嘟电影网在线观看| 欧美+日韩+精品| 91aial.com中文字幕在线观看| 9色porny在线观看| 欧美3d第一页| 国产成人freesex在线| 美女主播在线视频| 高清不卡的av网站| 久久久久久人妻| 亚洲av二区三区四区| 草草在线视频免费看| 美女福利国产在线| 亚洲图色成人| 高清午夜精品一区二区三区| 国产成人免费观看mmmm| 亚洲精品国产av成人精品| 国产一区有黄有色的免费视频| 亚洲欧美一区二区三区黑人 | 中文字幕精品免费在线观看视频 | 秋霞在线观看毛片| 尾随美女入室| 免费观看a级毛片全部| 黄色毛片三级朝国网站| 久久精品国产亚洲av天美| 亚洲av二区三区四区| 日日撸夜夜添| 五月天丁香电影| 赤兔流量卡办理| 久久99蜜桃精品久久| 精品少妇久久久久久888优播| 国产爽快片一区二区三区| 亚洲人成77777在线视频| 最近最新中文字幕免费大全7| 日韩亚洲欧美综合| 少妇的逼水好多| 亚洲精品日韩av片在线观看| 卡戴珊不雅视频在线播放| 成年美女黄网站色视频大全免费 | 欧美丝袜亚洲另类| 日韩,欧美,国产一区二区三区| 日韩不卡一区二区三区视频在线| 成人无遮挡网站| 久久狼人影院| 国产精品国产三级专区第一集| 日韩成人av中文字幕在线观看| 欧美成人午夜免费资源| 日本-黄色视频高清免费观看| 午夜久久久在线观看| 日本黄色片子视频| 国产免费现黄频在线看| 大码成人一级视频| 99热全是精品| 九草在线视频观看| 免费不卡的大黄色大毛片视频在线观看| av女优亚洲男人天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av欧美aⅴ国产| 观看美女的网站| 新久久久久国产一级毛片| 中文天堂在线官网| 观看美女的网站| 国产日韩欧美亚洲二区| 国产永久视频网站| 综合色丁香网| videos熟女内射| 视频在线观看一区二区三区| 国精品久久久久久国模美| 丝袜喷水一区| 国产精品偷伦视频观看了| 久久午夜福利片| 成人午夜精彩视频在线观看| 日韩中字成人| 老司机亚洲免费影院| 欧美人与善性xxx| 一本色道久久久久久精品综合| 亚洲精品成人av观看孕妇| 人人妻人人澡人人爽人人夜夜| 久久人人爽人人爽人人片va| 婷婷色av中文字幕| 免费观看av网站的网址| 美女xxoo啪啪120秒动态图| 极品少妇高潮喷水抽搐| 亚洲av中文av极速乱| 18禁在线无遮挡免费观看视频| 中文字幕久久专区| 少妇被粗大猛烈的视频| 中文乱码字字幕精品一区二区三区| 在线精品无人区一区二区三| 国产一区亚洲一区在线观看| 免费看av在线观看网站| 你懂的网址亚洲精品在线观看| 国产日韩一区二区三区精品不卡 | 少妇被粗大的猛进出69影院 | 在线观看一区二区三区激情| 国产精品一国产av| 女性被躁到高潮视频| 看免费成人av毛片| 有码 亚洲区| 熟妇人妻不卡中文字幕| 日韩av在线免费看完整版不卡| 久久午夜福利片| 亚洲性久久影院| 午夜激情久久久久久久| 亚洲人成网站在线播| 男女边摸边吃奶| 99九九在线精品视频| 人人妻人人爽人人添夜夜欢视频| xxx大片免费视频| 高清毛片免费看| 国产黄频视频在线观看| 一本一本综合久久| 日韩熟女老妇一区二区性免费视频| 男女高潮啪啪啪动态图| av在线老鸭窝| 国产极品粉嫩免费观看在线 | 乱人伦中国视频| 国产一区二区三区av在线| 简卡轻食公司| 在线天堂最新版资源| 亚洲,一卡二卡三卡| kizo精华| 一级a做视频免费观看| videos熟女内射| 久久精品国产自在天天线| 国产日韩欧美视频二区| 永久网站在线| 嫩草影院入口| 亚洲精品美女久久av网站| 制服人妻中文乱码| 日韩一区二区视频免费看| 一本色道久久久久久精品综合| 精品99又大又爽又粗少妇毛片| 蜜臀久久99精品久久宅男| 热99国产精品久久久久久7| 大香蕉久久网| 人人妻人人添人人爽欧美一区卜| 亚洲人成网站在线观看播放| 欧美性感艳星| 亚洲av二区三区四区| a级毛片黄视频| 亚洲激情五月婷婷啪啪| 麻豆精品久久久久久蜜桃| 午夜av观看不卡| 纯流量卡能插随身wifi吗| 国产精品秋霞免费鲁丝片| 亚洲少妇的诱惑av| 国产一区亚洲一区在线观看| 一级毛片aaaaaa免费看小|