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

    2020年江蘇主汛期短期暴雨預(yù)報檢驗*

    2022-04-07 10:52:04康志明
    氣象 2022年3期
    關(guān)鍵詞:預(yù)報員漏報對流

    蘇 翔 劉 梅 康志明 李 昕

    1 中國氣象局交通氣象重點開放實驗室,南京 210041 2 江蘇省氣象臺,南京 210041 3 南京氣象科技創(chuàng)新研究院,南京 210041

    提 要: 基于江蘇預(yù)報業(yè)務(wù)常用的三個全球模式(ECMWF、NCEP-GFS、CMA-GFS)、三個區(qū)域模式(CMA-MESO、CMA-SH9、PWAFS)、本地客觀預(yù)報和預(yù)報員主觀預(yù)報,對2020年江蘇主汛期(6—9月)中的短期暴雨預(yù)報總體性能進行檢驗,并按降水性質(zhì)分穩(wěn)定性暴雨和對流性暴雨分別進行了檢驗和個例展示。結(jié)果表明:從總體檢驗來看,預(yù)報員的暴雨主觀預(yù)報僅在24 h時效內(nèi)相比于ECMWF有TS評分正技巧,本地客觀和預(yù)報員暴雨預(yù)報的TS評分>區(qū)域模式>除ECMWF以外的全球模式,且無論是全球模式還是區(qū)域模式,并非分辨率越高預(yù)報性能越好;區(qū)域模式的暴雨預(yù)報大多面積偏大,空報率較高;全球模式除ECMWF以外的暴雨預(yù)報大多面積偏小,漏報率較高;20時起報的CMA-MESO的預(yù)報性能明顯好于08時起報的預(yù)報;對于穩(wěn)定性暴雨預(yù)報,應(yīng)重點參考ECMWF、48 h時效內(nèi)的本地客觀預(yù)報和20時起報的CMA-MESO;對于對流性暴雨預(yù)報,應(yīng)重點參考區(qū)域模式CMA-SH9、PWAFS和24 h時效內(nèi)的本地客觀預(yù)報。

    引 言

    2020年6月9日江蘇淮河以南地區(qū)自南向北先后入梅,江蘇正式進入主汛期(6—9月)。期間經(jīng)歷了長達51 d的超長梅雨期(6月9日至7月29日),暴雨頻發(fā),全省平均梅雨量達615.4 mm,是常年同期(236 mm)的2.61倍,為有氣象記錄以來第二多(僅次于1991年的738.8 mm)。不少氣象工作者通過診斷分析(杜小玲等,2016;蔣駿等,2016;金瓊等,2020;毛程燕等,2019)和數(shù)值試驗(劉建勇等,2011;閔錦忠等,2018;梅疏影和閔錦忠,2018;張舒陽和閔錦忠,2018)研究了梅雨鋒暴雨的環(huán)流特征和發(fā)展機制,而且在實際暴雨預(yù)報業(yè)務(wù)中,數(shù)值模式提供的降水預(yù)報是預(yù)報員最重要的參考依據(jù)。

    一些學(xué)者對部分業(yè)務(wù)模式的暴雨預(yù)報性能進行了檢驗評估。孫素琴等(2017)基于歐洲中期天氣預(yù)報中心(ECMWF)、日本氣象廳(JMA)和T639模式降水預(yù)報產(chǎn)品,對2015年江西省汛期區(qū)域性暴雨預(yù)報進行了檢驗,發(fā)現(xiàn)ECMWF模式的降水預(yù)報分布總體偏北,三個模式都存在暴雨漏報的問題。黨英娜(2018)基于ECMWF和華東區(qū)域氣象中心模式(CMA-SH9),對山東半島2016—2017年汛期暴雨預(yù)報能力進行了檢驗,發(fā)現(xiàn)ECMWF暴雨預(yù)報偏弱,CMA-SH9空報率高。宮宇等(2018)對比了中國氣象局全球同化預(yù)報系統(tǒng)(CMA-GFS)、ECMWF和T639在2013—2015年的20次南方暴雨、6次北方暴雨以及12次臺風(fēng)暴雨預(yù)報中的性能,發(fā)現(xiàn)CMA-GFS模式明顯優(yōu)于T639,預(yù)報穩(wěn)定性較好。然而上述研究與江蘇預(yù)報業(yè)務(wù)的實際情況存在一定差距,一方面是區(qū)域氣候特征上的差異,另一方面是所用的業(yè)務(wù)模式有所不同,例如JMA和T639模式在江蘇已被其他的數(shù)值模式所取代。目前江蘇預(yù)報業(yè)務(wù)上常用的三個全球模式為ECMWF、美國國家環(huán)境預(yù)報中心全球預(yù)報系統(tǒng)(NCEP-GFS)、CMA-GFS,常用的三個區(qū)域模式為中國氣象局中尺度天氣數(shù)值預(yù)報系統(tǒng)(CMA-MESO)、CMA-SH9和江蘇本地研發(fā)的精確化天氣分析與預(yù)報系統(tǒng)(PWAFS)(Li et al,2016;2020),與之相關(guān)的本地化檢驗評估研究仍相對缺乏。

    此外,基于統(tǒng)計學(xué)的客觀訂正和基于預(yù)報員經(jīng)驗的主觀訂正相比于數(shù)值模式是否有預(yù)報技巧也是預(yù)報業(yè)務(wù)亟需了解的問題。為了改進預(yù)報水平,除了大力發(fā)展更加先進的數(shù)值模式以外,模式統(tǒng)計后處理技術(shù)也在暴雨預(yù)報中得到了廣泛應(yīng)用(代刊等,2018;蘇翔和袁慧玲,2020)。例如,江蘇采用改進的頻率匹配法(Zhu and Luo,2015;蘇翔等,2021),對ECMWF確定性預(yù)報進行訂正,生成本地客觀降水預(yù)報業(yè)務(wù)產(chǎn)品。2020年主汛期,江蘇還引進中央氣象臺主客觀融合算法(唐健等,2018),并在本地化部署時開發(fā)了集合融合和物理量融合等多種落區(qū)與格點預(yù)報融合技術(shù),讓預(yù)報員在網(wǎng)格預(yù)報的基礎(chǔ)上制作編輯落區(qū)預(yù)報,實現(xiàn)落區(qū)反演格點預(yù)報產(chǎn)品入庫,形成統(tǒng)一的主客觀融合產(chǎn)品,解決了傳統(tǒng)預(yù)報“兩張皮”的問題,即預(yù)報員的降水落區(qū)預(yù)報與上傳的站點考核預(yù)報不一致的問題。

    除了需要考察江蘇2020年主汛期暴雨預(yù)報的整體性能以外,由于不同天氣系統(tǒng)引發(fā)的暴雨存在物理機制上的差異(傅云飛等,2020),可進一步將暴雨樣本按降水性質(zhì)分為穩(wěn)定性降水和對流性降水分別進行檢驗。例如梅雨鋒中不同尺度天氣系統(tǒng)的相互作用使得梅雨鋒暴雨的形成機制較為復(fù)雜(陸漢城,2019),既有層狀云降水,也可產(chǎn)生較強的積云對流(壽紹文,2019)。按降水性質(zhì)分類統(tǒng)計不同預(yù)報產(chǎn)品的暴雨預(yù)報性能對于業(yè)務(wù)具有更加實際的指導(dǎo)意義。

    本文基于江蘇預(yù)報業(yè)務(wù)常用的數(shù)值模式、本地客觀、預(yù)報員主觀的短期降水預(yù)報產(chǎn)品和站點觀測資料,首先對比不同預(yù)報產(chǎn)品在2020年江蘇主汛期短期暴雨預(yù)報中的整體性能,其次檢驗不同預(yù)報產(chǎn)品在穩(wěn)定性暴雨和對流性暴雨預(yù)報中的不同表現(xiàn),最后選取典型暴雨個例進行展示說明。這些分析有助于預(yù)報員了解不同預(yù)報產(chǎn)品在江蘇主汛期暴雨預(yù)報中的性能,總結(jié)預(yù)報經(jīng)驗,以期提升江蘇主汛期暴雨預(yù)報水平。

    1 資 料

    1.1 觀測資料

    所用的觀測資料來源于中國氣象局考核的259個江蘇基本站和優(yōu)質(zhì)加密站的08時至次日08時和20時至次日20時(北京時,下同)累積降水。相比于僅使用基本站和使用所有加密站,可在保證較高數(shù)據(jù)質(zhì)量的同時提供更加豐富的降水細(xì)節(jié)信息。需要注意的是,使用不同的觀測資料進行檢驗所得的預(yù)報評分會有差異。

    1.2 模式資料

    江蘇預(yù)報業(yè)務(wù)常用的三個全球模式分別為:(1)歐洲中期天氣預(yù)報中心ECMWF模式,水平分辨率為0.125°,每日進行2次0~15 d預(yù)報;(2)美國國家環(huán)境預(yù)報中心NCEP-GFS模式,水平分辨率為0.5°,每日進行4次0~16 d預(yù)報;(3)我國自主研發(fā)的CMA-GFS模式,水平分辨率為0.25°,每日進行2次0~10 d預(yù)報。

    江蘇預(yù)報業(yè)務(wù)常用的三個區(qū)域模式分別為:(1)我國自主研發(fā)的中尺度模式CMA-MESO,水平分辨率為3 km,對流可解析(無積云對流參數(shù)化),每日進行2次0~36 h預(yù)報,采用NCEP同時刻分析場和云分析生成模式初值(許晨璐等,2017;張小雯等,2020);(2)華東區(qū)域氣象中心模式CMA-SH9,水平分辨率為9 km,每日進行4次0~72 h預(yù)報,同化了中國高低空常規(guī)觀測資料、雷達基數(shù)據(jù)以及風(fēng)云衛(wèi)星等多源資料,采用復(fù)雜云分析技術(shù)改善了模式初始場,并對模式的對流參數(shù)化方案和微物理過程進行了優(yōu)化(呂林宜等,2019);(3)江蘇本地PWAFS模式,由江蘇省氣象局與美國俄克拉何馬大學(xué)CAPS實驗室合作搭建,內(nèi)核使用ARW-WRF V3.5.1,采用單向嵌套雙重網(wǎng)格,外層區(qū)域水平分辨率為15 km,采用Kain-Fritsch積云對流方案,內(nèi)層區(qū)域水平分辨率為3 km(本文使用3 km分辨率),對流可解析,每日進行8次預(yù)報,其中08時和20時起報的預(yù)報時效均為0~72 h(Li et al,2016)。

    1.3 本地客觀預(yù)報

    目前江蘇預(yù)報業(yè)務(wù)上采用改進的頻率匹配法(Zhu and Luo,2015;蘇翔等,2021),對ECMWF確定性降水預(yù)報進行實時訂正,生成與原始預(yù)報具有相同水平分辨率的本地客觀預(yù)報業(yè)務(wù)產(chǎn)品。對不同起報時間和預(yù)報時效單獨進行滑動建模,所用的訓(xùn)練樣本為去年同期前后各一個月(共61 d)的歷史預(yù)報和觀測資料。本地客觀預(yù)報算法的目的是通過消除ECMWF模式的頻率偏差對降水強度進行訂正,其降水分布形態(tài)通常與ECMWF模式預(yù)報具有較高的相似性??紤]到ECMWF模式數(shù)據(jù)計算和傳輸?shù)臅r間滯后性,本地客觀預(yù)報的起報時間與其所訂正的ECMWF模式降水預(yù)報產(chǎn)品相比落后12 h。例如本地客觀08時起報的24、48、72 h時效產(chǎn)品所訂正的對象為ECMWF模式前一日20時起報的36、60、84 h時效產(chǎn)品。

    1.4 預(yù)報員主觀預(yù)報

    江蘇2020年主汛期開始實行基于主客觀融合的預(yù)報業(yè)務(wù)流程。預(yù)報員自由選取不同的數(shù)值模式和客觀預(yù)報算法作為網(wǎng)格預(yù)報背景場,并在此基礎(chǔ)上制作編輯降水落區(qū)預(yù)報。系統(tǒng)平臺會通過主客觀融合算法將預(yù)報員繪制的降水落區(qū)預(yù)報轉(zhuǎn)換為連續(xù)變化的格點預(yù)報,形成統(tǒng)一的主客觀融合產(chǎn)品庫。每日06時之前由值班預(yù)報員制作完成08時起報的24、48和72 h逐日降水預(yù)報產(chǎn)品,每日17時之前由首席預(yù)報員制作完成20時起報的24、48和72 h逐日降水預(yù)報產(chǎn)品。

    2 檢驗方法

    2.1 主汛期暴雨樣本的選取

    從2020年江蘇主汛期(6—9月)08時至次日08時和20時至次日20時觀測出現(xiàn)暴雨的日期中剔除暴雨站點數(shù)小于總站點數(shù)2%的局地暴雨日期,分別得到35個和39個暴雨觀測樣本用于江蘇主汛期暴雨預(yù)報的總體性能檢驗。需要注意的是,多樣本總體檢驗可能存在部分降水較多日樣本評分覆蓋降水較少日樣本評分的情況,檢驗結(jié)果僅代表統(tǒng)計時段內(nèi)的總體性能,不代表單日樣本性能。

    2.2 預(yù)報時效的對應(yīng)

    模式資料由于計算和傳輸?shù)仍?,實際的產(chǎn)品可用時間滯后于模式起報時間。例如預(yù)報員早上制作08時起報的24 h預(yù)報時,實際可參考的資料為同一起報時間和預(yù)報時效的本地客觀產(chǎn)品以及各個模式前一日20時起報的36 h時效產(chǎn)品。為了考察預(yù)報員的主觀訂正相比于模式和客觀預(yù)報產(chǎn)品是否有預(yù)報技巧,將預(yù)報員和本地客觀08時起報的24、48、72 h時效產(chǎn)品與各個模式前一日20時起報的36、60、84 h時效產(chǎn)品進行對比檢驗,并將預(yù)報員和本地客觀20時起報的24、48、72 h時效產(chǎn)品與各個模式當(dāng)日08時起報的36、60、84 h時效產(chǎn)品進行對比檢驗。根據(jù)不同模式預(yù)報時效的長度,CMA-MESO模式僅檢驗36 h時效,CMA-SH9和PWAFS模式僅檢驗36 h和60 h時效。

    2.3 穩(wěn)定性暴雨和對流性暴雨的劃分

    對于24 h累積降水,目前并沒有一個明確的方法可以劃分穩(wěn)定性降水和對流性降水。層云穩(wěn)定性降水和積云對流性降水既可能存在水平和垂直空間上的疊加,也可能存在時間上更替,可根據(jù)兩者出現(xiàn)的比例,按一定客觀標(biāo)準(zhǔn)劃分為以穩(wěn)定性降水為主的暴雨樣本和以對流性降水為主的暴雨樣本??紤]到穩(wěn)定性降水具有持續(xù)時間長的特點,對流性降水具有強度大的特點,可采用逐小時站點降水觀測數(shù)據(jù)進行客觀劃分。根據(jù)經(jīng)驗和測試,設(shè)計了如下劃分標(biāo)準(zhǔn):計算將暴雨(≥50 mm·d-1)站點上的24 h累積降水拆分為逐小時降水,計算其中短時強降水(20 mm·h-1)小時數(shù)與降水(>0 mm·h-1)小時數(shù)的比值r;若r≥5%,則判斷為對流性降水(為主);若r<5%,則判斷為穩(wěn)定性降水(為主)。按此標(biāo)準(zhǔn),08時至次日08時和20時至次日20時暴雨觀測樣本可分別劃分為9個和12個穩(wěn)定性暴雨樣本,26個和27個對流性暴雨樣本。

    2.4 主汛期暴雨的總體檢驗評分

    先將所有模式、本地客觀和預(yù)報員主觀網(wǎng)格預(yù)報產(chǎn)品統(tǒng)一雙線性插值到觀測站點上,再分別計算暴雨(≥50 mm·d-1)的頻率偏差(FB)、TS評分、漏報率(PO)、空報率(FAR)(Jolliffe and Stephenson,2003):

    (1)

    (2)

    (3)

    (4)

    式中:H、M、F分別為暴雨命中樣本數(shù)、漏報樣本數(shù)和空報樣本數(shù)。這里的樣本數(shù)是時間樣本數(shù)和空間樣本(站點)數(shù)的乘積。式(1)中的頻率偏差為統(tǒng)計時段內(nèi)預(yù)報暴雨總樣本數(shù)與觀測暴雨總樣本數(shù)的比值,越接近1代表預(yù)報的系統(tǒng)偏差越小,大于1代表暴雨預(yù)報面積偏大,小于1代表暴雨預(yù)報面積偏小。此外,頻率偏差和TS評分還可以分別寫成漏報率和空報率的函數(shù):

    (5)

    (6)

    由式(5)可知,當(dāng)漏報率高于空報率時,暴雨預(yù)報面積偏??;當(dāng)空報率高于漏報率時,暴雨預(yù)報面積偏大。由式(6)可知,TS評分的數(shù)值由漏報率和空報率共同決定。選取ECMWF模式為參考模式,將不同預(yù)報產(chǎn)品的TS評分與之進行比較,若TS評分高于ECMWF,則為TS評分正技巧,反之則為TS評分負(fù)技巧。

    2.5 小樣本檢驗的注意事項

    小樣本檢驗可能存在評分不穩(wěn)定的情況。在大樣本總體檢驗時,通常默認(rèn)式(1)~式(4)的分母較大,不會出現(xiàn)評分不穩(wěn)定的情況;但在小樣本檢驗時,例如個例檢驗或按降水性質(zhì)分類檢驗時,式(1)~式(4)的分母可能會變得很小,導(dǎo)致評分?jǐn)?shù)值不穩(wěn)定,引起錯誤理解。例如,NCEP-GFS和CMA-GFS模式2020年6月17日08時起報的12~36 h累積降水相比于2020年6月17日20時至18日20時的降水觀測都沒有報出沿淮分布的暴雨雨帶,兩個模式預(yù)報的暴雨雨帶幾乎都在江蘇區(qū)域以外(3.3節(jié)中將對此進行詳細(xì)討論)。兩個模式預(yù)報的暴雨雨帶僅存在細(xì)微的差異,暴雨命中站點數(shù)H分別為1和0,暴雨空報站點數(shù)F分別為0和2,都遠(yuǎn)低于暴雨漏報站點數(shù)37和38,但卻導(dǎo)致兩者的空報率FAR的數(shù)值出現(xiàn)了0和1的巨大差別,與預(yù)報員的主觀感受偏差較大,這是由于暴雨預(yù)報樣本數(shù)(H+F的值)太小導(dǎo)致的評分?jǐn)?shù)值不穩(wěn)定??紤]到文中的個例檢驗或按降水性質(zhì)分類檢驗所用的暴雨樣本都是基于觀測選取的,即H+M的值不會太小,因此在上述小樣本檢驗時改用空報面積比(FAAR)代替空報率(FAR)以獲得更加符合預(yù)報員主觀感受的檢驗評分結(jié)果:

    (7)

    即:空報面積比是暴雨空報樣本數(shù)F與暴雨觀測樣本數(shù)H+M的比值,與漏報率PO具有相同的分母,是空報率FAR與頻率偏差FB的乘積。上述例子中兩個模式的空報面積比FAAR分別為0和0.05,代表模式空報的站點數(shù)相對于暴雨觀測站點數(shù)都很低,更加符合預(yù)報員的主觀認(rèn)知。在本文的穩(wěn)定性暴雨和對流性暴雨檢驗中將使用空報面積比代替空報率。

    3 檢驗結(jié)果

    3.1 總體性能檢驗

    3.1.1 24 h時效暴雨預(yù)報檢驗

    圖1為2020年江蘇主汛期預(yù)報員和本地客觀08時起報的24 h(對應(yīng)模式20時起報的36 h)暴雨預(yù)報檢驗評分,并已按評分?jǐn)?shù)值大小對預(yù)報產(chǎn)品進行排序。全球模式NCEP-GFS和CMA-GFS的暴雨預(yù)報面積最小且暴雨漏報率最高,區(qū)域模式CMA-MESO和CMA-SH9的暴雨預(yù)報面積最大且暴雨空報率較高。從暴雨預(yù)報的TS評分上看,本地客觀、預(yù)報員和CMA-MESO相比于ECMWF都具有一定正技巧。本地客觀和預(yù)報員較好地訂正了ECMWF模式暴雨預(yù)報面積偏小的系統(tǒng)偏差,訂正以后的頻率偏差值更接近1,兩者暴雨預(yù)報的TS評分也最高。CMA-MESO暴雨預(yù)報的TS評分正技巧主要得益于較低的暴雨漏報率。區(qū)域模式PWAFS由于暴雨漏報率和空報率均高于ECMWF,因此暴雨TS評分相比于ECMWF為負(fù)技巧。

    圖2為2020年江蘇主汛期預(yù)報員和本地客觀20時起報的24 h(對應(yīng)模式08時起報的36 h)暴雨預(yù)報檢驗評分。部分結(jié)論與圖1一致,即全球模式CMA-GFS和NCEP-GFS的暴雨預(yù)報面積最小且具有最大的暴雨漏報率,區(qū)域模式CMA-SH9和CMA-MESO的暴雨預(yù)報面積最大且具有最大的暴雨空報率。而在暴雨TS評分方面,預(yù)報員和本地客觀預(yù)報相比于ECMWF具有正技巧,這主要得益于較低的暴雨漏報率。區(qū)域模式PWAFS的暴雨漏報率略低于ECMWF,但空報率顯著高于ECMWF,因此其相比于ECMWF模式的暴雨TS評分為負(fù)技巧。值得注意的是,20時和08時起報的CMA-MESO的暴雨TS評分差異較大。雖然兩個時次起報的CMA-MESO的暴雨空報率均較高,但前者的暴雨漏報率低,因此暴雨TS評分較高(圖1),后者的暴雨漏報率與ECMWF相當(dāng),暴雨TS評分較低(圖2)。

    3.1.2 48 h時效暴雨預(yù)報檢驗

    圖3為2020年江蘇主汛期預(yù)報員和本地客觀08時起報的48 h(對應(yīng)模式20時起報的60 h)暴雨預(yù)報檢驗評分。全球模式NCEP-GFS和CMA-GFS的暴雨預(yù)報面積最小且暴雨漏報率最高。僅有區(qū)域模式PWAFS和CMA-SH9的暴雨預(yù)報面積輕微偏大,其他產(chǎn)品的暴雨預(yù)報面積均偏小。從暴雨TS評分上看,所有預(yù)報產(chǎn)品相比于ECMWF均為負(fù)技巧。雖然區(qū)域模式PWAFS和CMA-SH9的暴雨漏報率略低于ECMWF,但暴雨空報率較高,因此暴雨TS評分相比于ECMWF為負(fù)技巧。本地客觀預(yù)報和NCEP-GFS雖然空報率低于ECMWF,但暴雨漏報率較高,因此其暴雨TS評分也為負(fù)技巧。

    圖4為2020年江蘇主汛期預(yù)報員和本地客觀20時起報的48 h(對應(yīng)模式08時起報的60 h)暴雨預(yù)報檢驗評分。部分結(jié)論與圖3一致,即全球模式CMA-GFS和NCEP-GFS的暴雨預(yù)報面積最小且具有最大的暴雨漏報率,區(qū)域模式PWAFS和CMA-SH9的暴雨預(yù)報面積最大且具有最大的暴雨空報率。所有預(yù)報產(chǎn)品相比于ECMWF的暴雨預(yù)報TS評分均為負(fù)技巧,其中區(qū)域模式CMA-SH9和PWAFS的暴雨預(yù)報TS評分負(fù)技巧主要源于較高的空報率,全球模式NCEP-GFS和CMA-GFS的暴雨預(yù)報TS評分負(fù)技巧主要來源于較高的漏報率。

    3.1.3 72 h時效暴雨預(yù)報檢驗

    圖5為2020年江蘇主汛期預(yù)報員和本地客觀08時起報的72 h(對應(yīng)模式20時起報的84 h)暴雨預(yù)報檢驗評分。全球模式CMA-GFS的暴雨TS評分最差,原因是其暴雨漏報率和暴雨空報率都是最高的。NCEP-GFS暴雨預(yù)報面積最小,暴雨漏報率很高,但暴雨空報率最低,因此暴雨TS評分略高于CMA-GFS。預(yù)報員和本地客觀的暴雨漏報率和暴雨空報率均略高于ECMWF,因此暴雨TS評分均略低于ECMWF。

    圖1 2020年江蘇主汛期預(yù)報員和本地客觀08時起報的24 h(對應(yīng)模式20時起報的36 h)暴雨預(yù)報的 (a)TS評分,(b)頻率偏差,(c)漏報率和(d)空報率Fig.1 The 24 h forecasts of forecasters and local objective method initialized at 08:00 BT (corresponding to the 36 h forecast of models initialized at 20:00 BT) in terms of torrential rain during the 2020 Jiangsu main flood season (a) TS, (b) frequency bias, (c) miss rate, (d) false alarm ratio

    圖2 同圖1,但為預(yù)報員和本地客觀20時起報 (對應(yīng)模式08時起報)Fig.2 Same as Fig.1, but for forecast of forecasters and local objective method initialized at 20:00 BT (corresponding to the forecast of models initialized at 08:00 BT)

    圖3 2020年江蘇主汛期預(yù)報員和本地客觀08時起報的48 h(對應(yīng)模式20時起報的60 h)暴雨預(yù)報的 (a)TS評分,(b)頻率偏差,(c)漏報率和(d)空報率Fig.3 The 48 h forecast of forecasters and local objective method initialized at 08:00 BT (corresponding to the 60 h forecast of models initialized at 20:00 BT) in terms of torrential rain during the 2020 Jiangsu main flood season (a) TS, (b) frequency bias, (c) miss rate, (d) false alarm ratio

    圖4 同圖3,但為預(yù)報員和本地客觀20時起報 (對應(yīng)模式08時起報)Fig.4 Same as Fig.3, but for forecast of forecasters and local objective method initialized at 20:00 BT (corresponding to the forecast of models initialized at 08:00 BT)

    圖5 2020年江蘇主汛期預(yù)報員和本地客觀08時起報的72 h(對應(yīng)模式20時起報的84 h)暴雨預(yù)報的 (a)TS評分,(b)頻率偏差,(c)漏報率和(d)空報率Fig.5 The 72 h forecast of forecasters and local objective method initialized at 08:00 BT (corresponding to the 84 h forecast of models initialized at 20:00 BT) in terms of torrential rain during the 2020 Jiangsu main flood season (a) TS, (b) frequency bias, (c) miss rate, (d) false alarm ratio

    圖6為2020年江蘇主汛期預(yù)報員和本地客觀20時起報的72 h(對應(yīng)模式08時起報的84 h)暴雨預(yù)報檢驗評分。結(jié)論與圖5基本一致,主要區(qū)別在于08時起報的CMA-GFS暴雨預(yù)報面積偏小的程度相比于20時起報的更大,其暴雨空報率也明顯下降。

    從不同預(yù)報時效的總體檢驗來看:本地客觀預(yù)報和預(yù)報員的暴雨預(yù)報TS評分 > 區(qū)域模式 > 除ECMWF以外的全球模式,僅有24 h時效的本地客觀預(yù)報、預(yù)報員以及20時起報的CMA-MESO相比于ECMWF模式具有暴雨TS評分正技巧,而08時起報的CMA-MESO則性能較差;區(qū)域模式大多暴雨預(yù)報面積偏大,空報率較高;全球模式除ECMWF以外大多暴雨預(yù)報面積偏小,暴雨漏報率較高。

    3.2 按降水性質(zhì)檢驗

    3.2.1 穩(wěn)定性暴雨預(yù)報檢驗

    表1和表2展示了不同起報時間和預(yù)報時效下,不同預(yù)報產(chǎn)品的穩(wěn)定性暴雨TS評分。6個模式里面ECMWF的總體表現(xiàn)最好,僅有20時起報的36 h時效的CMA-MESO和CMA-SH9以及08時起報的84 h時效的NCEP-GFS相比于ECMWF具有暴雨TS評分正技巧。預(yù)報員不同起報時間和預(yù)報時效的暴雨預(yù)報TS評分與ECMWF都相差不大,說明預(yù)報員的主觀預(yù)報參考ECMWF模式較多。然而預(yù)報員僅有08時起報的24 h預(yù)報和20時起報的48 h預(yù)報相比于ECMWF具有暴雨TS評分正技巧,且訂正技巧不如本地客觀預(yù)報。值得注意的是20時起報的36 h時效的CMA-MESO和CMA-SH9相比于ECMWF也具有暴雨TS評分正技巧,這主要得益于兩者具有較小的暴雨漏報率(表略),但也伴隨著較大的暴雨空報面積比(表略)。此外,08時起報的84 h時效的NCEP-GFS暴雨預(yù)報的TS評分也高于ECMWF,這是由于其暴雨空報面積比很低(表略)。總體來看,對于穩(wěn)定性暴雨,ECMWF模式的整體表現(xiàn)較好,48 h時效內(nèi)的本地客觀預(yù)報對ECMWF的訂正大多具有暴雨TS評分正技巧,可供預(yù)報員參考。

    圖6 同圖5,但為預(yù)報員和本地客觀20時起報 (對應(yīng)模式08時起報)Fig.6 Same as Fig.5, but for forecast of forecasters and local objective method initialized at 20:00 BT (corresponding to the forecast of models initialized at 08:00 BT)

    表1 2020年江蘇主汛期預(yù)報員和本地客觀08時起報的24, 48, 72 h (對應(yīng)模式20時起報的36, 60, 84 h) 暴雨預(yù)報中穩(wěn)定性降水樣本的TS評分Table 1 TSs of stable precipitation samples of the 24, 48, 72 h torrential rain forecast of forecasters and local objective method initialized at 08:00 BT during the 2020 Jiangsu main flood season (corresponding to the 36, 60, 84 h torrential rain forecast of models initialized at 20:00 BT)

    表2 同表1,但為預(yù)報員和本地客觀20時起報(對應(yīng)模式08時起報)Table 2 Same as Table 1, but for forecast of forecasters and local objective method initialized at 20:00 BT (corresponding to the forecast of models initialized at 08:00 BT)

    3.2.2 對流性暴雨預(yù)報檢驗

    表3和表4展示了不同起報時間和預(yù)報時效下,不同預(yù)報產(chǎn)品的對流性暴雨TS評分。本地客觀預(yù)報和預(yù)報員僅在24 h時效下表現(xiàn)較好,區(qū)域模式CMA-SH9和PWAFS相比于ECMWF大多具有暴雨TS評分正技巧,主要得益于較低的暴雨漏報率(表略),但也存在暴雨預(yù)報面積明顯偏大的問題。值得注意的是,20時起報的48 h和72 h時效的NCEP-GFS相比于ECMWF具有暴雨TS評分正技巧,這是由于其空報面積比很低(表略),暴雨預(yù)報面積明顯偏小。此外,08時起報的CMA-MESO表現(xiàn)較差的原因在于其空報面積比超過ECMWF的2倍,但暴雨漏報率相比于ECMWF卻沒有降低,虛假暴雨落區(qū)較多??傮w來看,對于對流性暴雨,區(qū)域模式CMA-SH9、PWAFS以及24 h時效的本地客觀預(yù)報總體表現(xiàn)較好,具有參考價值。

    3.3 暴雨個例檢驗

    3.3.1 穩(wěn)定性暴雨個例檢驗

    圖7為不同模式2020年6月17日08時起報的36 h降水預(yù)報、預(yù)報員和本地客觀2020年6月17日20時起報的24 h降水預(yù)報以及對應(yīng)的2020年6月17日20時至18日20時降水觀測,表5為該穩(wěn)定性降水個例中不同預(yù)報產(chǎn)品暴雨等級的不同評分。下面以此個例為例探討不同類型模式在穩(wěn)定性暴雨預(yù)報中的偏差原因。此次暴雨過程中,觀測的暴雨呈帶狀分布在沿淮一帶。而NCEP-GFS和CMA-GFS漏報了這次暴雨過程,與其存在較大的暴雨漏報率且暴雨預(yù)報面積明顯偏小的統(tǒng)計結(jié)論相一致,這可能是由于這兩個全球模式的分辨率較粗,與降水相關(guān)的物理參數(shù)化方案不夠準(zhǔn)確。全球模式ECMWF的分辨率相對較高,對于穩(wěn)定性暴雨雨帶的強度和位置都把握較好,暴雨TS評分最高。本地客觀預(yù)報在該個例中對于ECMWF模式的訂正效果不理想,暴雨落區(qū)有所減小,暴雨TS評分相比于ECMWF為負(fù)技巧。區(qū)域模式由于分辨率較高,均報出了沿淮一帶的暴雨雨帶,但在強度和位置上仍存在一定的差異。PWAFS預(yù)報的暴雨雨帶偏窄且略偏北,CMA-MESO模式預(yù)報的暴雨雨帶偏寬偏南,因此兩者暴雨預(yù)報的TS評分均較低。CMA-SH9模式預(yù)報的暴雨雨帶的位置預(yù)報較好,暴雨TS評分相對較高。這可能是由于PWAFS和CMA-MESO在3 km分辨率上對于對流過程的解析不夠準(zhǔn)確,而分辨率為9 km的CMA-SH9采用優(yōu)化的對流參數(shù)化方案和微物理過程,獲得了較好的預(yù)報效果。預(yù)報員對ECMWF的暴雨進行了增強和落區(qū)擴展,出現(xiàn)了較多的空報,因此暴雨TS評分略低于ECMWF。

    表3 同表1,但為對流性降水Table 3 Same as Table 1, but for convective precipitation samples

    表4 同表2,但為對流性降水Table 4 Same as Table 2, but for convective precipitation samples

    圖7 2020年6月17日08時起報的12~36 h降水預(yù)報:(a)ECMWF、(b)NCEP-GFS、(c)CMA-GFS、 (d)CMA-MESO、(e)CMA-SH9、(f)PWAFS, 17日20時起報的0~24 h降水預(yù)報: (g)本地客觀、(h)預(yù)報員,以及(i)對應(yīng)的17日20時至18日20時降水觀測Fig.7 (a-f) the 12-36 h precipitation forecast initialized at 08:00 BT 17, (g, h) the 0-24 h precipitation forecast initialized at 20:00 BT 17 and (i) the corresponding precipitation observation from 20:00 BT 17 to 20:00 BT 18 June 2020 (a) ECMWF, (b) NCEP-GFS, (c) CMA-GFS, (d) CMA-MESO, (e) CMA-SH9, (f) PWAFS, (g) local objective method, (h) forecaster

    表5 預(yù)報員和本地客觀2020年6月17日20時起報的24 h(對應(yīng)模式2020年6月 17日08時起報的36 h)暴雨預(yù)報的TS評分、頻率偏差、漏報率和空報面積比Table 5 TS, frequency bias, miss rate and false alarm area rate for the 24 h torrential rain forecast of forecasters and local objective method initialized at 20:00 BT 17 June 2020 (corresponding to the 36 h torrential rain forecast initialized at 08:00 BT June 17 2020)

    圖8 同圖7,但為2020年7月18日08時對流性暴雨個例 (a-f)08時起報,(g,h)20時起報,(i)18日20時至19日20時觀測Fig.8 Same as Fig.7, but for the case on 18 July 2020 (a-f) 12-36 h precipitation forecast initialized at 08:00 BT 18, (g, h) 0-24 h precipitation forecast initialized at 20:00 BT 18, (i) observation from 20:00 BT 18 to 20:00 BT 19

    3.3.2 對流性暴雨個例檢驗

    圖8為不同模式2020年7月18日08時起報的36 h降水預(yù)報、預(yù)報員和本地客觀預(yù)報18日20時起報的24 h降水預(yù)報以及對應(yīng)的18日20時至19日20時降水觀測,表6為該對流性降水個例中不同預(yù)報產(chǎn)品暴雨等級的不同評分。以該個例為例探討不同類型模式在對流性暴雨預(yù)報中的偏差原因。此次暴雨過程的雨帶主體分為兩塊,分別位于江淮之間東南部和沿江東部地區(qū)。三個全球模式都沒能報出這次對流性暴雨過程的雨帶主體部分,其中ECMWF的暴雨雨帶位于沿淮和淮北地區(qū),而NCEP-GFS和CMA-GFS暴雨預(yù)報面積明顯偏小,沒有明顯的暴雨雨帶,這可能是由于全球模式分辨率較粗,對于中小尺度天氣系統(tǒng)的預(yù)報能力有限所致。區(qū)域模式均在不同程度上報出了觀測的暴雨雨帶,其中CMA-SH9模式的暴雨落區(qū)與實況最為接近,PWAFS也在一定程度上報出了兩塊暴雨主體,但雨帶較為分散,CMA-MESO模式僅報出了江淮之間東部的暴雨主體。CMA-SH9模式預(yù)報好于PWAFS和CMA-MESO的原因可能也是前者采用了積云對流參數(shù)化方案,雖然預(yù)報的雨帶相對平滑但位置較為準(zhǔn)確,后兩者采用了3 km高分辨率的對流解析取代了積云對流參數(shù)化方案,雖然預(yù)報的雨帶結(jié)構(gòu)較為精細(xì),但暴雨落區(qū)的準(zhǔn)確度不如CMA-SH9模式。本地客觀預(yù)報和預(yù)報員均參考了ECMWF模式預(yù)報的雨帶,落區(qū)偏差較大,暴雨TS評分較低。

    表6 預(yù)報員和本地客觀預(yù)報2020年7月18日20時起報的24 h(對應(yīng)模式2020年7月 18日08時起報的36 h)暴雨預(yù)報的TS評分、頻率偏差、漏報率和空報面積比Table 6 TS, frequency bias, miss rate and false alarm area rate for the 24 h torrential rain forecast of forecasters and local objective method initialized at 20:00 BT 18 July 2020 (corresponding to the 36 h torrential rain forecast initialized at 08:00 BT 18 July 2020)

    4 結(jié)論與討論

    本文對2020年主汛期暴雨預(yù)報中江蘇預(yù)報業(yè)務(wù)常用的三個全球模式(ECMWF、NCEP-GFS、CMA-GFS)、三個區(qū)域模式(CMA-MESO、CMA-SH9、PWAFS)、本地客觀預(yù)報、預(yù)報員主觀預(yù)報進行了總體檢驗評估,并按降水性質(zhì)分穩(wěn)定性暴雨和對流性暴雨樣本分別進行了預(yù)報性能檢驗,最后展示了典型的穩(wěn)定性暴雨和對流性暴雨個例的檢驗分析結(jié)果,得出以下結(jié)論:

    (1)從總體檢驗來看,預(yù)報員的暴雨主觀預(yù)報僅在24 h時效內(nèi)相比于ECMWF有TS評分正技巧,本地客觀預(yù)報和預(yù)報員暴雨預(yù)報的TS評分 > 區(qū)域模式 > 除ECMWF以外的全球模式,且無論是全球模式還是區(qū)域模式,并非分辨率越高預(yù)報性能越好。

    (2)區(qū)域模式的暴雨預(yù)報大多面積偏大,空報率較高;全球模式除ECMWF以外的暴雨預(yù)報大多面積偏小,漏報率較高。

    (3)20時起報的CMA-MESO的預(yù)報性能明顯好于08時起報的預(yù)報。

    (4)對穩(wěn)定性暴雨預(yù)報,應(yīng)重點參考ECMWF、48 h時效內(nèi)的本地客觀預(yù)報和20時起報的區(qū)域模式CMA-MESO。

    (5)對于對流性暴雨預(yù)報,應(yīng)重點參考區(qū)域模式CMA-SH9、PWAFS和24 h時效內(nèi)的本地客觀預(yù)報。

    以上結(jié)論對于江蘇主汛期的實際業(yè)務(wù)預(yù)報具有一定的指示意義。下面就暴雨預(yù)報檢驗中需要注意的幾點問題進行討論:

    (1)本文的研究結(jié)論來源于江蘇區(qū)域2020年主汛期(6—9月)的統(tǒng)計結(jié)果,暴雨樣本數(shù)量有限,尤其是穩(wěn)定性暴雨的數(shù)量相對不足,在未來的工作中還需加入更多年份的暴雨樣本進行檢驗分析。考慮到數(shù)值模式發(fā)展過程中存在不定期的系統(tǒng)升級,模式的預(yù)報偏差性能亦可能發(fā)生變化,還需要跟蹤檢驗其預(yù)報性能的變化。此外,未來還可以考慮將檢驗區(qū)域由江蘇拓展到更大的華東區(qū)域,有助于揭示江蘇范圍模式預(yù)報性能相對于周邊地區(qū)的異同之處,亦可為華東區(qū)域的暴雨預(yù)報業(yè)務(wù)提供參考。

    (2)雖然三個全球模式里面,分辨率較高的ECMWF性能明顯優(yōu)于分辨率相對較低的NCEP-GFS和CMA-GFS,但在三個區(qū)域模式中,分辨率高達3 km 的CMA-MESO和PWAFS總體上不如分辨率9 km的CMA-SH9模式。雖然CMA-MESO和PWAFS采用對流解析取代了積云對流參數(shù)化方案,但并未取得更準(zhǔn)確的預(yù)報效果。此外,20時和08時起報的CMA-MESO模式在預(yù)報性能上存在的較大差異,還需要模式研發(fā)人員根據(jù)模式的設(shè)計研究其中的具體原因。

    (3)基于頻率匹配的本地客觀訂正在江蘇區(qū)域內(nèi)效果有時并不理想。主要原因包括:首先,模型訓(xùn)練期與檢驗期的頻率偏差不可避免地存在一定的差異;其次,江蘇地區(qū)的空間訓(xùn)練樣本數(shù)較小導(dǎo)致模型并不穩(wěn)定;此外,當(dāng)原始模式的暴雨頻率偏差接近1時,頻率匹配的訂正效果也會受到影響。因此,未來還需要進一步加強客觀預(yù)報技術(shù)的研發(fā),增加多模式集成與人工智能等先進的暴雨訂正技術(shù)方法的研究。

    (4)在一些暴雨過程中,穩(wěn)定性降水和對流性降水同時存在,如何科學(xué)合理地對這兩類降水類型進行劃分以及不同的劃分標(biāo)準(zhǔn)是否會對研究結(jié)論產(chǎn)生影響也是需要進一步研究的問題。此外,本文基于觀測的逐小時降水對穩(wěn)定性暴雨和對流性暴雨進行劃分,而在實際業(yè)務(wù)中預(yù)報員還需考慮如何根據(jù)模式預(yù)報資料對暴雨降水性質(zhì)進行預(yù)判以選擇合適的預(yù)報策略。

    (5)在暴雨誤差分類檢驗中,僅僅根據(jù)降水性質(zhì)對暴雨樣本進行劃分還不夠細(xì)致。若對暴雨樣本分類太粗,可能無法準(zhǔn)確刻畫不同天氣背景引發(fā)的暴雨預(yù)報的偏差特征;若分類太細(xì),又可能導(dǎo)致檢驗樣本量不足而失去了統(tǒng)計意義。未來可考慮增加更多年份的暴雨樣本,并根據(jù)不同環(huán)流背景對不同暴雨預(yù)報產(chǎn)品進行分類檢驗。

    猜你喜歡
    預(yù)報員漏報對流
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    江西省天氣預(yù)報員現(xiàn)狀調(diào)查分析與對策建議
    預(yù)報員培訓(xùn)的“復(fù)盤”方法
    航空氣象預(yù)報人為因素的分析初探
    各類氣體報警器防誤報漏報管理系統(tǒng)的應(yīng)用
    基于ANSYS的自然對流換熱系數(shù)計算方法研究
    黔東南州預(yù)報員隊伍建設(shè)的實踐與思考
    二元驅(qū)油水界面Marangoni對流啟動殘余油機理
    傳染病漏報原因分析及對策
    日本廠商在美漏報事故千余起被指管理疏漏
    丝袜美腿诱惑在线| 国产亚洲精品第一综合不卡| 90打野战视频偷拍视频| 国产熟女午夜一区二区三区| 在线播放国产精品三级| 手机成人av网站| 精品亚洲乱码少妇综合久久| 黄色视频不卡| 国产成人影院久久av| 法律面前人人平等表现在哪些方面| 热re99久久国产66热| 激情视频va一区二区三区| 午夜福利一区二区在线看| 性色av乱码一区二区三区2| 老司机午夜十八禁免费视频| 亚洲av成人一区二区三| 国产野战对白在线观看| 99国产综合亚洲精品| 国产成人av教育| 看免费av毛片| 好男人电影高清在线观看| 99re6热这里在线精品视频| 人成视频在线观看免费观看| 狂野欧美激情性xxxx| 亚洲欧美一区二区三区久久| 国产在线观看jvid| 亚洲一码二码三码区别大吗| 午夜激情久久久久久久| 免费在线观看视频国产中文字幕亚洲| 涩涩av久久男人的天堂| 精品久久久久久久毛片微露脸| 99riav亚洲国产免费| 色尼玛亚洲综合影院| av不卡在线播放| 狠狠精品人妻久久久久久综合| 国产91精品成人一区二区三区 | 亚洲精品成人av观看孕妇| 国产日韩欧美在线精品| 欧美精品人与动牲交sv欧美| 首页视频小说图片口味搜索| 下体分泌物呈黄色| 国产xxxxx性猛交| 黄频高清免费视频| 日韩中文字幕视频在线看片| 五月天丁香电影| 99国产极品粉嫩在线观看| 丰满人妻熟妇乱又伦精品不卡| 我要看黄色一级片免费的| 自线自在国产av| 日韩三级视频一区二区三区| 热99re8久久精品国产| 国产精品电影一区二区三区 | 久久精品aⅴ一区二区三区四区| 久久中文字幕人妻熟女| 亚洲国产av影院在线观看| 一区二区三区国产精品乱码| 免费黄频网站在线观看国产| 51午夜福利影视在线观看| 一本—道久久a久久精品蜜桃钙片| av国产精品久久久久影院| 精品国产亚洲在线| 国产亚洲欧美在线一区二区| 亚洲精品乱久久久久久| 五月天丁香电影| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 国产精品.久久久| 最黄视频免费看| 黑人欧美特级aaaaaa片| 国产高清视频在线播放一区| 成在线人永久免费视频| 777米奇影视久久| 丝袜美足系列| 国产日韩欧美在线精品| 桃红色精品国产亚洲av| 50天的宝宝边吃奶边哭怎么回事| 999久久久精品免费观看国产| 欧美老熟妇乱子伦牲交| 丝袜人妻中文字幕| 99精品在免费线老司机午夜| av天堂在线播放| 老司机影院毛片| bbb黄色大片| 中文字幕高清在线视频| 国产在线观看jvid| 人成视频在线观看免费观看| 丰满人妻熟妇乱又伦精品不卡| 岛国在线观看网站| 亚洲成人国产一区在线观看| tocl精华| 亚洲国产欧美网| videos熟女内射| 免费高清在线观看日韩| 久久精品国产99精品国产亚洲性色 | 天堂俺去俺来也www色官网| 两个人免费观看高清视频| a级毛片在线看网站| 色尼玛亚洲综合影院| 亚洲欧美激情在线| 十八禁网站网址无遮挡| 成人国产av品久久久| 99国产综合亚洲精品| 国产一区二区 视频在线| 国产无遮挡羞羞视频在线观看| 黄片小视频在线播放| 亚洲熟女精品中文字幕| 国产免费视频播放在线视频| 黑人巨大精品欧美一区二区蜜桃| 91成年电影在线观看| 夜夜夜夜夜久久久久| 不卡av一区二区三区| 日韩精品免费视频一区二区三区| 嫁个100分男人电影在线观看| 正在播放国产对白刺激| 欧美日韩黄片免| 日韩三级视频一区二区三区| 亚洲成a人片在线一区二区| 一区二区三区乱码不卡18| 免费看a级黄色片| 狠狠婷婷综合久久久久久88av| 18禁裸乳无遮挡动漫免费视频| 无限看片的www在线观看| 啦啦啦免费观看视频1| 国产精品 国内视频| 久久影院123| 久久久久久免费高清国产稀缺| 美女高潮喷水抽搐中文字幕| 亚洲人成电影免费在线| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 如日韩欧美国产精品一区二区三区| 91老司机精品| 乱人伦中国视频| 51午夜福利影视在线观看| 国产精品欧美亚洲77777| 在线观看66精品国产| 国产一区二区三区综合在线观看| 国产成人免费观看mmmm| 麻豆成人av在线观看| 十八禁网站网址无遮挡| 亚洲综合色网址| 国产视频一区二区在线看| 黄色视频不卡| 国产精品免费大片| 变态另类成人亚洲欧美熟女 | 啦啦啦在线免费观看视频4| 99久久人妻综合| 国产福利在线免费观看视频| 国产在线精品亚洲第一网站| 日韩欧美国产一区二区入口| 在线看a的网站| 男女高潮啪啪啪动态图| 精品亚洲乱码少妇综合久久| 亚洲av第一区精品v没综合| 成人影院久久| 天天躁夜夜躁狠狠躁躁| 精品免费久久久久久久清纯 | 欧美日韩av久久| 午夜福利,免费看| 午夜老司机福利片| 少妇 在线观看| 亚洲精品美女久久久久99蜜臀| 国产欧美日韩综合在线一区二区| 女性被躁到高潮视频| 精品一品国产午夜福利视频| 成人黄色视频免费在线看| 男男h啪啪无遮挡| 女人高潮潮喷娇喘18禁视频| 久久久久久久久久久久大奶| 国产99久久九九免费精品| 91成人精品电影| 亚洲国产欧美一区二区综合| 亚洲av成人一区二区三| 人妻一区二区av| 国产欧美日韩综合在线一区二区| 成年版毛片免费区| 国产免费现黄频在线看| 欧美久久黑人一区二区| 搡老岳熟女国产| 免费人妻精品一区二区三区视频| 99九九在线精品视频| 80岁老熟妇乱子伦牲交| 热99re8久久精品国产| 手机成人av网站| 99在线人妻在线中文字幕 | 国产成人精品久久二区二区免费| 国产真人三级小视频在线观看| 国产伦人伦偷精品视频| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 国产精品偷伦视频观看了| 啦啦啦在线免费观看视频4| 久久国产精品男人的天堂亚洲| 亚洲黑人精品在线| 午夜日韩欧美国产| av天堂在线播放| 欧美精品一区二区大全| 熟女少妇亚洲综合色aaa.| 777久久人妻少妇嫩草av网站| 亚洲成人国产一区在线观看| 99热网站在线观看| 日本黄色视频三级网站网址 | 视频区图区小说| 久久精品亚洲精品国产色婷小说| 亚洲,欧美精品.| 性少妇av在线| 国产精品一区二区精品视频观看| 国产男女超爽视频在线观看| 亚洲国产中文字幕在线视频| 国产男靠女视频免费网站| 丝袜人妻中文字幕| 色老头精品视频在线观看| 一二三四社区在线视频社区8| 两人在一起打扑克的视频| 国产视频一区二区在线看| 天天添夜夜摸| 精品国产乱子伦一区二区三区| 女人高潮潮喷娇喘18禁视频| 99香蕉大伊视频| 老熟妇乱子伦视频在线观看| 黄色视频,在线免费观看| 伦理电影免费视频| 午夜福利视频在线观看免费| av福利片在线| 精品国产亚洲在线| 欧美精品一区二区大全| 高清欧美精品videossex| av视频免费观看在线观看| 欧美日韩一级在线毛片| 激情视频va一区二区三区| 男女之事视频高清在线观看| 午夜激情久久久久久久| a级毛片黄视频| 十八禁网站网址无遮挡| 美女扒开内裤让男人捅视频| 日韩成人在线观看一区二区三区| 免费观看a级毛片全部| 欧美日韩av久久| 露出奶头的视频| 久久午夜综合久久蜜桃| 男女无遮挡免费网站观看| 国产欧美日韩一区二区精品| 侵犯人妻中文字幕一二三四区| 大香蕉久久网| 亚洲欧美日韩高清在线视频 | 9色porny在线观看| 精品一品国产午夜福利视频| 两人在一起打扑克的视频| 在线观看www视频免费| 久久久久久久大尺度免费视频| 男女之事视频高清在线观看| 国产成人精品久久二区二区免费| 一本一本久久a久久精品综合妖精| 日本vs欧美在线观看视频| 成人特级黄色片久久久久久久 | av国产精品久久久久影院| 日本欧美视频一区| 99精品久久久久人妻精品| 在线观看66精品国产| 精品国产乱码久久久久久男人| 国产高清国产精品国产三级| 精品视频人人做人人爽| 女人精品久久久久毛片| 国产精品影院久久| 黄色片一级片一级黄色片| 亚洲精品成人av观看孕妇| 黄色视频,在线免费观看| 淫妇啪啪啪对白视频| 午夜福利乱码中文字幕| 久久久久久亚洲精品国产蜜桃av| 久久久久久人人人人人| 国产黄色免费在线视频| 欧美国产精品va在线观看不卡| 亚洲久久久国产精品| 色婷婷av一区二区三区视频| 动漫黄色视频在线观看| 高清毛片免费观看视频网站 | 色婷婷av一区二区三区视频| 亚洲国产欧美日韩在线播放| 一区二区三区精品91| av电影中文网址| 大香蕉久久网| 人人澡人人妻人| 精品一区二区三卡| 操美女的视频在线观看| 欧美+亚洲+日韩+国产| 大香蕉久久成人网| 麻豆成人av在线观看| 最近最新中文字幕大全免费视频| 咕卡用的链子| 亚洲av电影在线进入| 91老司机精品| 深夜精品福利| 国产日韩一区二区三区精品不卡| 国产主播在线观看一区二区| 成人三级做爰电影| 精品国产一区二区三区四区第35| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| 久久精品成人免费网站| 侵犯人妻中文字幕一二三四区| 国产有黄有色有爽视频| xxxhd国产人妻xxx| 精品人妻熟女毛片av久久网站| 亚洲avbb在线观看| 日韩精品免费视频一区二区三区| 久久精品国产99精品国产亚洲性色 | 99国产精品免费福利视频| 精品午夜福利视频在线观看一区 | 十八禁人妻一区二区| 欧美日韩av久久| 男女之事视频高清在线观看| 99国产精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 久久热在线av| 国产精品免费一区二区三区在线 | 黄片播放在线免费| 美女高潮喷水抽搐中文字幕| 精品一区二区三区av网在线观看 | 国产单亲对白刺激| 亚洲国产毛片av蜜桃av| 精品一区二区三区四区五区乱码| 久久久久久人人人人人| 亚洲人成电影免费在线| 男人舔女人的私密视频| 一进一出抽搐动态| 亚洲 国产 在线| 大香蕉久久成人网| 丰满少妇做爰视频| 在线观看免费视频网站a站| 国产精品偷伦视频观看了| 美女国产高潮福利片在线看| 欧美在线一区亚洲| 一边摸一边抽搐一进一出视频| 精品久久久久久久毛片微露脸| 亚洲欧美日韩高清在线视频 | 精品高清国产在线一区| 亚洲av电影在线进入| 91精品国产国语对白视频| 国产一卡二卡三卡精品| 男女高潮啪啪啪动态图| 1024香蕉在线观看| 欧美久久黑人一区二区| 一二三四社区在线视频社区8| 正在播放国产对白刺激| 国产精品免费一区二区三区在线 | 国精品久久久久久国模美| 久久久精品免费免费高清| 亚洲成a人片在线一区二区| 亚洲伊人色综图| 无人区码免费观看不卡 | 久久久精品国产亚洲av高清涩受| 国产日韩欧美视频二区| 青草久久国产| 999久久久精品免费观看国产| 午夜福利在线观看吧| 亚洲色图av天堂| 亚洲av日韩在线播放| 欧美乱码精品一区二区三区| 免费观看av网站的网址| 一本大道久久a久久精品| 国产精品自产拍在线观看55亚洲 | 久久人人97超碰香蕉20202| 午夜精品久久久久久毛片777| 亚洲中文字幕日韩| tocl精华| 好男人电影高清在线观看| bbb黄色大片| 男人舔女人的私密视频| 免费不卡黄色视频| 丝袜喷水一区| av网站在线播放免费| 精品国产一区二区三区四区第35| 可以免费在线观看a视频的电影网站| 国产在线一区二区三区精| 欧美日韩一级在线毛片| 乱人伦中国视频| 性少妇av在线| 天天影视国产精品| www.精华液| 一区二区日韩欧美中文字幕| 久久久久精品国产欧美久久久| 一区福利在线观看| 欧美大码av| 欧美黄色片欧美黄色片| 精品国内亚洲2022精品成人 | 久久久久精品人妻al黑| 91精品国产国语对白视频| 日韩三级视频一区二区三区| 亚洲av片天天在线观看| 欧美日韩亚洲综合一区二区三区_| 精品卡一卡二卡四卡免费| 99国产综合亚洲精品| 国产成+人综合+亚洲专区| 亚洲av第一区精品v没综合| 五月天丁香电影| 亚洲欧美一区二区三区黑人| 黄色视频在线播放观看不卡| 久久人妻福利社区极品人妻图片| 久久99热这里只频精品6学生| 深夜精品福利| 久久 成人 亚洲| 亚洲精品美女久久av网站| 国产免费视频播放在线视频| 人人妻人人澡人人爽人人夜夜| 在线观看免费高清a一片| 亚洲精品粉嫩美女一区| 精品卡一卡二卡四卡免费| 国内毛片毛片毛片毛片毛片| 欧美日韩av久久| 亚洲一区二区三区欧美精品| 中文亚洲av片在线观看爽 | 夫妻午夜视频| 他把我摸到了高潮在线观看 | 色94色欧美一区二区| 黑人操中国人逼视频| 狠狠婷婷综合久久久久久88av| 不卡av一区二区三区| 日本a在线网址| 成人高潮视频无遮挡免费网站| 成人午夜高清在线视频| h日本视频在线播放| 在线观看午夜福利视频| 999精品在线视频| 狂野欧美白嫩少妇大欣赏| 久久久久亚洲av毛片大全| 首页视频小说图片口味搜索| 国产精品av视频在线免费观看| 精品久久久久久成人av| 亚洲av美国av| 成人欧美大片| 观看美女的网站| 亚洲色图av天堂| 一本久久中文字幕| 国产爱豆传媒在线观看| 国产69精品久久久久777片 | 亚洲人成网站高清观看| 久久久国产精品麻豆| 亚洲国产中文字幕在线视频| 日本 av在线| 久久国产精品人妻蜜桃| 一二三四社区在线视频社区8| 亚洲在线观看片| 99在线人妻在线中文字幕| 国产亚洲精品综合一区在线观看| 国产精品av久久久久免费| 国产亚洲精品一区二区www| 成人午夜高清在线视频| 亚洲电影在线观看av| 国产美女午夜福利| 91av网一区二区| 国产av一区在线观看免费| 一二三四社区在线视频社区8| 欧美最黄视频在线播放免费| 国产精品av久久久久免费| 一个人看视频在线观看www免费 | 蜜桃久久精品国产亚洲av| 免费电影在线观看免费观看| 99国产精品99久久久久| 少妇裸体淫交视频免费看高清| 欧美精品啪啪一区二区三区| www.www免费av| 日韩欧美 国产精品| 日本黄色片子视频| 久久99热这里只有精品18| 夜夜躁狠狠躁天天躁| 露出奶头的视频| 亚洲国产精品合色在线| 亚洲成av人片免费观看| 波多野结衣高清无吗| 精品日产1卡2卡| 黄片小视频在线播放| 啦啦啦韩国在线观看视频| 叶爱在线成人免费视频播放| 亚洲午夜精品一区,二区,三区| 中亚洲国语对白在线视频| 老汉色∧v一级毛片| 三级毛片av免费| 国产真人三级小视频在线观看| 欧美国产日韩亚洲一区| 麻豆国产97在线/欧美| 免费在线观看视频国产中文字幕亚洲| a级毛片在线看网站| 999久久久精品免费观看国产| 亚洲成人久久爱视频| 搡老岳熟女国产| 成人18禁在线播放| 国产精品日韩av在线免费观看| 国产午夜精品论理片| 麻豆av在线久日| 麻豆久久精品国产亚洲av| 国产69精品久久久久777片 | 免费观看精品视频网站| 在线观看美女被高潮喷水网站 | 蜜桃久久精品国产亚洲av| 波多野结衣高清作品| 91在线观看av| 这个男人来自地球电影免费观看| 午夜视频精品福利| 性欧美人与动物交配| 国产私拍福利视频在线观看| 人人妻人人看人人澡| 免费搜索国产男女视频| 一进一出抽搐gif免费好疼| 老司机午夜福利在线观看视频| 久久婷婷人人爽人人干人人爱| 黑人欧美特级aaaaaa片| 久久中文字幕一级| h日本视频在线播放| 亚洲人成电影免费在线| 精品久久久久久,| 一本精品99久久精品77| 听说在线观看完整版免费高清| 91九色精品人成在线观看| 日韩av在线大香蕉| 1000部很黄的大片| 亚洲精品456在线播放app | 精品久久久久久久久久免费视频| 午夜福利免费观看在线| 久久久久国产精品人妻aⅴ院| 一本久久中文字幕| 国产野战对白在线观看| www.熟女人妻精品国产| 国产一区二区在线av高清观看| 香蕉久久夜色| 免费观看的影片在线观看| 亚洲美女黄片视频| 国产aⅴ精品一区二区三区波| 欧美一区二区国产精品久久精品| 久久久久免费精品人妻一区二区| 国产精品 国内视频| 悠悠久久av| 两人在一起打扑克的视频| 久久精品aⅴ一区二区三区四区| 夜夜躁狠狠躁天天躁| 搡老妇女老女人老熟妇| 欧美日韩黄片免| 日本黄色片子视频| 国产成年人精品一区二区| 老司机深夜福利视频在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美+亚洲+日韩+国产| 午夜福利成人在线免费观看| 九色成人免费人妻av| 亚洲最大成人中文| 亚洲欧美日韩高清专用| 97人妻精品一区二区三区麻豆| 国产高清videossex| 国产一区二区在线av高清观看| 天天一区二区日本电影三级| 后天国语完整版免费观看| 人人妻人人澡欧美一区二区| xxx96com| 狠狠狠狠99中文字幕| 两个人的视频大全免费| 精品久久久久久成人av| 久久久久久大精品| 久久香蕉精品热| 性欧美人与动物交配| 十八禁网站免费在线| 少妇熟女aⅴ在线视频| 99在线视频只有这里精品首页| 亚洲国产精品999在线| 免费在线观看日本一区| 一夜夜www| 日韩欧美 国产精品| 日韩欧美在线二视频| 露出奶头的视频| 亚洲成人免费电影在线观看| www国产在线视频色| 长腿黑丝高跟| 国产黄片美女视频| 免费在线观看成人毛片| 亚洲av免费在线观看| 亚洲最大成人中文| 亚洲人与动物交配视频| 国产精品久久久人人做人人爽| 波多野结衣巨乳人妻| 99riav亚洲国产免费| 少妇人妻一区二区三区视频| 成人一区二区视频在线观看| 最近最新免费中文字幕在线| 国产精品影院久久| 亚洲av成人精品一区久久| 99国产精品一区二区蜜桃av| 999久久久国产精品视频| 亚洲男人的天堂狠狠| 舔av片在线| 亚洲国产高清在线一区二区三| 国产97色在线日韩免费| 亚洲av免费在线观看| 婷婷亚洲欧美| 在线看三级毛片| 亚洲熟妇中文字幕五十中出| 国语自产精品视频在线第100页| 狂野欧美白嫩少妇大欣赏| 免费电影在线观看免费观看| 在线十欧美十亚洲十日本专区| 久久久国产成人免费| 国产精品av视频在线免费观看| 中文字幕高清在线视频| 免费看十八禁软件| 国产男靠女视频免费网站| 中文字幕久久专区| 久久天躁狠狠躁夜夜2o2o| 无遮挡黄片免费观看| 精品一区二区三区av网在线观看| 51午夜福利影视在线观看| 12—13女人毛片做爰片一| 黄色丝袜av网址大全| 18禁美女被吸乳视频| 国产三级中文精品| 国产在线精品亚洲第一网站| 免费搜索国产男女视频| 国产三级在线视频|