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

    多普勒雷達(dá)中氣旋判據(jù)及算法的發(fā)展與應(yīng)用

    2016-11-16 03:16:32費(fèi)海燕周小剛王秀明
    關(guān)鍵詞:龍卷徑向速度氣旋

    費(fèi)海燕 周小剛 王秀明

    (中國(guó)氣象局氣象干部培訓(xùn)學(xué)院,北京 100081)

    多普勒雷達(dá)中氣旋判據(jù)及算法的發(fā)展與應(yīng)用

    費(fèi)海燕 周小剛 王秀明

    (中國(guó)氣象局氣象干部培訓(xùn)學(xué)院,北京 100081)

    由于超過90%的中氣旋伴隨著龍卷、強(qiáng)雹、大風(fēng)等強(qiáng)天氣發(fā)生,因此中氣旋的識(shí)別與應(yīng)用對(duì)于準(zhǔn)確、實(shí)時(shí)發(fā)布強(qiáng)天氣預(yù)警具有十分重要的意義。目前我國(guó)僅CINRAD WSR—98D SA/SB多普勒雷達(dá)內(nèi)置了美國(guó)強(qiáng)風(fēng)暴實(shí)驗(yàn)室開發(fā)的中氣旋算法,并有相應(yīng)的中氣旋產(chǎn)品,其他型號(hào)的多普勒雷達(dá)則只能從徑向速度圖上人工識(shí)別中氣旋?;仡櫫酥袣庑袚?jù)演變的三個(gè)階段,對(duì)當(dāng)前業(yè)務(wù)上使用的第三階段的中氣旋判據(jù)存在的問題進(jìn)行了分析;依據(jù)中氣旋算法發(fā)展歷程,介紹了中氣旋切變算法、中氣旋算法及中氣旋探測(cè)算法,并對(duì)當(dāng)前業(yè)務(wù)上使用的中氣旋算法得到的產(chǎn)品應(yīng)用及其存在問題作了討論。

    中氣旋判據(jù),中氣旋切變算法,中氣旋算法,中氣旋探測(cè)算法

    0 引言

    在20世紀(jì)中期,Brooks[1]就已認(rèn)識(shí)到在雷暴單體中可能存在尺度比龍卷要大的氣旋式渦旋。1963年,F(xiàn)ujita[2]將其正式命名為中氣旋(mesocyclone)。根據(jù)國(guó)外的統(tǒng)計(jì),90%以上的中氣旋會(huì)伴隨著龍卷、強(qiáng)雹、大風(fēng)等強(qiáng)天氣發(fā)生,因此中氣旋探測(cè)對(duì)于準(zhǔn)確、實(shí)時(shí)發(fā)布強(qiáng)天氣預(yù)警具有十分重要的意義。

    在多普勒雷達(dá)應(yīng)用之前,由于缺乏對(duì)速度場(chǎng)的探測(cè)資料,中氣旋研究進(jìn)展緩慢。1971年,美國(guó)國(guó)家海洋和大氣管理局(NOAA)的強(qiáng)風(fēng)暴實(shí)驗(yàn)室(NSSL)擁有了位于諾曼(Norman)的第一臺(tái)用于研究的S波段多普勒雷達(dá)。其后,依據(jù)Norman多普勒雷達(dá)所收集的數(shù)據(jù)集,NSSL對(duì)中氣旋的基本特征進(jìn)行了深入研究,先后提出了利用徑向速度對(duì)來識(shí)別中氣旋的判據(jù)及中氣旋客觀算法。

    中氣旋的判據(jù)發(fā)展經(jīng)歷了三個(gè)階段:使用徑向速度方位切變閾值和厚度閾值;使用與距雷達(dá)站距離有關(guān)的徑向速度方位切變閾值和厚度閾值;使用與距雷達(dá)站距離有關(guān)的旋轉(zhuǎn)速度閾值和厚度閾值。目前,國(guó)內(nèi)業(yè)務(wù)中最常使用的中氣旋判據(jù)是使用與距雷達(dá)站距離有關(guān)的旋轉(zhuǎn)速度和厚度閾值。

    中氣旋客觀算法的發(fā)展也經(jīng)歷了三個(gè)階段:中氣旋切變算法(mesocyclone shear algorithm)、中氣旋算法(mesocyclone algorithm)、中氣旋探測(cè)算法(mesocyclone detection algorithm)。目前在全美布網(wǎng)的165部WSR-88D S波段雷達(dá)build 9.0版本中內(nèi)置的是中氣旋算法,預(yù)期下一階段將會(huì)更新為中氣旋探測(cè)算法。我國(guó)CINRAD WSR-98D SA/SB目前使用的也是中氣旋算法。

    國(guó)內(nèi)在中氣旋判據(jù)和中氣旋算法得到的中氣旋產(chǎn)品的應(yīng)用上,業(yè)務(wù)人員也積累了一些相關(guān)經(jīng)驗(yàn)。個(gè)例分析[3]表明,與非超級(jí)單體龍卷相比,導(dǎo)致強(qiáng)龍卷的中氣旋底高明顯偏低,基本在1km以下。臺(tái)風(fēng)前部龍卷個(gè)例研究[4]發(fā)現(xiàn),中氣旋算法產(chǎn)品沒有報(bào)警,認(rèn)為中氣旋閾值應(yīng)結(jié)合當(dāng)?shù)靥攸c(diǎn)作調(diào)整。由中氣旋文本產(chǎn)品給出的最強(qiáng)切變值統(tǒng)計(jì)發(fā)現(xiàn),當(dāng)中氣旋最強(qiáng)切變達(dá)到15m·s-1·km-1以上極有可能出現(xiàn)25m·s-1以上的強(qiáng)風(fēng)或者龍卷天氣[5]。根據(jù)對(duì)中氣旋的旋轉(zhuǎn)速度統(tǒng)計(jì),64.3%以上的冰雹天氣的中氣旋達(dá)中等強(qiáng)度以上,50%以上的雷雨大風(fēng)天氣的中氣旋達(dá)中等強(qiáng)度以上,而強(qiáng)降水天氣達(dá)中等以上強(qiáng)度僅為36.4%[6]。

    國(guó)內(nèi)多普勒雷達(dá)應(yīng)用研究起步較晚,且國(guó)內(nèi)大部分C波段多普勒雷達(dá)應(yīng)用軟件中無中氣旋算法,因此業(yè)務(wù)上由中氣旋判據(jù)和由中氣旋算法產(chǎn)品來識(shí)別中氣旋是并存的。

    1 中氣旋判據(jù)發(fā)展的三個(gè)階段

    由于多普勒雷達(dá)僅能探測(cè)到一維的徑向風(fēng)場(chǎng),因此必須對(duì)中氣旋流場(chǎng)作一定假設(shè),才能在徑向速度場(chǎng)上給予合理解釋。故通常假設(shè)中氣旋流場(chǎng)滿足Rankine渦旋條件,即:

    式中,vt表示切向速度,Vt表示核半徑r=Rt時(shí)的最大切向速度。在核區(qū)內(nèi)(r≤Rt),切向速度隨半徑線性增加;在核區(qū)外,速度與半徑成反比。對(duì)氣旋式旋轉(zhuǎn),Vt為正;對(duì)反氣旋式旋轉(zhuǎn),Vt為負(fù)。

    中氣旋為氣旋式旋轉(zhuǎn)。由于多普勒雷達(dá)平面掃描是以順時(shí)針方向進(jìn)行的(正北為0°方位角,順時(shí)針方位角增加到360°)。因此,在多普勒雷達(dá)平面掃描徑向速度場(chǎng)上,中氣旋就表示為同一距離圈上,相隔一定方位角距離(即直徑)的最大負(fù)速度與最大正速度(簡(jiǎn)稱最大正負(fù)速度對(duì)或徑向速度對(duì))。傳統(tǒng)意義上,中氣旋可以定義為具有風(fēng)暴尺度(2~10km)的Rankine 渦旋特征的多普勒徑向速度對(duì)。徑向速度對(duì)的值越大或徑向速度對(duì)之間的距離越小,說明中氣旋的強(qiáng)度越強(qiáng)。

    1.1 第一階段的中氣旋判據(jù)

    從傳統(tǒng)意義上的中氣旋定義到中氣旋定量識(shí)別判據(jù),主要是考慮到中氣旋為三維旋轉(zhuǎn)系統(tǒng),加入了對(duì)中氣旋厚度要求的限制。中氣旋識(shí)別判據(jù)的第一階段即是使用徑向速度方位切變值和厚度值作為中氣旋識(shí)別判據(jù)。

    該判據(jù)是Donaldson[7]在1970年提出的。Donaldson[7]使用美國(guó)空軍5.4cm多普勒氣象雷達(dá),分析了1968年春季美國(guó)Massachusetts東部地區(qū)一次伴隨大范圍強(qiáng)雹和龍卷過程的中氣旋演變,首次提出中氣旋定量識(shí)別判據(jù)閾值為:徑向速度方位切變值≥5×10-3s-1;其伸展厚度大于2Rt。1976年,Burgess[8]利用NSSL的Norman多普勒雷達(dá)在1971—1975年的每年春季收集到的37個(gè)Oklahoma中氣旋(62%伴有龍卷)進(jìn)行特征統(tǒng)計(jì)研究時(shí),降低了Donaldson判據(jù)的厚度值要求,采用的中氣旋識(shí)別判據(jù)閾值為:徑向速度方位切變值≥5×10-3s-1;其伸展厚度大于3km。

    1.2 第二階段的中氣旋判據(jù)

    第二階段的中氣旋判據(jù)是使用與距雷達(dá)距離有關(guān)的徑向速度方位切變閾值和厚度閾值識(shí)別中氣旋。這是1977年NSSL研究人員在Norman多普勒雷達(dá)(注:Norman多普勒雷達(dá)徑向速度探測(cè)范圍為345km)進(jìn)行的春季聯(lián)合多普勒業(yè)務(wù)試驗(yàn)(Joint Doppler Operational Project,JDOP)中提出來的。主要是考慮到受雷達(dá)采樣方式的影響,距雷達(dá)越遠(yuǎn)分辨率越低,因此,采用的中氣旋識(shí)別判據(jù)閾值為:距雷達(dá)230km范圍內(nèi),方位切變值≥5×10-3s-1,其伸展厚度大于3km;距雷達(dá)230km范圍外,切變值≥1×10-3s-1,無厚度要求[9]。此后,針對(duì)中氣旋的研究一直使用該識(shí)別判據(jù)。

    1.3 第三階段的中氣旋判據(jù)

    第三階段的中氣旋判據(jù)是使用與距雷達(dá)距離有關(guān)的旋轉(zhuǎn)速度閾值和厚度閾值識(shí)別中氣旋。這是1987年NSSL研究人員在多普勒/閃電(Doppler/Lightning,DOPLIGHT)試驗(yàn)中首次提出來的,其使用旋轉(zhuǎn)速度(即Vt)代替方位切變值作為判據(jù)。中氣旋識(shí)別判據(jù)閾值為:距雷達(dá)150km范圍內(nèi),旋轉(zhuǎn)速度≥15m·s-1,其伸展厚度大于3km;距雷達(dá)150~230km范圍,旋轉(zhuǎn)速度≥11m·s-1,無厚度要求[10]。此后,美國(guó)預(yù)報(bào)業(yè)務(wù)人員開始使用旋轉(zhuǎn)速度作為定量識(shí)別中氣旋的重要參量。

    1997年,Andra[11]根據(jù)美國(guó)俄克拉荷馬(Oklahoma)中部伴有強(qiáng)雹和龍卷的中氣旋進(jìn)行統(tǒng)計(jì),給出了距雷達(dá)不同距離圈上觀測(cè)到的旋轉(zhuǎn)速度和中氣旋強(qiáng)度的關(guān)系,將中氣旋判據(jù)進(jìn)行了細(xì)化,根據(jù)旋轉(zhuǎn)速度和距雷達(dá)距離,將其分成了弱切變、弱中氣旋、中等強(qiáng)度中氣旋、強(qiáng)中氣旋四個(gè)區(qū)域(圖1)。

    圖1 中氣旋強(qiáng)度閾值(改自Andra[11])Fig. 1 Mesocyclone strength threshold

    圖1 是目前我國(guó)業(yè)務(wù)中最常使用的中氣旋判據(jù)依據(jù)。對(duì)此圖的使用要注意以下幾點(diǎn):1)圖1是根據(jù)美國(guó)Oklahoma中部伴有強(qiáng)雹和龍卷的中氣旋統(tǒng)計(jì)結(jié)果繪出的,按照Burgess等[12]的統(tǒng)計(jì)結(jié)論,Oklahoma中部地區(qū)伴有強(qiáng)雹或龍卷的中氣旋成熟階段時(shí)的核半徑平均為2.75km,旋轉(zhuǎn)速度平均為23.3m·s-1。由圖1可知,該旋轉(zhuǎn)速度對(duì)應(yīng)為強(qiáng)中氣旋;而根據(jù)我國(guó)的統(tǒng)計(jì)[13],中氣旋成熟階段時(shí)的平均核半徑為3.085km,平均旋轉(zhuǎn)速度為18.945m·s-1。由圖1可知,在距雷達(dá)150km范圍內(nèi),僅對(duì)應(yīng)為中等強(qiáng)度中氣旋。因此在我國(guó)使用圖1時(shí),要考慮到我國(guó)中氣旋的平均強(qiáng)度比Oklahoma中部地區(qū)伴有強(qiáng)雹或龍卷的中氣旋要弱的事實(shí),適當(dāng)降低對(duì)旋轉(zhuǎn)速度值的要求;2)圖1是主觀估計(jì)中氣旋強(qiáng)度的工具,對(duì)無經(jīng)驗(yàn)的預(yù)報(bào)員,可以起指導(dǎo)作用。但要注意圖中中氣旋的特定直徑是6.5km,對(duì)直徑明顯小于或大于6.5km的中氣旋,該圖會(huì)過低或過高地估計(jì)了中氣旋的強(qiáng)度;3)根據(jù)NSSL中氣旋第三階段判據(jù),中氣旋要有一定的垂直伸展厚度,而圖1僅考慮了旋轉(zhuǎn)速度,并沒有考慮厚度。由于單層中氣旋有可能是沒有意義的,因此,在雷達(dá)探測(cè)單仰角上旋轉(zhuǎn)速度滿足閾值后,還要確定一下所探測(cè)到的旋轉(zhuǎn)速度區(qū)是否滿足垂直伸展厚度大于3km的要求。

    2 中氣旋客觀算法的發(fā)展

    由第一節(jié)可知,Donaldson[7]首次提出了用單多普勒雷達(dá)識(shí)別中氣旋的判據(jù),即徑向速度方位切變值≥5×10-3s-1,垂直伸展厚度等于或大于其直徑的2倍。其后,NSSL算法開發(fā)人員對(duì)Donaldson的判據(jù)進(jìn)行了改進(jìn),并開發(fā)了中氣旋切變算法來自動(dòng)識(shí)別是否存在滿足閾值的中氣旋切變區(qū)域。目前我國(guó)CINRAD WSR-98D SA/SB使用的中氣旋算法是NSSL對(duì)中氣旋切變算法的改進(jìn)。因此,本節(jié)按照NSSL對(duì)中氣旋客觀算法改進(jìn)的過程,分別介紹中氣旋切變算法、中氣旋算法、中氣旋探測(cè)算法。

    2.1 中氣旋切變算法

    圖2為中氣旋切變算法流程圖,圖中一些術(shù)語定義如下:

    1)類型矢徑(pattern vector):指在常定距離圈上尋找到的速度增加的長(zhǎng)度。

    2)角動(dòng)量:類型矢徑的速度差和方位角距離的乘積。

    3)切變:類型矢徑的速度差除以方位角距離。

    4)2D特征:一個(gè)類型矢徑與其他類型矢徑在方位角和距離圈距離小于閾值時(shí),則可以組成一個(gè)2D特征。

    5)對(duì)稱特征:2D特征切向和徑向長(zhǎng)度之比在閾值范圍內(nèi),則認(rèn)為是對(duì)稱特征。

    圖2 中氣旋切變算法流程圖(引自Zrnic等[14])Fig. 2 Flow chart of mesocyclone shear algorithm

    1985年,Zrnic等[14]依據(jù)1977年春季美國(guó)聯(lián)合多普勒業(yè)務(wù)項(xiàng)目(JDOP)觀測(cè),利用位于Norman的多普勒雷達(dá)所觀測(cè)到的Oklahoma中部40個(gè)中氣旋例子(伴隨強(qiáng)雹和強(qiáng)龍卷天氣的中氣旋),給出了中氣旋切變算法中參數(shù)的閾值:如低切變?nèi)?m·s-1·km-1,高切變?nèi)?m·s-1·km-1,低角動(dòng)量取為50m·s-1·km-1,高角動(dòng)量取為150m·s-1·km-1。這些值是根據(jù)這40個(gè)中氣旋主觀確定的,原則上使用這些閾值,可以保證40個(gè)中氣旋沒有一個(gè)被遺漏。

    對(duì)比圖1發(fā)現(xiàn),低切變?nèi)?m·s-1·km-1,意味著強(qiáng)度在弱切變以上的中氣旋,在中氣旋切變算法中都給予保留。這樣算法識(shí)別的中氣旋能完全涵蓋由圖1識(shí)別的弱切變以上的中氣旋。

    2.2 中氣旋算法及產(chǎn)品應(yīng)用

    2.2.1 中氣旋算法

    NSSL對(duì)中氣旋切變算法進(jìn)行改進(jìn)后,將其命名為中氣旋算法,1996年秋季開始正式用于全美布網(wǎng)的165部S波段多普勒雷達(dá)。目前我國(guó)CINRAD WSR-98D SA/SB使用的即是中氣旋算法。

    中氣旋切變算法考慮的僅是單層中氣旋,中氣旋算法增加了對(duì)垂直連續(xù)性的檢驗(yàn),即考慮中氣旋有一定的伸展厚度。中氣旋算法第一步是由中氣旋切變算法探測(cè)2D特征(對(duì)稱或不對(duì)稱);第二步是通過對(duì)垂直連續(xù)性的檢驗(yàn)將中氣旋切變分成三種類型。其流程圖見圖3。

    圖3 中氣旋算法流程圖Fig.3 Flow chart of mesocyclone algorithm

    2.2.2 中氣旋產(chǎn)品應(yīng)用

    由圖3可見,中氣旋算法通過對(duì)2D特征做垂直相關(guān)性檢驗(yàn),最終得到三種類型:非相關(guān)切變、三維切變、中氣旋。在目前的CINRAD WSR-98D SA/SB中氣旋圖形產(chǎn)品(M,60號(hào))中,中氣旋用深黃色的圓圈(線的寬度為4個(gè)像素)顯示,并帶有離它最近的被識(shí)別的風(fēng)暴單體的標(biāo)識(shí)號(hào);三維切變用淺黃色的圓圈(線的寬度為1個(gè)像素)顯示,不帶有風(fēng)暴單體的標(biāo)識(shí)號(hào);非相關(guān)切變?cè)趫D形產(chǎn)品上不顯示,但在中氣旋文本產(chǎn)品(M,60號(hào);中氣旋文本產(chǎn)品是中氣旋圖形產(chǎn)品的匹配產(chǎn)品,共用一個(gè)產(chǎn)品號(hào))中有顯示。

    中氣旋圖形產(chǎn)品可以疊加到反射率因子或徑向速度圖形產(chǎn)品上,與強(qiáng)單體相伴隨的中氣旋伴隨著龍卷、強(qiáng)雹、大風(fēng)等強(qiáng)天氣發(fā)生,要充分引起注意。圖4出示了2013年3月23日14:53(北京時(shí)間)1.5°仰角百色雷達(dá)反射率因子圖疊加中氣旋圖形產(chǎn)品(黃色圓圈),此超級(jí)單體風(fēng)暴維持了較長(zhǎng)時(shí)間的中氣旋,并造成沿途經(jīng)過的田林縣定安鎮(zhèn)、八渡鄉(xiāng)、八桂鄉(xiāng)出現(xiàn)了直徑為20~30mm的冰雹、強(qiáng)風(fēng)和雷電等強(qiáng)對(duì)流天氣 。

    圖4 2013年3月23日百色雷達(dá)強(qiáng)反射率因子疊加中氣旋(黃色圓圈)Fig.4 Strong reflectivity with mesocyclone product (yellow circle) in Baise radar on 23 March 2013

    在中氣旋圖形產(chǎn)品上并不能確定中氣旋最強(qiáng)2D特征的高度,但在中氣旋文本產(chǎn)品中會(huì)顯示最強(qiáng)切變值及其高度??梢岳梦谋井a(chǎn)品中的最強(qiáng)切變值及其高度進(jìn)一步區(qū)分強(qiáng)天氣類型。

    2013年3月23日午后,百色市田林縣、右江區(qū)、田陽縣、田東縣、平果縣的部分地區(qū)先后出現(xiàn)雷暴、冰雹、雷雨大風(fēng)等災(zāi)害性天氣。根據(jù)百色SB雷達(dá)探測(cè),導(dǎo)致田林降雹的強(qiáng)風(fēng)暴維持了一個(gè)持續(xù)較長(zhǎng)時(shí)間的中氣旋,圖5即是中氣旋最強(qiáng)切變高度演變圖。由圖5可見,強(qiáng)中氣旋的最強(qiáng)切變維持在0℃層左右,這里稱為中層中氣旋,意味著旋轉(zhuǎn)上升氣流超過0℃層。由冰雹形成和增長(zhǎng)的概念模型,上升氣流超過0℃層高度且持續(xù)時(shí)間足夠長(zhǎng),說明冰雹胚胎在過冷卻水累積區(qū)生長(zhǎng),冰雹則可以增長(zhǎng)到很大的尺寸。

    圖5 2013年3月23日14:40—15:43百色雷達(dá)探測(cè)降雹單體中氣旋最強(qiáng)切變高度Fig.5 The height of the strongest shear in the hail cell detected by Baise radar at 14:40 to 15:43 BT on 23 March 2013

    2013年7月7日午后,一個(gè)伴有短時(shí)強(qiáng)降水和龍卷產(chǎn)生的超級(jí)單體風(fēng)暴維持了一個(gè)持續(xù)較長(zhǎng)時(shí)間的中氣旋,圖6是中氣旋最強(qiáng)切變高度演變圖。由圖6可見,強(qiáng)中氣旋的最強(qiáng)切變維持在0℃層距地高度一半左右,可以稱為低層中氣旋,可以排除降雹可能。在環(huán)境中邊界層相對(duì)濕度較高和 0~1 km垂直風(fēng)切變較強(qiáng)時(shí),如果低層中氣旋強(qiáng)度強(qiáng)、維持時(shí)間長(zhǎng),龍卷發(fā)生的可能則較大。

    圖6 2013年7月7日15:16—16:59南京雷達(dá)探測(cè)的龍卷單體中氣旋最強(qiáng)切變高度Fig.6 The height of the strongest shear in the tornado cell detected by Nanjing radar at 15:16 to 16:59 BT on 7 J uly, 2013

    2.2.3 中氣旋產(chǎn)品應(yīng)用中需要注意的問題

    1)注意剔除退速度模糊引起的虛假中氣旋。在中氣旋算法中,沒有去除風(fēng)暴運(yùn)動(dòng),因此對(duì)退模糊算法有很大的依賴性,不正確的退模糊算法會(huì)導(dǎo)致不正確的中氣旋識(shí)別。如圖7所示的2006年4月9日湖南永州的例子,由于退速度模糊的影響,圖中標(biāo)識(shí)號(hào)為Do的深黃色圓圈(算法識(shí)別的中氣旋)是虛假的中氣旋。

    2)目前算法中的大部分缺省參數(shù)來源于中氣旋方位切變算法開發(fā)使用的參數(shù)。由于資料集只有幾個(gè)風(fēng)暴過程時(shí)間只有幾小時(shí),且是20世紀(jì)70年代Oklahoma中部春季伴有強(qiáng)冰雹和龍卷的中氣旋參數(shù)。因此在我國(guó)使用時(shí)可以適當(dāng)對(duì)缺省參數(shù)做修改。如對(duì)于我國(guó)臺(tái)風(fēng)前部龍卷,可以適當(dāng)調(diào)低缺省參數(shù)中的類型矢徑(TPV)值。

    圖7 2006年4月9日湖南永州疊加中氣旋產(chǎn)品的徑向速度圖Fig.7 Radial velocity and mesocyclone product in Yongzhou, Hunan on 9 April 2006

    2.3 中氣旋探測(cè)算法

    中氣旋算法參數(shù)來源于中氣旋切變算法參數(shù),是以Norman多普勒雷達(dá)觀測(cè)到的中氣旋資料集建立的。隨著美國(guó)多普勒雷達(dá)的布網(wǎng),NSSL發(fā)現(xiàn)一些地區(qū)產(chǎn)生了龍卷等強(qiáng)天氣的風(fēng)暴尺度(2~10km)渦旋會(huì)被中氣旋算法漏報(bào)。因此,NSSL又設(shè)計(jì)開發(fā)出了一個(gè)新算法,稱為中氣旋探測(cè)算法[15]。根據(jù)算法過程,給出下列流程圖(圖8)。

    圖8 中氣旋探測(cè)算法流程圖Fig.8 The flow chart of mesocyclone detection algorithm

    與目前使用的中氣旋算法相比,中氣旋探測(cè)算法使用了中氣旋強(qiáng)度閾值來代替切變閾值;并使用多個(gè)不同距離的強(qiáng)度閾值,以識(shí)別更多的2D特征;增加了時(shí)間連續(xù)性測(cè)試,能夠跟蹤并預(yù)報(bào)中氣旋的移動(dòng)路徑。圖形產(chǎn)品中增加了隨時(shí)間變化的最大切變高度顯示等特征(圖5、圖6是根據(jù)目前中氣旋文本產(chǎn)品的最大切變高度值和探空的0℃層高度做的,中氣旋探測(cè)算法的圖形產(chǎn)品則將會(huì)直接顯示隨時(shí)間變化的最大切變高度),這樣可以使預(yù)報(bào)員更快地判斷中氣旋的強(qiáng)度值和高度值,為相應(yīng)強(qiáng)天氣預(yù)警提供寶貴的提前時(shí)間。

    3 結(jié)語

    中氣旋探測(cè)與強(qiáng)天氣預(yù)警密切相關(guān)。中氣旋的判據(jù)發(fā)展經(jīng)歷了三個(gè)階段,目前國(guó)內(nèi)業(yè)務(wù)中最常使用的中氣旋判據(jù)是其中的第三階段,即使用與距雷達(dá)站距離有關(guān)的旋轉(zhuǎn)速度和厚度閾值。對(duì)沒有內(nèi)置中氣旋算法的多普勒雷達(dá),預(yù)報(bào)業(yè)務(wù)人員依據(jù)圖1來識(shí)別中氣旋及其強(qiáng)度,要注意其閾值不一定完全適用于國(guó)內(nèi)中氣旋的實(shí)際狀況,因此需有更多的統(tǒng)計(jì)結(jié)果,才能由權(quán)威部門對(duì)圖1的閾值做修訂。此外,在雷達(dá)探測(cè)單仰角上旋轉(zhuǎn)速度滿足閾值后,還要確定所探測(cè)到的旋轉(zhuǎn)速度區(qū)是否滿足垂直伸展厚度大于3km的要求。對(duì)滿足垂直伸展厚度的真實(shí)中氣旋要立即發(fā)布強(qiáng)天氣警報(bào)。

    中氣旋客觀算法的發(fā)展也經(jīng)歷了三個(gè)階段:中氣旋切變算法、中氣旋算法和中氣旋探測(cè)算法。對(duì)目前有中氣旋算法的CINRAD WSR-98D SA/SB多普勒雷達(dá),對(duì)與強(qiáng)單體相伴隨的真實(shí)中氣旋要引起高度注意。由于目前我國(guó)多普勒雷達(dá)業(yè)務(wù)中使用的速度退模糊算法效果不夠顯著,因此對(duì)算法識(shí)別的中氣旋要人工剔除由退速度模糊算法引起的虛假中氣旋。同時(shí)要充分關(guān)注文本產(chǎn)品中切變值大小和高度。根據(jù)中氣旋文本產(chǎn)品中的最大2D切變高度與環(huán)境0℃層高度的比較,可以將中氣旋分成中層中氣旋和低層中氣旋,通常在環(huán)境邊界層相對(duì)濕度較高和 0~1km強(qiáng)垂直風(fēng)切變時(shí),如果低層中氣旋強(qiáng)度強(qiáng)、維持時(shí)間長(zhǎng),則要優(yōu)先發(fā)龍卷預(yù)警;對(duì)直徑較大的低層中氣旋,優(yōu)先發(fā)短時(shí)強(qiáng)降水預(yù)警;對(duì)中層強(qiáng)中氣旋,優(yōu)先發(fā)冰雹預(yù)警。

    [1] Brooks E M. T e tornado cyclone. Weatherwise, 1949, 2: 32-33.

    [2] Fujita T T. Analytical mesometeorology: a review. Meteor Monogr,1963, 5: 77-125.

    [3] 鄭媛媛, 朱紅芳, 方翔, 等. 強(qiáng)龍卷超級(jí)單體風(fēng)暴特征分析與預(yù)警研究. 高原氣象, 2009, 28(3): 617-625.

    [4] 蔣義芳, 吳海英, 沈樹勤, 等. 0808號(hào)臺(tái)風(fēng)鳳凰前部龍卷的環(huán)境場(chǎng)和雷達(dá)回波分析. 氣象, 2009, 35(4): 68-75.

    [5] 方翀,鄭媛媛. 新一代天氣雷達(dá)中氣旋產(chǎn)品特征值統(tǒng)計(jì)和個(gè)例分析.氣象, 2007, 33(11): 16-20.

    [6] 馮晉勤, 湯達(dá)章, 俞小鼎, 等. 新一代天氣雷達(dá)中氣旋識(shí)別產(chǎn)品的統(tǒng)計(jì)分析. 氣象, 2010, 36(8): 47-52.

    [7] Donaldson R J. Vortex signature recognition by a Doppler radar. J Appl Meteor, 1970, 9: 661-670.

    [8] Burgess D W. Single Doppler radar vortex recognition. Part I: Mesocyclone signatures. Preprints, 17th Conf on Radar Meteorology, Seattle, WA. Amer Meteor Soc, 1976.

    [9] JDOP Staff . Final report on the Joint Doppler Operational Project(JDOP), 1979.

    [10] Forsyth D E, Burgess D W, Mooney L E, et al. DOPLIGHT 87 program summary. NOAA Tech. Memo. ERL NSSL-101. National Severe Storms Laboratory, Norman, OK, 1989.

    [11] Andra D L. The origin and evolution of the WSR-88 D mesocyclone recognition nomogram. Preprints, 28th Conf on Radar Meteorology, Austin, TX. Amer Meteor Soc, 1997.

    [12] Burgess D W, Wood V T, Brown R A. Mesocyclone evolution statistics. Preprints, 12th Conf on Severe Local Storms, San Antonio, TX. Amer Meteor Soc, 1982.

    [13] 周小剛, 王秀明, 俞小鼎, 等. 逾量旋轉(zhuǎn)動(dòng)能在區(qū)分我國(guó)龍卷與非龍卷中氣旋中的應(yīng)用. 高原氣象, 2012, 31(1): 137-143.

    [14] Zrnic D S, Burgess D W, Hennington L D. Automatic detection of mesocyclonic shear with Doppler radar. J Atmos Oceanic Technol,1985, 2: 425-438.

    [15] Stumpf G J, Witt A, Mitchell E D W, et al. T e National Severe Storms Laboratory Mesocyclone Detection Algorithm for the WSR-88D. Wea Forecasting, 1998, 13: 304-326.

    The Development and Application of Doppler Radar Mesocyclone Criterion and Algorithm

    Fei Haiyan, Zhou Xiaogang, Wang Xiuming
    (China Meteorological Administration Training Centre, Beijing 100081)

    Since more than 90 % of the mesocyclones are accompanied by severe weather such as tornadoes, severe hails, strong winds etc., it is a great signifi cance to identify the mesocyclones and apply to publishing the severe weather warning in real-time for more accurate. At present, only the CINRAD WSR-98D SA / SB in China has mesocyclone algorithm, which was developed by US National Severe Storms Laboratory, and mesocyclone products. While other types of Doppler radar have to recognize mesocyclone by using radial velocity maps artifi cially only. This paper reviews the evolution of the mesocyclone criterion for three stages in the operational application, and analyzes some issues in the third stage. Based on the development of the mesocyclone algorithm, this paper introduces the mesocyclone shear algorithm, mesocyclone algorithm and mesocyclone detection algorithm,and discusses the operational applications and the issues of products from the mesocyclone algorithm.

    mesocyclone, criterion, algorithm, shear, detection

    10.3969/j.issn.2095-1973.2016.05.003

    2015年6月3日;

    2015年9月11日

    費(fèi)海燕(1983—),Email: feihy2009@163.com

    資助信息:國(guó)家自然科學(xué)基金項(xiàng)目(41475042);江蘇省氣象科學(xué)研究所北極閣基金(BJG201305)

    猜你喜歡
    龍卷徑向速度氣旋
    2003年3月北大西洋上兩個(gè)爆發(fā)性氣旋的“吞并”過程及發(fā)展機(jī)制分析?
    氨還原劑噴槍中保護(hù)氣旋流噴射作用的研究
    能源工程(2021年1期)2021-04-13 02:05:56
    海南龍卷統(tǒng)計(jì)特征
    我國(guó)龍卷預(yù)警試驗(yàn)的前沿陣地
    ——訪佛山市龍卷風(fēng)研究中心炎利軍主任
    龍卷的雷達(dá)探測(cè)研究進(jìn)展
    非圓形光纖研究進(jìn)展
    航空兵器(2017年6期)2018-01-24 15:00:10
    龍卷過程中雷達(dá)合成切變產(chǎn)品的分析
    臺(tái)風(fēng)威馬遜造成云南文山州強(qiáng)降水天氣雷達(dá)回波分析
    北太平洋上一個(gè)爆發(fā)性氣旋族的結(jié)構(gòu)分析?
    2014年3月大西洋上一個(gè)爆發(fā)性氣旋的研究
    免费看av在线观看网站| 日日啪夜夜撸| 国产一区二区激情短视频| 六月丁香七月| 美女高潮的动态| 天堂av国产一区二区熟女人妻| 小说图片视频综合网站| 日日撸夜夜添| 岛国在线免费视频观看| 国产亚洲欧美98| 亚洲av不卡在线观看| av专区在线播放| 国产白丝娇喘喷水9色精品| 亚洲欧美日韩无卡精品| 国产一区二区三区在线臀色熟女| 1024手机看黄色片| 看片在线看免费视频| 国产午夜福利久久久久久| 亚洲欧洲国产日韩| 好男人视频免费观看在线| 如何舔出高潮| 一个人看视频在线观看www免费| 午夜福利高清视频| 亚洲精品国产成人久久av| 国产真实乱freesex| 亚洲人成网站高清观看| 日本黄色视频三级网站网址| 高清日韩中文字幕在线| 欧美xxxx黑人xx丫x性爽| 一卡2卡三卡四卡精品乱码亚洲| 麻豆成人av视频| 欧美又色又爽又黄视频| 舔av片在线| 欧美xxxx性猛交bbbb| 国产探花在线观看一区二区| 男的添女的下面高潮视频| 99热只有精品国产| 国产精品嫩草影院av在线观看| 国产在视频线在精品| 国语自产精品视频在线第100页| eeuss影院久久| 一个人免费在线观看电影| 成年免费大片在线观看| 日本熟妇午夜| 男女做爰动态图高潮gif福利片| 亚洲av熟女| 国产一级毛片在线| h日本视频在线播放| 亚洲自拍偷在线| 淫秽高清视频在线观看| 久久精品国产亚洲网站| 国产成人freesex在线| 亚洲欧美日韩无卡精品| 免费搜索国产男女视频| 最近最新中文字幕大全电影3| 国产精品一区二区在线观看99 | 国产成人午夜福利电影在线观看| 亚洲乱码一区二区免费版| 日韩强制内射视频| 干丝袜人妻中文字幕| 97超碰精品成人国产| 日本免费一区二区三区高清不卡| 最近视频中文字幕2019在线8| 三级毛片av免费| 美女被艹到高潮喷水动态| 精品一区二区三区人妻视频| 综合色av麻豆| 国产日本99.免费观看| 国产激情偷乱视频一区二区| 久久精品人妻少妇| 免费大片18禁| 国产久久久一区二区三区| 婷婷精品国产亚洲av| 丰满的人妻完整版| 国产老妇伦熟女老妇高清| 欧美一级a爱片免费观看看| 精品一区二区三区人妻视频| 成人漫画全彩无遮挡| 国产亚洲欧美98| 国产精品99久久久久久久久| 热99在线观看视频| 99久久精品国产国产毛片| 国产色爽女视频免费观看| 精品熟女少妇av免费看| 日韩强制内射视频| 99热这里只有精品一区| 永久网站在线| 最近视频中文字幕2019在线8| 美女高潮的动态| 日本熟妇午夜| 美女国产视频在线观看| 好男人在线观看高清免费视频| 欧美最黄视频在线播放免费| 国产精品无大码| 日日干狠狠操夜夜爽| 国产91av在线免费观看| 成年女人看的毛片在线观看| 有码 亚洲区| 国产成人aa在线观看| 国产成人午夜福利电影在线观看| 男人舔女人下体高潮全视频| 午夜免费激情av| 日韩视频在线欧美| 2021天堂中文幕一二区在线观| 欧美一区二区亚洲| 直男gayav资源| 美女脱内裤让男人舔精品视频 | 内地一区二区视频在线| 亚洲第一区二区三区不卡| 天堂√8在线中文| 一进一出抽搐动态| 麻豆成人av视频| 久久精品久久久久久噜噜老黄 | 国产综合懂色| 欧美+亚洲+日韩+国产| 国产一区二区三区在线臀色熟女| 亚洲精品乱码久久久v下载方式| 国产老妇女一区| 精品午夜福利在线看| 在线免费十八禁| 亚洲精华国产精华液的使用体验 | 尤物成人国产欧美一区二区三区| 身体一侧抽搐| 精品免费久久久久久久清纯| 久久人妻av系列| 日韩欧美 国产精品| 久久精品综合一区二区三区| 国产69精品久久久久777片| 欧美激情久久久久久爽电影| 国产精品一区二区三区四区久久| 一进一出抽搐动态| 成年版毛片免费区| 床上黄色一级片| 夜夜夜夜夜久久久久| 大香蕉久久网| 哪里可以看免费的av片| 极品教师在线视频| 插逼视频在线观看| 久久精品影院6| 99热全是精品| 一级黄色大片毛片| 亚洲精品国产成人久久av| 中出人妻视频一区二区| 美女脱内裤让男人舔精品视频 | 精品少妇黑人巨大在线播放 | 亚洲精品影视一区二区三区av| 最近最新中文字幕大全电影3| 一进一出抽搐gif免费好疼| a级毛片a级免费在线| 国产单亲对白刺激| 1000部很黄的大片| 一个人观看的视频www高清免费观看| 一本精品99久久精品77| av在线观看视频网站免费| 两个人视频免费观看高清| 久久久a久久爽久久v久久| 日韩av不卡免费在线播放| 午夜亚洲福利在线播放| 最近最新中文字幕大全电影3| 日韩人妻高清精品专区| 亚洲婷婷狠狠爱综合网| 亚洲三级黄色毛片| 少妇被粗大猛烈的视频| 人妻少妇偷人精品九色| 亚洲色图av天堂| 黑人高潮一二区| 免费在线观看成人毛片| 国产毛片a区久久久久| 久久中文看片网| 国产美女午夜福利| 日本成人三级电影网站| 舔av片在线| 别揉我奶头 嗯啊视频| 寂寞人妻少妇视频99o| 亚洲精品国产成人久久av| 欧美另类亚洲清纯唯美| 91午夜精品亚洲一区二区三区| 亚洲av第一区精品v没综合| 天堂影院成人在线观看| 亚洲精品日韩在线中文字幕 | 国内精品美女久久久久久| 激情 狠狠 欧美| 日本三级黄在线观看| 亚洲成人久久性| 免费看美女性在线毛片视频| 日韩大尺度精品在线看网址| 久久国产乱子免费精品| 亚洲四区av| 国产伦理片在线播放av一区 | 亚洲成a人片在线一区二区| 欧美不卡视频在线免费观看| 精品午夜福利在线看| 日本黄色片子视频| av在线蜜桃| 免费观看人在逋| 国产在视频线在精品| 欧美又色又爽又黄视频| 18禁在线无遮挡免费观看视频| 99热精品在线国产| www日本黄色视频网| 亚洲av免费高清在线观看| 国产精品99久久久久久久久| 国产一区二区激情短视频| 一卡2卡三卡四卡精品乱码亚洲| 色哟哟哟哟哟哟| 日本爱情动作片www.在线观看| 国产日韩欧美在线精品| 精华霜和精华液先用哪个| 97热精品久久久久久| 中文字幕av在线有码专区| 男人舔女人下体高潮全视频| 国产视频内射| 12—13女人毛片做爰片一| 搡老妇女老女人老熟妇| a级一级毛片免费在线观看| 3wmmmm亚洲av在线观看| 少妇人妻精品综合一区二区 | 综合色丁香网| av专区在线播放| 最新中文字幕久久久久| 夜夜夜夜夜久久久久| av.在线天堂| 免费大片18禁| 欧美变态另类bdsm刘玥| 精品一区二区免费观看| 国产精品蜜桃在线观看 | 国产成人影院久久av| 亚洲欧美成人精品一区二区| 国产成人精品婷婷| 久久精品国产鲁丝片午夜精品| 日韩一本色道免费dvd| 一个人观看的视频www高清免费观看| a级一级毛片免费在线观看| 99久久精品热视频| 久久精品国产99精品国产亚洲性色| 精品熟女少妇av免费看| 色播亚洲综合网| 久久这里有精品视频免费| 久久婷婷人人爽人人干人人爱| 非洲黑人性xxxx精品又粗又长| 在线a可以看的网站| 久久久精品欧美日韩精品| 天堂影院成人在线观看| 男女下面进入的视频免费午夜| 国产麻豆成人av免费视频| h日本视频在线播放| 久久久精品欧美日韩精品| 国产色爽女视频免费观看| 深夜a级毛片| 精品欧美国产一区二区三| 少妇熟女aⅴ在线视频| 国产国拍精品亚洲av在线观看| 亚洲av不卡在线观看| 亚洲精品乱码久久久久久按摩| 亚洲婷婷狠狠爱综合网| 成人高潮视频无遮挡免费网站| 国产高清视频在线观看网站| 乱码一卡2卡4卡精品| 国产三级在线视频| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| 久久久精品欧美日韩精品| 亚洲在线观看片| 3wmmmm亚洲av在线观看| 久久久国产成人精品二区| 亚洲欧美精品自产自拍| 一本久久中文字幕| 免费av不卡在线播放| 国产精品女同一区二区软件| 麻豆一二三区av精品| 国产私拍福利视频在线观看| 亚洲欧美精品综合久久99| 久久久久久九九精品二区国产| 午夜福利在线观看吧| 亚洲在久久综合| 哪里可以看免费的av片| 欧美一级a爱片免费观看看| 久久这里只有精品中国| 中出人妻视频一区二区| 亚洲成人精品中文字幕电影| 日本撒尿小便嘘嘘汇集6| 国产一级毛片七仙女欲春2| 欧美zozozo另类| 国产精品人妻久久久久久| 五月玫瑰六月丁香| 国产又黄又爽又无遮挡在线| 国产精品av视频在线免费观看| 久久久国产成人精品二区| 国产乱人视频| 久99久视频精品免费| 精品不卡国产一区二区三区| 毛片一级片免费看久久久久| 国内精品一区二区在线观看| 国产黄色小视频在线观看| 久久久久久久久中文| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产欧洲综合997久久,| 99久久成人亚洲精品观看| 亚洲国产精品sss在线观看| 麻豆成人午夜福利视频| 寂寞人妻少妇视频99o| 国产高清不卡午夜福利| 亚洲五月天丁香| 国产免费男女视频| 午夜亚洲福利在线播放| 特级一级黄色大片| 美女xxoo啪啪120秒动态图| 大又大粗又爽又黄少妇毛片口| 久久精品国产亚洲av香蕉五月| 成人鲁丝片一二三区免费| 人妻夜夜爽99麻豆av| 久久久久九九精品影院| 亚洲人成网站在线播放欧美日韩| 午夜福利在线观看吧| 少妇高潮的动态图| 99国产精品一区二区蜜桃av| 成人高潮视频无遮挡免费网站| 国产精品精品国产色婷婷| 99热精品在线国产| 亚洲精品粉嫩美女一区| 欧美潮喷喷水| 床上黄色一级片| 国产成人精品一,二区 | 日韩在线高清观看一区二区三区| 国国产精品蜜臀av免费| 国产精品永久免费网站| 亚洲欧洲国产日韩| 国内久久婷婷六月综合欲色啪| 看非洲黑人一级黄片| 干丝袜人妻中文字幕| 日本五十路高清| 日韩 亚洲 欧美在线| 我的老师免费观看完整版| 插阴视频在线观看视频| 中文亚洲av片在线观看爽| 一区二区三区免费毛片| 九九在线视频观看精品| 国模一区二区三区四区视频| 国产老妇女一区| 亚洲婷婷狠狠爱综合网| 日本成人三级电影网站| 午夜免费激情av| 在线观看免费视频日本深夜| 一本久久精品| 亚洲电影在线观看av| 插阴视频在线观看视频| 日本-黄色视频高清免费观看| 一区二区三区免费毛片| 国产精品1区2区在线观看.| 成人午夜高清在线视频| 少妇高潮的动态图| 精品久久久久久久末码| 夫妻性生交免费视频一级片| 中国美女看黄片| 婷婷精品国产亚洲av| 一个人观看的视频www高清免费观看| 黄色配什么色好看| 少妇猛男粗大的猛烈进出视频 | 久久久色成人| 成人毛片a级毛片在线播放| 99在线视频只有这里精品首页| 你懂的网址亚洲精品在线观看 | 又黄又爽又刺激的免费视频.| 久久中文看片网| 国产大屁股一区二区在线视频| 日日撸夜夜添| 国产日韩欧美在线精品| 精品久久国产蜜桃| 成人一区二区视频在线观看| 欧美性猛交黑人性爽| 天堂网av新在线| 午夜福利在线在线| 欧美色视频一区免费| 亚洲欧美精品自产自拍| 六月丁香七月| 久久久久久久久久成人| 一级毛片久久久久久久久女| www.av在线官网国产| 亚洲av第一区精品v没综合| 一个人看视频在线观看www免费| 嫩草影院新地址| 内射极品少妇av片p| 99九九线精品视频在线观看视频| 一个人看视频在线观看www免费| 一个人免费在线观看电影| 天天一区二区日本电影三级| 卡戴珊不雅视频在线播放| 黄色一级大片看看| 亚洲精品影视一区二区三区av| 在线观看66精品国产| 18禁裸乳无遮挡免费网站照片| 日本成人三级电影网站| 老司机影院成人| 12—13女人毛片做爰片一| 乱系列少妇在线播放| 日本免费一区二区三区高清不卡| a级毛色黄片| 亚洲乱码一区二区免费版| 男人狂女人下面高潮的视频| 亚洲欧美成人综合另类久久久 | 亚洲成人av在线免费| 亚洲精品久久国产高清桃花| 久久精品91蜜桃| 最近手机中文字幕大全| av在线蜜桃| 国产黄色视频一区二区在线观看 | 乱人视频在线观看| 天天一区二区日本电影三级| 欧美xxxx黑人xx丫x性爽| 国产精品久久久久久久电影| 亚洲国产色片| 日日啪夜夜撸| 精品日产1卡2卡| 日韩精品有码人妻一区| 男女那种视频在线观看| av又黄又爽大尺度在线免费看 | 亚洲欧美日韩东京热| 色播亚洲综合网| 国产高清不卡午夜福利| 久久久久网色| 神马国产精品三级电影在线观看| 免费看av在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 日本三级黄在线观看| 国产精品国产三级国产av玫瑰| 黄色欧美视频在线观看| 久久久精品欧美日韩精品| 久久精品国产亚洲av涩爱 | 国产一级毛片七仙女欲春2| 国产人妻一区二区三区在| 天天躁夜夜躁狠狠久久av| 日韩一区二区三区影片| 五月伊人婷婷丁香| 天堂影院成人在线观看| 色哟哟·www| 人人妻人人看人人澡| av在线天堂中文字幕| 亚洲最大成人中文| 一本久久精品| 久久精品人妻少妇| 久久久国产成人精品二区| 亚洲欧洲国产日韩| 午夜精品一区二区三区免费看| 欧美一区二区国产精品久久精品| 日韩国内少妇激情av| 久久九九热精品免费| av福利片在线观看| 国产亚洲精品久久久com| 国内少妇人妻偷人精品xxx网站| 99热全是精品| 午夜a级毛片| 日本一本二区三区精品| 噜噜噜噜噜久久久久久91| 国产免费男女视频| 一级毛片我不卡| 一进一出抽搐gif免费好疼| 国产精品一二三区在线看| 免费观看a级毛片全部| 又爽又黄a免费视频| 一边亲一边摸免费视频| 国产免费男女视频| 一本久久精品| 亚洲精品乱码久久久久久按摩| 国产一级毛片在线| 国产精品免费一区二区三区在线| 欧美变态另类bdsm刘玥| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久久电影| 99热6这里只有精品| 日本色播在线视频| 国产色爽女视频免费观看| 久久久精品94久久精品| 国产精品一区二区三区四区免费观看| 日本欧美国产在线视频| 久久99蜜桃精品久久| 99热这里只有精品一区| 国产亚洲91精品色在线| 欧美人与善性xxx| 97热精品久久久久久| 国内少妇人妻偷人精品xxx网站| 久久久国产成人精品二区| 日韩欧美一区二区三区在线观看| 啦啦啦韩国在线观看视频| 免费看美女性在线毛片视频| 人体艺术视频欧美日本| 国产麻豆成人av免费视频| 婷婷色综合大香蕉| 99热这里只有是精品50| 欧美一区二区亚洲| 国产熟女欧美一区二区| 美女 人体艺术 gogo| 国产不卡一卡二| 亚洲在线自拍视频| 亚洲成人av在线免费| 国产成人aa在线观看| eeuss影院久久| 亚洲av成人精品一区久久| or卡值多少钱| 三级男女做爰猛烈吃奶摸视频| 能在线免费观看的黄片| 成人三级黄色视频| 国产成人a区在线观看| 欧美性感艳星| 小说图片视频综合网站| 99热全是精品| 国产成人91sexporn| 国产欧美日韩精品一区二区| 嫩草影院新地址| 亚洲第一区二区三区不卡| a级毛色黄片| 国产成人91sexporn| 亚洲人成网站在线观看播放| 日韩av在线大香蕉| 22中文网久久字幕| 蜜桃久久精品国产亚洲av| 国产成人福利小说| 黄色配什么色好看| 在线免费十八禁| av国产免费在线观看| 精品久久久久久久久久免费视频| 久久综合国产亚洲精品| 欧美三级亚洲精品| 干丝袜人妻中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av天堂在线播放| 亚洲精品乱码久久久v下载方式| 国产av不卡久久| 久久精品国产亚洲av涩爱 | 久久99蜜桃精品久久| 乱码一卡2卡4卡精品| 我要看日韩黄色一级片| 国产亚洲精品久久久久久毛片| 国产真实乱freesex| 黄片wwwwww| 午夜精品一区二区三区免费看| 一夜夜www| 三级国产精品欧美在线观看| 欧美成人免费av一区二区三区| 国产成年人精品一区二区| 中国美白少妇内射xxxbb| videossex国产| 国产极品天堂在线| 国产av一区在线观看免费| 噜噜噜噜噜久久久久久91| 小说图片视频综合网站| 桃色一区二区三区在线观看| www.色视频.com| 国产欧美日韩精品一区二区| 毛片一级片免费看久久久久| 亚洲av中文av极速乱| 日本在线视频免费播放| 最近视频中文字幕2019在线8| 国产亚洲精品久久久久久毛片| 日本黄色片子视频| 亚洲国产欧美在线一区| 日韩国内少妇激情av| 麻豆成人午夜福利视频| 久久精品夜夜夜夜夜久久蜜豆| 波多野结衣巨乳人妻| videossex国产| 人人妻人人澡欧美一区二区| 99热这里只有是精品50| 成人午夜高清在线视频| 日韩av不卡免费在线播放| 99久久成人亚洲精品观看| 欧美性感艳星| 精品一区二区三区人妻视频| 国产91av在线免费观看| 亚洲无线在线观看| 你懂的网址亚洲精品在线观看 | 国产毛片a区久久久久| 日韩 亚洲 欧美在线| 色综合亚洲欧美另类图片| 国产午夜精品论理片| 老司机影院成人| 久久亚洲精品不卡| 人人妻人人看人人澡| 午夜免费激情av| 亚洲精品成人久久久久久| a级一级毛片免费在线观看| 一边摸一边抽搐一进一小说| 国产精品无大码| 国产精品电影一区二区三区| 一区二区三区免费毛片| 中文欧美无线码| 国产日韩欧美在线精品| 免费人成视频x8x8入口观看| 亚洲精品日韩在线中文字幕 | 成人漫画全彩无遮挡| 寂寞人妻少妇视频99o| 免费黄网站久久成人精品| 天堂√8在线中文| 人人妻人人澡人人爽人人夜夜 | 国产黄片视频在线免费观看| 午夜福利高清视频| 乱码一卡2卡4卡精品| 狠狠狠狠99中文字幕| 亚洲激情五月婷婷啪啪| 一本精品99久久精品77| 69人妻影院| 久久热精品热| 国产 一区 欧美 日韩| 精品免费久久久久久久清纯| 中文字幕av成人在线电影| 久久久a久久爽久久v久久| 女人十人毛片免费观看3o分钟| 欧美又色又爽又黄视频| 成人美女网站在线观看视频| 69av精品久久久久久| 此物有八面人人有两片| 亚洲在线观看片| 亚洲精品影视一区二区三区av| 91aial.com中文字幕在线观看| 成人一区二区视频在线观看| 日韩欧美精品免费久久|