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

    三種數(shù)值模式對四川省汛期降水預(yù)報性能的檢驗(yàn)

    2020-03-30 02:17:30屠妮妮何光碧衡志煒吳蓬萍
    高原山地氣象研究 2020年4期
    關(guān)鍵詞:中雨實(shí)況小雨

    屠妮妮,何光碧,衡志煒,吳蓬萍

    (中國氣象局成都高原氣象研究所/高原與盆地暴雨旱澇災(zāi)害四川省重點(diǎn)實(shí)驗(yàn)室, 成都 610072)

    引言

    ECMWF、GRAPES和WRF模式預(yù)報產(chǎn)品被廣泛的應(yīng)用在我國天氣預(yù)報業(yè)務(wù)中[1-7]。張宏芳等[8]利用不同檢驗(yàn)方法對比分析了ECMWF模式和日本高分辨率模式對降水的預(yù)報能力,通過對檢驗(yàn)結(jié)果進(jìn)行分析,發(fā)現(xiàn)調(diào)整閾值改變預(yù)報偏差能夠在一定程度上提高降水的預(yù)報技巧。熊秋芬[9]利用經(jīng)典的TS、預(yù)報偏差、去除隨機(jī)事件后的公平T評分(ETS)及真實(shí)技巧評分TSS等統(tǒng)計檢驗(yàn)方法,分析了GRAPES模式在中國不同地區(qū)降水的預(yù)報能力,GRAPES模式預(yù)報5mm以上降水在長江以南多于觀測,尤其是在西南地區(qū)東部預(yù)報次數(shù)明顯偏多。SWCWARMS模式在成都區(qū)域中心業(yè)務(wù)運(yùn)行多年,其多時段的降水預(yù)報產(chǎn)品是業(yè)務(wù)預(yù)報的重要參考依據(jù),與ECMWF模式相比,SWCWARMS高能暖區(qū)型暴雨和大暴雨ETS評分高于ECMWF模式,斜壓鋒生型大暴雨優(yōu)于ECMWF模式[10]。張琪等[11]分析SWCWARMS模式對四川盆地兩次強(qiáng)降水開展天氣學(xué)檢驗(yàn),指出SWCWARMS模式在強(qiáng)降水預(yù)報中有3大優(yōu)勢。屠妮妮等[12]每年對比分析SWCWARMS模式和GRAPES_MESO模式在西南地區(qū)的預(yù)報能力,SWCWARMS模式在多個量級降水預(yù)報上優(yōu)于GRAPES模式。

    本文選取SWCWARMS模式、歐洲中心高分辨率模式和GRAPES_MESO模式,對2017~2019年四川省5~10月降水預(yù)報進(jìn)行檢驗(yàn),有助于科研人員和業(yè)務(wù)工作者了解三個模式預(yù)報性能的差異。

    1 資料與方法

    數(shù)值模式資料來源2017、2018和2019年5~10月每日00時起報的降水場預(yù)報結(jié)果, SWCWARMS模式分辨率為9km,ECMWF模式分辨率為0.125°×0.125°,GRAPES_MESO(簡稱為GRAPES)模式分辨率為0.1°×0.1°,三個模式為每日00時起報的0~24h、12~36h、24~48h、36~60h、48~72h的24h累積雨量。實(shí)況降水選取2017年、2018年和2019年連續(xù)三年5~10月每日00~00時、12~12時(文中時間均為世界時)四川省內(nèi)國家站24h累積降水。

    統(tǒng)計檢驗(yàn)基于MET檢驗(yàn)系統(tǒng)的站點(diǎn)檢驗(yàn)方法,插值方法為Unweighted mean方法,該插值方法與DW_MEAN近似,但忽略距離權(quán)重。24h累積降水評分等級分別為≥0.1mm(小雨)、≥10mm(中雨)、≥25mm(大雨)、≥50mm(暴雨)、≥100mm(大暴雨),統(tǒng)計量包括TS評分、BS偏差、ETS、FAR和TSS,TS評分對小概率事件的評分趨向于氣候概率,忽略了ND的影響;BS評分用來測量事件的頻率而不考慮預(yù)報的準(zhǔn)確率,當(dāng)偏差>1說明該量級降水預(yù)報面積大于實(shí)況即出現(xiàn)空報,反之則為漏報,偏差=1時與實(shí)況一致;ETS評分表示除去隨機(jī)偶然事件后,事件發(fā)生且正確的樣本數(shù)與事件發(fā)生總樣本數(shù)的比率;TSS評分能反映成功預(yù)報與空報的差異,可用來全面評價某級降水的總效果,即它不但包括降水事件的預(yù)報效果,還包括對非降水事件的預(yù)報效果,但對小概率事件趨向于POD(命中率),當(dāng)TSS=0時表示無預(yù)報能力,計算公式分別如下:

    TS評分:

    (1)

    BS偏差:

    (2)

    空報率:

    (3)

    ETS評分:

    (4)

    其中

    (5)

    (6)

    (7)

    公式(1)~(7)中,NA、NB、NC和ND分別為預(yù)報正確,空報,漏報和預(yù)報及觀測都未出現(xiàn)降水,NA可表示為預(yù)報時段內(nèi)預(yù)報和實(shí)況降雨量都出現(xiàn)的站點(diǎn)數(shù),即預(yù)報命中的站數(shù);NB表示預(yù)報出現(xiàn)降水而實(shí)況未出現(xiàn)的站點(diǎn)數(shù),即空報站數(shù);NC表示未預(yù)報而實(shí)況出現(xiàn)降水的站點(diǎn)數(shù),即漏報站數(shù);ND表示觀測和預(yù)報都未出現(xiàn)降水的站點(diǎn)數(shù)。

    2 結(jié)果分析

    2.1 總檢驗(yàn)

    按照<0.1mm、0.1~9.99mm、10~24.99mm、25~49.99mm、50~99.99mm、≥100mm上述6個標(biāo)準(zhǔn)將三個模式預(yù)報測站雨量及各模式對應(yīng)實(shí)況雨量分別進(jìn)行統(tǒng)計,連續(xù)552天統(tǒng)計結(jié)果顯示,SWCWARMS模式5個預(yù)報時效上累計樣本數(shù)分別為79439個(0~24h預(yù)報時效上)、79018個(12~36h)、79440個(24~48h)、79018個(36~60h)、79437個(48~72h),GRAPES模式分別為79968個、79092個、79517個、78870個、79271個,ECMWF模式分別為79595個、78521個、79486個、78428個、79288個,樣本總數(shù)均在7萬以上,SWCWARMS模式24h和48h樣本數(shù)略少于GRAPES及ECMWF模式,ECMWF模式36h和60h樣本數(shù)略偏少;格點(diǎn)數(shù)據(jù)插值至站點(diǎn)后出現(xiàn)3個模式在同一預(yù)報時效上樣本數(shù)不一致情況,可能與3個模式分辨率有關(guān)。由三個模式5個預(yù)報時效上分級降水樣本數(shù)占總樣本數(shù)比重可見(表1),小雨比重最大,其次為無降水或微量降水,SWCWARMS、GRAPES和ECMWF模式預(yù)報小雨和中雨都多于實(shí)況,SWCWARMS、GRAPES和ECMWF模式小雨比重較實(shí)況分別偏多3%~8%、1%~4%和18%~20%,中雨較實(shí)況偏多3%~5%、2%~7%和4%~6%;SWCWARMS和GRAPES模式預(yù)報大雨和暴雨多于實(shí)況、ECMWF模式與實(shí)況相當(dāng)、偏差在-0.1%~0.21%,SWCWARMS模式預(yù)報大暴雨多于實(shí)況、GRAPES模式和ECMWF模式較實(shí)況偏少。綜上可見,SWCWARMS模式小到大暴雨降水范圍大于實(shí)況,GRAPES模式小到暴雨降水范圍大于實(shí)況、大暴雨多漏報,ECMWF模式小雨和大雨降水范圍大于實(shí)況、大到大暴雨多漏報,三個模式無降水或微量降水均少于實(shí)況,即模式預(yù)報降水較實(shí)況偏多。

    表1 5個預(yù)報時效上分級降水樣本占比最小值和最大值分布

    Q-Q圖用于檢驗(yàn)數(shù)據(jù)分布的相似性,橫坐標(biāo)與縱坐標(biāo)組成的散點(diǎn)圖代表同一個累計概率所對應(yīng)的分位數(shù),如果兩組數(shù)據(jù)分布相似,則該Q-Q圖散度趨近于落在y=x線對角線上。比較三個模式0~250mm降水的Q-Q圖分布可見(圖1),在0~24h預(yù)報時效上降水分布顯示,SWCWARMS模式0~250mm降水大于實(shí)況,0~200mm降水偏差逐漸增大,該模式預(yù)報雨強(qiáng)大于實(shí)況;GRAPES模式50mm以下降水強(qiáng)度大于實(shí)況,暴雨及以上量級降水強(qiáng)度以小于實(shí)況為主;ECMWF模式小雨和中雨與實(shí)況接近,大雨及以上量級降水以小于實(shí)況為主,150~200mm降水相比前兩模式與實(shí)況最接近。12~36h預(yù)報時效上降水分布顯示,SWCWARMS模式0~250mm降水強(qiáng)度大于實(shí)況,暴雨和大暴雨偏大較多;GRAPES模式小到大雨大于實(shí)況、暴雨及以上量級降水強(qiáng)度小于實(shí)況,尤其是150mm以上量級降水;ECMWF模式小雨和中雨降水強(qiáng)度與實(shí)況接近,大雨及以上量級降水較實(shí)況偏小,200mm以上降水大于實(shí)況。在24~72h預(yù)報時效上(圖略),SWCWARMS模式預(yù)報雨強(qiáng)偏強(qiáng),GRAPES模式暴雨及以上量級降水偏弱,ECMWF模式大雨及以上量級降水偏弱,且強(qiáng)度弱于GRAPES模式。

    2.2 分級降水檢驗(yàn)

    由四川省2017~2019年5~10月降水TS評分可見(圖2a),3個模式隨降水量級增加TS評分逐漸降低,SWCWARMS模式小雨評分分布在0.605~0.656,ECMWF模式小雨評分分布在0.628~0.658,GRAPES模式小雨評分在0.57~0.663。小雨到暴雨TS評分上,ECMWF模式評分最高(除0~24h小雨和12~36h暴雨外),SWCWARMS模式次之,GRAPES模式評分最低,0~24h預(yù)報時效上GRAPES模式小雨評分最高,12~36h預(yù)報時效上SWCWARMS模式暴雨評分最高;對于大暴雨預(yù)報,24h、48h和60h預(yù)報時效上SWCWARMS模式評分最高,其它預(yù)報時效上ECMWF模式評分最高,GRAPES模式評分以低于其它兩個模式為主,其中36~60h大暴雨預(yù)報上,GRAPES模式評分低于SWCWARMS模式但略高于ECMWF模式。

    由3個模式降水偏差分布可見(圖2b), SWCWARMS模式預(yù)報五個量級降水范圍均大于實(shí)況,GRAPES模式預(yù)報小到暴雨預(yù)報降水范圍大于實(shí)況,大暴雨以漏報為主,ECMWF模式預(yù)報小雨、中雨和48~72h大雨降水范圍大于實(shí)況,大雨(0~60h)、暴雨和大暴雨預(yù)報降水范圍小于實(shí)況。在小雨預(yù)報上,ECMWF模式預(yù)報降水落區(qū)范圍偏大最多,SWCWARMS模式次之,GEAPES相對偏差略小。對于中雨而言,12~36h時效上GRAPES模式偏差最大,其它預(yù)報時段SWCWARMS模式偏差最大,ECMWF模式中雨空報少于區(qū)域模式。對大雨和暴雨而言,SWCWARMS模式預(yù)報降水落區(qū)偏大最多,GRAPES模式次之,ECMWF模式預(yù)報落區(qū)偏小;對于大暴雨預(yù)報,SWCWARMS模式空報最多,ECMWF模式在36h和48h漏報最多,GRAPES模式在24h、60h和72h漏報最多。SWCWARMS模式各級降水多空報,5個預(yù)報時效上分別有中雨(60h)、大雨(36h、48h、72h)和暴雨(24h)空報最多,GRAPES模式對小到暴雨多空報而大暴雨多漏報,中雨(60h)和大雨(24h、36h、48h、72h)空報最多,ECMWF模式對小雨和中雨表現(xiàn)為空報多而其它量級降水多漏報,小雨空報最多。

    由3個模式ETS評分分布可見(圖2c),隨降水強(qiáng)度增強(qiáng)評分逐漸降低,SWCWARMS模式小雨ETS評分分布在0.135~0.215,小雨評分高于其它降水量級,ECMWF模式小雨評分分布在0.105~0.142,中雨評分高于其它降水量級,GRAPES模式小雨評分在0.139~0.207,小雨評分高于其它降水量級。在小雨預(yù)報上,SWCWARMS模式24h評分最高,GRAPES模式次之,ECMWF模式評分最低,在36~72h預(yù)報時效上,GRAPES模式評分最高,SWCWARMS模式次之,ECMWF模式評分最低;GRAPES模式小雨優(yōu)于SWCWARMS和ECMWF模式,SWCWARMS模式優(yōu)于ECMWF模式。在中雨、大雨和暴雨預(yù)報上,ECMWF模式評分最高,SWCWARMS模式次之,GRAPES模式評分最低,但在60h暴雨預(yù)報上,SWCWARMS模式EST評分最低;ECMWF模式中到暴雨預(yù)報優(yōu)于其它兩個模式,SWCWARMS模式略優(yōu)于GRAPES模式。在大暴雨預(yù)報上,SWCWARMS模式在24h、48h、60h評分最高,ECMWF在36h和72h評分最高,GRAPES模式在24h、36h、72h評分最低,ECMWF在48h和60h評分最低;SWCWARMS模式大暴雨預(yù)報在部分時次上略優(yōu)于其它兩個模式。GRAPES模式小雨評分較高,ECMWF模式中到暴雨評分較高,SWCWARMS模式大暴雨部分時次上優(yōu)于其它兩模式,SWCWARMS模式中到大暴雨ETS評分以高于GRAPES模式為主。

    由3個模式TSS評分分布可見(圖2d),對小雨而言,SWCWARMS模式在24h和36h預(yù)報時效上評分最高,小雨評分分布在0.23~0.344,小雨評分高于其它量級降水;GRAPES模式在48h、60h和72h預(yù)報時效上小雨評分最高,介于0.234~0.311,小雨評分高于其它量級降水;ECMWF模式小雨評分最低,分布在0.181~0.234。對中雨而言,ECMWF模式評分最高,該模式中雨評分高于其它量級降水,SWCWARMS模式次之,GRAPES模式評分最低。對大雨而言,SWCWARMS模式在24h、36h和48h評分最高,ECMWF模式在60h和72h評分最高,GRAPES模式評分最低。對暴雨而言,在24h、36h、48h和60h預(yù)報時效上SWCWARMS模式評分最高,ECMWF模式在72h評分最高,除24h時效上ECMWF模式評分最低,其它時效上GRAPES模式評分最低。對大暴雨而言,SWCWARMS模式在24h、36h、48h、60h評分最高,GRAPES模式72h預(yù)報時效上評分最高,ECMWF模式在48h、60h、72h評分最低,GRAPES模式24h和36h預(yù)報時效上評分最低。SWCWARMS模式在多個預(yù)報時效上大到大暴雨TSS評分高于ECMWF模式,中到大暴雨評分高于GRAPES模式。

    綜上可見,ECMWF模式對小到大雨具有較好的預(yù)報能力,SWCWARMS模式對暴雨和大暴雨在部分量級降水預(yù)報上優(yōu)于ECMWF模式,GRAPES模式相對上述兩模式偏差。SWCWARMS模式預(yù)報降水范圍大于實(shí)況,該模式中到暴雨空報多于ECMWF和GRAPES模式,ECMWF模式小雨和中雨表現(xiàn)為空報多而其它量級降水多漏報,GRAPES模式對小到暴雨多空報而大暴雨表現(xiàn)為漏報,GRAPES模式中到暴雨空報多于ECMWF模式,大暴雨在24h、60h和72h預(yù)報時效上漏報多于ECMWF模式。

    2.3 年際變化

    比較2017~2019年5~10月三個模式對四川省降水預(yù)報≥0.1mm檢驗(yàn)結(jié)果的年變化(圖3),由TS評分可見(圖3a),SWCWARMS模式24h預(yù)報時效上2018年評分最低,2019年評分最高,36~60h預(yù)報時效上2017年評分最低,2019年評分最高,72h預(yù)報時效上2017年評分最低,2018年評分最高;ECMWF模式5個預(yù)報時效上2017年評分最低,2019年評分最高;GRAPES模式36~72h預(yù)報時效上2017年評分最低,2019年評分最高,24h預(yù)報時效上2018年評分最低,2019年評分最高。2017年,SWCWARMS模式評分在24~60h評分高于其它兩模式,ECMWF模式72h評分高于其它兩模式,GRAPES模式在36~72h評分偏低,ECMWF模式在24h評分偏低,在2017年SWCWARMS模式評分以高于其它兩模式為主;在2018年24h、36h、48h預(yù)報時效上,SWCWARMS模式評分最高、GRAPES模式評分最低,60h和72h預(yù)報時效上ECMWF模式評分高于SWCWARMS和GRAPES模式,GRAPES模式評分最低,2018年SWCWARMS模式優(yōu)勢在前48h降水預(yù)報上,GRAPES模式預(yù)報始終差于ECMWF和SWCWARMS模式;在2019年24h和36h時效上GRAPES模式評分最高、SWCWARMS模式評分最低,48~72h時效上ECMWF模式評分最高、SWCWARMS模式評分最低,2019年ECMWF模式較2017年TS評分提高7%左右,GRAPES模式TS評分提高8%~11%,而SWCWARMS模式TS評分僅提高1%~2%。

    由BS預(yù)報偏差分布可見(圖3b),3個模式預(yù)報降水范圍均大于實(shí)況,誤差分布在1.1~2.2。SWCWARMS模式在24h預(yù)報時效上,2017年空報多于2018年和2019年而2019年空報較前兩年略有減少,36~72h預(yù)報時效上,2018年空報最多而2019年空報最少;ECMWF及GRAPES模式在24h、48h、60h和72h預(yù)報時效上,2019年空報多于前兩年而2017年空報少于后兩年,ECMWF及GRAPES模式在36h預(yù)報時效上,2018年空報最多而2017年空報最少。2017~2019年這3年期間SWCWARMS模式降水空報先增大后減小,ECMWF和GRAPES模式降水空報逐年增大的趨勢。ECMWF模式BS值始終大于SWCWARMS和GRAPES模式,2017年和2018年GRAPES模式偏差最小,2018年24h預(yù)報時效上及2019年SWCWARMS模式BS值小于其它兩模式;這3年間SWCWARMS模式誤差先增大后減小而ECMWF和GRAPES模式空報逐年增大。

    由ETS評分分布可見(圖3c),3個模式ETS評分隨預(yù)報時效延長評分降低,僅ECMWF模式在2018年72h時效上ETS評分略高于60h,3個模式ETS評分均表現(xiàn)為2018年評分低于2017年和2019年,其中2018年72h時效上ETS評分大于2017年,2019年評分高于2017年。對于2017年而言,SWCWARMS模式在60h和72h評分最高,GRAPES模式在24h、36h、48h評分最高,ECMWF模式評分最低;對于2018年和2019年,SWCWARMS模式在24h評分最高,GRAPES模式在36h、48h、60h、72h評分最高,ECMWF模式評分最低。3個模式2019年ETS評分均高于前兩年,GRAPES模式評分以高于SWCWARMS模式為主,ECMWF模式ETS評分始終低于GRAPES和SWCWARMS模式。

    由TSS評分分布可見(圖3d),3個模式TSS評分隨預(yù)報時效延長評分整體呈降低趨勢,SWCWARMS模式24~72h、ECMWF模式24~60h、GRAPES模式24~60h預(yù)報時效上2019年評分高于前兩年、2018年評分低于2017年和2019年,ECMWF和GRAPES模式72h時效上2019年評分最高而2017年評分最低。ECMWF模式TSS評分低于其它兩個模式,SWCWARMS模式在2017年60h預(yù)報時效上、2018年24h預(yù)報時效上、2019年5個預(yù)報時效上TSS評分高于GRAPES模式。3個模式TSS評分整體呈2019年評分最高、2018年評分最低分布,GRAPES模式在2017年和2018年以優(yōu)于SWCWARMS和ECMWF模式為主,SWCWARMS模式2019年評分最高,ECMWF模式評分始終低于其它兩模式。

    由降水準(zhǔn)確率分布可見(圖3e),3個模式降水預(yù)報準(zhǔn)確率隨預(yù)報時效延長評分以降低為主,SWCWARMS模式降水準(zhǔn)確率在70%以上,ECMWF模式降水準(zhǔn)確率在68%以上,GRAPES模式降水準(zhǔn)確率在67%以上。SWCWARMS模式在24h和72h預(yù)報時效上、ECMWF模式5個預(yù)報時效上、GRAPES模式24h、48h、60h和72h預(yù)報時效上,2019年較前兩年準(zhǔn)確率最高而2017年最低;SWCWARMS模式在36h、48h和60h時效上及GRAPES模式36h預(yù)報時效上2019年準(zhǔn)確率最高而2018年最低。SWCWARMS模式在2017年48h、60h和72h及2018年24h、48h和72h準(zhǔn)確率高于其它兩模式,2019年48h和72h準(zhǔn)確率低于其它兩模式;ECMWF模式在2019年72h準(zhǔn)確率最高,2017年5個預(yù)報時效、2018年24~60h、2019年24h、36h和60h準(zhǔn)確率低于其它兩模式;GRAPES模式在2017年24h和36h、2018年36h和60h、2018年24~60h準(zhǔn)確率最高,2018年72h準(zhǔn)確率低于其它兩模式。SWCWARMS模式在2017年和2018年與GRAPES模式相當(dāng),2019年降水準(zhǔn)確率低于GRAPES模式,ECMWF模式準(zhǔn)確率以低于上述兩模式為主。

    由降水空報率分布可見(圖3f),3個模式空報率隨預(yù)報時效延長整體呈先增大再減小再增大再減小的趨勢,SWCWARMS模式24h、36h和72h、ECMWF模式24~72h及GRAPES模式60h和72h預(yù)報時效上2019年空報率最低而2017年最高,SWCWARMS模式48h和60h及GRAPES模式24h、36h和48h預(yù)報時效上2019年空報率最低而2018年空報最多。ECMWF模式空報多于其它兩模式,2017年24~72h、2018年36~72h預(yù)報時效上GRAPES模式空報為三個模式最少,2018年24h、2019年5個預(yù)報時效上SWCWARMS模式空報最少。3個模式2019年降水空報都少于前兩年,SWCWARMS模式在2019年相對2017年空報減少最多。

    對于小雨及以上量級降水統(tǒng)計檢驗(yàn)顯示,三個模式2019年預(yù)報效果在多個預(yù)報時效上好于2017年及2018年,準(zhǔn)確率提高而空報減小。ECMWF模式預(yù)報偏差、空報率大于SWCWARMS和GRAPES模式,ETS評分、TSS評分和預(yù)報準(zhǔn)確率低于SWCWARMS和GRAPES模式;在2017年及2018年SWCWARMS模式在多個預(yù)報時效上TS評分、ETS評分、TSS評分和預(yù)報準(zhǔn)確率高于GRAPES模式,但降水空報率也大于GRAPES模式,2019年GRAPES模式TS評分、ETS評分、預(yù)報準(zhǔn)確率、預(yù)報偏差及空報率大于SWCWARMS模式。

    比較三個模式對暴雨及以上量級降水TS評分可見(圖4),2017年和2018年GRAPES模式對暴雨預(yù)報能力差于SWCWARMS和ECMWF模式,2019年ECMWF模式對暴雨預(yù)報能力優(yōu)于SWCWARMS和GRAPES模式,SWCWARMS模式在2017年24h和72h、2018年24h、36h、48h和60h評分高于其它兩模式,SWCWARMS模式在2019年36h、60h和72h時效上評分低于其它兩模式;ECMWF模式在2017年36h、48h和60h、2018年72h及2019年5個預(yù)報時效上評分最高;GRAPES模式在2017年、2018年和2019年24h和48h預(yù)報時效上評分較其它模式偏低。對應(yīng)BS偏差值顯示,SWCWARMS模式BS值均大于1,即該模式預(yù)報暴雨范圍均大于實(shí)況,除2019年60h及72h時效外ECMWF模式預(yù)報暴雨范圍小于實(shí)況即該模式對暴雨及以上量級降水以漏報為主,除2018年72h時效外GRAPES模式預(yù)報暴雨范圍大于實(shí)況,SWCWARMS模式空報大于GRAPES和ECMWF模式。由此可見,GRAPES及ECMWF模式2019年暴雨及以上量級降水預(yù)報改善大于SWCWARMS模式,SWCWARMS模式在這3年期間TS評分波動較小,但2019年降水空報較2017年和2018年顯著降低。

    2.4 逐月降水檢驗(yàn)

    為了進(jìn)一步了解3個模式在汛期不同月份的預(yù)報效果,以月為單位對3個模式的降水預(yù)報進(jìn)行檢驗(yàn)對比分析(圖5),3個模式TS評分峰值出現(xiàn)月份不同,SWCWARMS模式小雨、中雨、大雨和暴雨TS評分最大值分別出現(xiàn)在201809(年月,下同)、201807、201807、201710,TS評分最小值分別出現(xiàn)在201808、201810、201810、201905;ECMWF模式評分最大值分別出現(xiàn)在201907、201907、201909、201807,評分最小值分別出現(xiàn)在201710、201706、201908、201906;GRAPES模式評分最大值分別出現(xiàn)在201806、201809、201809、201908,評分最小值分別出現(xiàn)在201705、201710、201709、201710和201805;SWCWARMS模式2018年10月及GRAPES模式2017年9月暴雨評分缺失,SWCWARMS模式2019年5月暴雨評分為0,GRAPES模式2017年9月大雨,201710、201805、201910暴雨評分為0。對小雨而言,SWCWARMS模式在2017年5~6月、9~10月、2018年9月、2019年5~6月、9月評分都高于GRAPES和ECMWF模式,ECMWF模式在2017年5~6月和9月、2018年5~8月和10月、2019年8月和10月評分大于SWCWARMS模式,GRAPES模式在2017年7~8月、2018年5~8月和10月、2019年8月和10月評分大于SWCWARMS和ECMWF模式。對中雨而言,SWCWARMS模式優(yōu)勢在2017年6月、8月、2018年6~7月、2019年8月,ECMWF模式優(yōu)勢在2017年5月、7月、9月和10月、2018年5月、2019年7月、9~10月,GRAPES模式優(yōu)勢在2018年8~10月、2019年5月和6月。對大雨而言,SWCWARMS模式在2017年6月、8~9月、2018年8月評分高于其它兩模式,ECMWF模式在2017年7月、10月、2018年5~7月、2019年9~10月評分高于其它兩模式,GRAPES模式在2017年5月、2018年9~10月、2019年5~8月評分高于其它兩模式。對暴雨而言,SWCWARMS模式優(yōu)勢在2017年7~10月、2018年6月、8月、2019年9~10月,ECMWF模式優(yōu)勢在2017年5~6月、2018年5月、7月、9月,GRAPES模式優(yōu)勢在2018年10月、2019年5~8月。SWCWARMS模式在2018年8月和2019年8月小雨評分較低,2018年10月中到暴雨評分較低,其它月份降水評分較為穩(wěn)定;ECMWF模式2017年和2018年評分穩(wěn)定,2019年各級降水評分振幅較大;GRAPES模式中到暴雨在2017年秋季和2018年第二季度評分偏低,2018年下半年開始中雨及以上量級降水評分顯著提高。在小雨和中雨預(yù)報上,三個模式各有所長;在大雨預(yù)報上,GRAPES模式評分逐漸增大并高于其它兩模式,SWCWARMS模式評分介于GRAPES和ECMWF模式之間;對于暴雨,ECMWF模式介于GRAPES和SWCWARMS模式之間,2017年SWCWARMS模式優(yōu)于GRAPES模式,2019年反之。

    2.5 四川省東部和西部

    沿103°E將四川省分為東部和西部,比較模式對川東和川西降水的預(yù)報能力。由表2可見,對于SWCWARMS模式及GRAPES模式川西TS、ETS評分高于川東,預(yù)報偏差小于川東;ECMWF模式川西TS評分高于川東,川西ETS和BS低于川東。三個模式均存在預(yù)報過度現(xiàn)象,川東空報大于川西。初步分析川西評分高于川東的原因如下:(1)東部測站多,西部測站數(shù)量少;(2)東部降水預(yù)報偏差大的測站數(shù)多(NB樣本數(shù)多),使得在計算TS和ETS時分母偏大,使得TS和ETS數(shù)值降低。3個模式0~24h降水TS和ETS評分高于其它時段,SWCWARMS模式在川西36h、48h和60h時效TS評分高于其它兩模式,ECMWF模式在川西72h時效、川東24~72h時效上TS評分高于其它兩模式,GRAPES模式在川西24h時效上TS評分高于其它兩模式,SWCWARMS模式在川東24h、36h時效評分低于ECMWF和GRAPES模式,ECMWF模式在川西24h時效上評分偏低,GRAPES模式在川西36h、48h、60h和72h時效、川東48h、60h和72h時效上評分偏低。SWCWARMS模式在川西24h時效、川東5個預(yù)報時效上ETS評分高于其它兩模式,ECMWF模式ETS評分始終低于其它兩模式。3個模式在預(yù)報偏差分布上川東與川西一致,ECMWF模式預(yù)報偏差最大,24h時效上SWCWARMS模式偏差最小,36h、48h、60h和72h時效GRAPES模式偏差最小。

    表2 2017~2019年四川東部和西部24~72h時效三個模式(≥0.1mm)TS、ETS和BS檢驗(yàn)結(jié)果

    3 結(jié)論與討論

    通過SWCWARMS、ECMWF和GRAPES三個模式對2017~2019年四川汛期(5~10月)24h累積降水進(jìn)行檢驗(yàn)分析,得出以下幾點(diǎn)結(jié)論:

    (1)3年平均TS評分顯示ECMWF模式對小到大雨預(yù)報能力優(yōu)于SWCWARMS和GRAPES模式,SWCWARMS模式部分時次上暴雨和大暴雨預(yù)報優(yōu)于ECMWF模式,GRAPES模式TS評分以低于上述兩模式為主。SWCWARMS模式預(yù)報降水范圍大于實(shí)況,ECMWF模式小雨和中雨降水范圍大于實(shí)況而其它量級降水偏少,GRAPES模式對小到暴雨多預(yù)報偏多而大暴雨預(yù)報偏少。

    (2)GRAPES模式在2019年中雨及以上量級降水預(yù)報上改善大于SWCWARMS和ECMWF模式,ECMWF和SWCWARMS模式在這3年期間TS評分波動較小,SWCWARMS模式2019年空報較2017年和2018年顯著降低。逐月降水檢驗(yàn)顯示,3個模式在小雨和中雨預(yù)報上不相上下,GRAPES模式優(yōu)勢在2019年大雨和暴雨預(yù)報上,ECMWF模式優(yōu)勢在2017年秋季和2018年初夏大雨預(yù)報上,SWCWARMS模式大雨和暴雨預(yù)報能力介于二者之間;SWCWARMS模式2018年第三季度評分整體偏低,ECMWF模式2019年中到暴雨評分振幅相對較大,GRAPES模式從2018年第三季度開始中到暴雨降水評分顯著提高。

    (3)在川東地區(qū),ECMWF模式TS評分高于SWCWARMS和GRAPES模式,SWCWARMS模式TS評分24h和36h低于GRAPES模式而其它時次高于GRAPES模式,SWCWARMS模式ETS評分高于其它兩模式,ECMWF模式ETS評分低于區(qū)域模式;在川西地區(qū),ECMWF模式ETS評分低于其它兩模式,SWCWARMS模式部分時次TS評分高于GRAPES模式而ETS評分低于GRAPES模式;三個模式存在不同程度空報,川東空報略多于川西,其中ECMWF模式空報最多,GRAPES模式空報略少于SWCWARMS和ECMWF模式。

    猜你喜歡
    中雨實(shí)況小雨
    起落架用30CrMnSiNi2A鋼表面Cd-Ti鍍層在海水和鹽霧環(huán)境中的腐蝕行為研究
    下 雨
    淺談跨座式單軌交通建設(shè)中雨污水管線遷改工作
    無題(6)
    鄉(xiāng)村小學(xué)的愿望與現(xiàn)實(shí)——宜君鄉(xiāng)村教育實(shí)況
    天舟一號貨運(yùn)飛船發(fā)射實(shí)況掠影
    臥薪嘗膽
    可愛潮咖們的獨(dú)門彩妝實(shí)況直播
    一箭之仇
    假途滅虢
    日韩,欧美,国产一区二区三区| 18+在线观看网站| av黄色大香蕉| 国产亚洲精品第一综合不卡 | 久久久久久人妻| 国产视频内射| 亚洲国产欧美日韩在线播放| 亚洲av不卡在线观看| 中文字幕制服av| 99久久人妻综合| 乱码一卡2卡4卡精品| 亚洲国产精品国产精品| 中文字幕精品免费在线观看视频 | 美女中出高潮动态图| 中文字幕免费在线视频6| 九色亚洲精品在线播放| 成年人午夜在线观看视频| 99久国产av精品国产电影| 曰老女人黄片| 少妇熟女欧美另类| 婷婷成人精品国产| 精品久久久精品久久久| 一区二区三区免费毛片| 狂野欧美激情性xxxx在线观看| 亚洲情色 制服丝袜| 曰老女人黄片| 免费人成在线观看视频色| 久久婷婷青草| 熟妇人妻不卡中文字幕| 69精品国产乱码久久久| 精品国产一区二区三区久久久樱花| 亚洲成人手机| 免费久久久久久久精品成人欧美视频 | 午夜激情av网站| 性色av一级| 久久精品国产亚洲av天美| 美女中出高潮动态图| 91精品一卡2卡3卡4卡| 五月开心婷婷网| 成人亚洲欧美一区二区av| 午夜福利在线观看免费完整高清在| www.色视频.com| av卡一久久| 中文字幕免费在线视频6| 国产在线免费精品| 亚洲国产av新网站| 春色校园在线视频观看| 久久这里有精品视频免费| 天天操日日干夜夜撸| 97超视频在线观看视频| 99久国产av精品国产电影| 国产免费又黄又爽又色| 亚洲av中文av极速乱| 日韩强制内射视频| 国产乱来视频区| 国产不卡av网站在线观看| 日韩人妻高清精品专区| 免费大片18禁| 亚洲av在线观看美女高潮| .国产精品久久| 亚洲综合色惰| 国产精品久久久久久精品电影小说| 免费人成在线观看视频色| 91久久精品电影网| 制服人妻中文乱码| 久久久久久久大尺度免费视频| 伊人久久精品亚洲午夜| 亚洲av男天堂| 欧美日韩av久久| 一区二区日韩欧美中文字幕 | 啦啦啦啦在线视频资源| 18禁裸乳无遮挡动漫免费视频| 亚洲美女视频黄频| 女人精品久久久久毛片| 午夜激情久久久久久久| 久久久精品94久久精品| 亚洲天堂av无毛| 3wmmmm亚洲av在线观看| 亚洲欧洲国产日韩| 18禁裸乳无遮挡动漫免费视频| 伦理电影免费视频| 亚洲婷婷狠狠爱综合网| 欧美精品亚洲一区二区| 晚上一个人看的免费电影| 精品亚洲成国产av| 青青草视频在线视频观看| 午夜久久久在线观看| av一本久久久久| 天美传媒精品一区二区| 十八禁高潮呻吟视频| 国产av码专区亚洲av| 久久久久视频综合| 人妻系列 视频| 色5月婷婷丁香| 国产熟女午夜一区二区三区 | 免费少妇av软件| 三级国产精品欧美在线观看| 国产日韩一区二区三区精品不卡 | 成人毛片a级毛片在线播放| 欧美国产精品一级二级三级| 欧美 日韩 精品 国产| 性高湖久久久久久久久免费观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲五月色婷婷综合| 赤兔流量卡办理| 国产黄频视频在线观看| 久久午夜福利片| 午夜福利视频在线观看免费| 婷婷色麻豆天堂久久| 欧美日韩视频高清一区二区三区二| 99热这里只有是精品在线观看| 国产探花极品一区二区| 久久精品久久久久久噜噜老黄| 一区二区三区免费毛片| 少妇被粗大的猛进出69影院 | 国产一区二区在线观看av| 99久久精品国产国产毛片| 亚洲精品久久成人aⅴ小说 | 丝袜脚勾引网站| 免费看不卡的av| 国产一区二区三区av在线| 91久久精品国产一区二区成人| 国产日韩欧美亚洲二区| 国产黄色免费在线视频| 嫩草影院入口| 精品人妻在线不人妻| av播播在线观看一区| 观看美女的网站| 亚洲人成网站在线播| 精品国产一区二区三区久久久樱花| 日韩中文字幕视频在线看片| 韩国高清视频一区二区三区| www.av在线官网国产| 亚洲国产精品专区欧美| 亚洲五月色婷婷综合| 日本黄色日本黄色录像| 国产男人的电影天堂91| 晚上一个人看的免费电影| 精品国产露脸久久av麻豆| 高清av免费在线| 欧美xxxx性猛交bbbb| 麻豆成人av视频| 亚洲国产av新网站| 麻豆乱淫一区二区| 在线天堂最新版资源| 欧美另类一区| 亚洲av不卡在线观看| 欧美变态另类bdsm刘玥| 国产av码专区亚洲av| 永久免费av网站大全| 日日啪夜夜爽| 亚洲第一av免费看| 人妻制服诱惑在线中文字幕| 女人久久www免费人成看片| 久久热精品热| 国产一区亚洲一区在线观看| 日本爱情动作片www.在线观看| 精品人妻熟女av久视频| 中文字幕人妻熟人妻熟丝袜美| 久久鲁丝午夜福利片| 永久网站在线| 亚洲综合色惰| 啦啦啦视频在线资源免费观看| av女优亚洲男人天堂| 成年av动漫网址| 一级毛片黄色毛片免费观看视频| 日韩成人av中文字幕在线观看| 国产日韩欧美亚洲二区| 欧美精品一区二区免费开放| 热99国产精品久久久久久7| 蜜桃在线观看..| 少妇被粗大猛烈的视频| 午夜激情av网站| 成年人免费黄色播放视频| 97超视频在线观看视频| 国产无遮挡羞羞视频在线观看| 国产高清不卡午夜福利| 在线观看美女被高潮喷水网站| 亚洲图色成人| 成人手机av| 大片免费播放器 马上看| 青春草视频在线免费观看| 国产精品成人在线| 人妻少妇偷人精品九色| 99精国产麻豆久久婷婷| 亚洲精品久久午夜乱码| 国产色婷婷99| 国产伦精品一区二区三区视频9| 一本一本综合久久| 国产精品.久久久| 午夜激情久久久久久久| 亚洲国产欧美日韩在线播放| 纵有疾风起免费观看全集完整版| 精品少妇内射三级| 久久这里有精品视频免费| 国产欧美日韩综合在线一区二区| 少妇熟女欧美另类| 三级国产精品欧美在线观看| 国产av国产精品国产| 成人综合一区亚洲| 九九久久精品国产亚洲av麻豆| 精品人妻熟女av久视频| av有码第一页| 亚洲欧洲精品一区二区精品久久久 | 国产亚洲av片在线观看秒播厂| 免费高清在线观看视频在线观看| 99久国产av精品国产电影| 欧美激情国产日韩精品一区| 色哟哟·www| 国产精品国产三级国产专区5o| 人妻少妇偷人精品九色| 少妇的逼好多水| 18禁裸乳无遮挡动漫免费视频| 国产高清有码在线观看视频| 99热全是精品| 亚洲美女搞黄在线观看| 国产在线一区二区三区精| 日日摸夜夜添夜夜添av毛片| 69精品国产乱码久久久| 国内精品宾馆在线| 欧美老熟妇乱子伦牲交| 亚洲综合色网址| 欧美 日韩 精品 国产| 我的老师免费观看完整版| av国产精品久久久久影院| 在线观看一区二区三区激情| 一区二区av电影网| 51国产日韩欧美| 国产高清不卡午夜福利| 亚洲,一卡二卡三卡| av专区在线播放| 极品人妻少妇av视频| 日韩一区二区三区影片| 母亲3免费完整高清在线观看 | 肉色欧美久久久久久久蜜桃| 大片免费播放器 马上看| 国产精品不卡视频一区二区| 欧美精品国产亚洲| 亚洲精品一区蜜桃| 日韩中字成人| 免费看av在线观看网站| 高清黄色对白视频在线免费看| 一二三四中文在线观看免费高清| 午夜视频国产福利| 桃花免费在线播放| 成人二区视频| 日韩av在线免费看完整版不卡| 亚洲精品av麻豆狂野| 久久婷婷青草| 一级毛片 在线播放| 九色亚洲精品在线播放| 一级毛片黄色毛片免费观看视频| 啦啦啦中文免费视频观看日本| 欧美一级a爱片免费观看看| 欧美日韩在线观看h| 国产在视频线精品| 久久鲁丝午夜福利片| 久久久国产一区二区| 18禁动态无遮挡网站| 永久免费av网站大全| 91在线精品国自产拍蜜月| 精品久久久噜噜| 欧美日韩在线观看h| 久久国内精品自在自线图片| 女性生殖器流出的白浆| av一本久久久久| 亚洲精品一区蜜桃| 街头女战士在线观看网站| 欧美日韩国产mv在线观看视频| 在线播放无遮挡| 99久久中文字幕三级久久日本| 天堂俺去俺来也www色官网| 99国产综合亚洲精品| 激情五月婷婷亚洲| 精品久久久噜噜| 欧美bdsm另类| 99热这里只有是精品在线观看| 高清午夜精品一区二区三区| 国产成人免费无遮挡视频| 成人午夜精彩视频在线观看| av免费在线看不卡| 纯流量卡能插随身wifi吗| 亚洲欧美一区二区三区黑人 | 中国美白少妇内射xxxbb| 少妇人妻久久综合中文| 国产片特级美女逼逼视频| 久久这里有精品视频免费| 日韩制服骚丝袜av| 丝袜在线中文字幕| 亚洲精品视频女| 久久 成人 亚洲| 亚洲国产成人一精品久久久| 亚洲精品成人av观看孕妇| 亚洲欧洲精品一区二区精品久久久 | 熟女电影av网| 97超视频在线观看视频| a 毛片基地| 欧美老熟妇乱子伦牲交| 亚洲熟女精品中文字幕| av国产精品久久久久影院| 伊人久久国产一区二区| 天堂中文最新版在线下载| 精品一区在线观看国产| 国产精品国产三级国产专区5o| 国产一区二区三区综合在线观看 | 夜夜骑夜夜射夜夜干| 国产成人freesex在线| 成人18禁高潮啪啪吃奶动态图 | av免费在线看不卡| 乱码一卡2卡4卡精品| 搡女人真爽免费视频火全软件| 丝袜在线中文字幕| 欧美精品一区二区免费开放| 免费看av在线观看网站| 国产精品一国产av| 丝袜在线中文字幕| 3wmmmm亚洲av在线观看| 中文字幕免费在线视频6| av免费观看日本| 欧美丝袜亚洲另类| 亚洲色图 男人天堂 中文字幕 | 久久精品国产a三级三级三级| 精品久久久久久久久亚洲| 韩国高清视频一区二区三区| 街头女战士在线观看网站| 少妇熟女欧美另类| 麻豆成人av视频| 日韩三级伦理在线观看| tube8黄色片| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av不卡在线观看| 欧美日韩视频精品一区| 欧美 日韩 精品 国产| 国产一区二区三区av在线| 国产精品国产三级国产专区5o| 久久久亚洲精品成人影院| 下体分泌物呈黄色| 男女边摸边吃奶| 精品一区二区三区视频在线| 国产永久视频网站| 少妇熟女欧美另类| 成人手机av| 国产亚洲一区二区精品| 亚洲精品日韩av片在线观看| 欧美日韩亚洲高清精品| 最黄视频免费看| 婷婷成人精品国产| 精品国产乱码久久久久久小说| 嘟嘟电影网在线观看| 国产午夜精品久久久久久一区二区三区| 久久精品国产亚洲网站| 一区二区三区精品91| 国产乱来视频区| 亚洲精品久久久久久婷婷小说| 狂野欧美白嫩少妇大欣赏| 欧美精品人与动牲交sv欧美| 亚洲精品第二区| 午夜日本视频在线| 99热这里只有是精品在线观看| 国产极品天堂在线| 久久精品久久久久久久性| 成人手机av| 另类精品久久| 99热国产这里只有精品6| 在线亚洲精品国产二区图片欧美 | tube8黄色片| 免费不卡的大黄色大毛片视频在线观看| 乱码一卡2卡4卡精品| 色吧在线观看| 欧美成人午夜免费资源| 夜夜看夜夜爽夜夜摸| 久久人妻熟女aⅴ| 免费高清在线观看视频在线观看| 大香蕉97超碰在线| 亚洲av成人精品一区久久| 久久久久久久久久久丰满| 国产男女超爽视频在线观看| 一边亲一边摸免费视频| 女的被弄到高潮叫床怎么办| 一级毛片黄色毛片免费观看视频| 国产片内射在线| 亚洲美女搞黄在线观看| 男的添女的下面高潮视频| 久久青草综合色| 又黄又爽又刺激的免费视频.| 日韩av不卡免费在线播放| 国语对白做爰xxxⅹ性视频网站| a级片在线免费高清观看视频| 草草在线视频免费看| 亚洲精品第二区| 夜夜看夜夜爽夜夜摸| 黑人欧美特级aaaaaa片| 亚洲精华国产精华液的使用体验| 视频中文字幕在线观看| 看十八女毛片水多多多| 夜夜骑夜夜射夜夜干| 欧美xxⅹ黑人| 亚洲欧美一区二区三区国产| 一区二区三区精品91| 婷婷色av中文字幕| 日韩精品有码人妻一区| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品久久成人aⅴ小说 | 久久精品久久久久久久性| 成人黄色视频免费在线看| 亚洲国产成人一精品久久久| 99热网站在线观看| 免费人成在线观看视频色| 激情五月婷婷亚洲| 精品人妻熟女毛片av久久网站| 国精品久久久久久国模美| videosex国产| 国产精品久久久久成人av| 亚洲精品aⅴ在线观看| 免费黄色在线免费观看| 日本午夜av视频| 亚洲精品乱码久久久v下载方式| 桃花免费在线播放| 亚洲精品久久成人aⅴ小说 | 80岁老熟妇乱子伦牲交| 国产综合精华液| 免费大片黄手机在线观看| 国产精品久久久久久av不卡| 日本免费在线观看一区| 在现免费观看毛片| 蜜臀久久99精品久久宅男| 亚洲av综合色区一区| 26uuu在线亚洲综合色| 精品人妻一区二区三区麻豆| 国产成人精品婷婷| 色5月婷婷丁香| 国产av国产精品国产| 黑丝袜美女国产一区| 老司机亚洲免费影院| 亚洲av中文av极速乱| 成人二区视频| 色婷婷av一区二区三区视频| 一级毛片aaaaaa免费看小| 亚洲精品国产色婷婷电影| 精品少妇黑人巨大在线播放| 午夜91福利影院| 我的老师免费观看完整版| 少妇熟女欧美另类| 成人国产麻豆网| 国产高清三级在线| 丰满乱子伦码专区| 免费大片黄手机在线观看| 性色avwww在线观看| 国产成人精品无人区| 国产日韩欧美视频二区| 成人国产麻豆网| 国产一区有黄有色的免费视频| 欧美日韩综合久久久久久| 久久久久久久精品精品| 99国产综合亚洲精品| 下体分泌物呈黄色| 日本色播在线视频| 丝袜喷水一区| 国产高清国产精品国产三级| 伦理电影大哥的女人| 国产精品久久久久成人av| 日韩av在线免费看完整版不卡| 亚洲精品456在线播放app| 国产伦精品一区二区三区视频9| 日本欧美视频一区| 久久午夜综合久久蜜桃| 国产成人精品在线电影| 国产精品久久久久久精品电影小说| 免费高清在线观看视频在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产精品成人久久小说| 中文欧美无线码| 久久久久久久久久久丰满| 一级片'在线观看视频| 色5月婷婷丁香| 大片免费播放器 马上看| 永久免费av网站大全| 性色av一级| 精品久久久久久电影网| 久久久精品免费免费高清| 国产综合精华液| 成人亚洲欧美一区二区av| 简卡轻食公司| 少妇被粗大猛烈的视频| 国产色婷婷99| 亚洲精品久久午夜乱码| 久久人妻熟女aⅴ| 三级国产精品片| 亚洲怡红院男人天堂| 一级毛片aaaaaa免费看小| 精品少妇久久久久久888优播| 大陆偷拍与自拍| 国产黄色视频一区二区在线观看| 国产欧美日韩综合在线一区二区| 男人添女人高潮全过程视频| 国产黄片视频在线免费观看| 日本-黄色视频高清免费观看| 麻豆精品久久久久久蜜桃| 日韩制服骚丝袜av| 少妇的逼水好多| 国产免费视频播放在线视频| 中文天堂在线官网| 中文字幕亚洲精品专区| 国产亚洲欧美精品永久| 99re6热这里在线精品视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品视频女| 亚洲美女搞黄在线观看| 欧美另类一区| 亚洲成色77777| 欧美精品亚洲一区二区| 亚洲欧洲国产日韩| 久久久精品免费免费高清| 久久 成人 亚洲| 男的添女的下面高潮视频| 亚洲国产成人一精品久久久| 2018国产大陆天天弄谢| 欧美激情极品国产一区二区三区 | 男女无遮挡免费网站观看| 最近中文字幕高清免费大全6| 亚洲第一区二区三区不卡| 欧美丝袜亚洲另类| 亚洲精品色激情综合| 欧美+日韩+精品| 最近手机中文字幕大全| 大香蕉久久成人网| 一个人免费看片子| 欧美日韩视频高清一区二区三区二| 在线天堂最新版资源| 国产精品一二三区在线看| 99久久中文字幕三级久久日本| 在线播放无遮挡| 一级毛片我不卡| 久久99精品国语久久久| av在线老鸭窝| 欧美日韩视频高清一区二区三区二| 肉色欧美久久久久久久蜜桃| 亚洲精品中文字幕在线视频| 欧美精品国产亚洲| 国产精品一区二区三区四区免费观看| 久久久久精品性色| 亚洲欧美成人综合另类久久久| 99久久精品一区二区三区| 在线观看免费高清a一片| 22中文网久久字幕| 成年av动漫网址| 一级毛片电影观看| 成人亚洲欧美一区二区av| 国产精品国产三级国产专区5o| 99久久中文字幕三级久久日本| 午夜福利,免费看| 99久国产av精品国产电影| 欧美人与性动交α欧美精品济南到 | 亚洲欧洲日产国产| 国产成人精品在线电影| 国产精品国产三级国产专区5o| 十分钟在线观看高清视频www| 亚洲性久久影院| 免费日韩欧美在线观看| 国产熟女午夜一区二区三区 | 久久精品久久精品一区二区三区| 交换朋友夫妻互换小说| av线在线观看网站| 在线观看国产h片| 狠狠婷婷综合久久久久久88av| 啦啦啦啦在线视频资源| 大片电影免费在线观看免费| 久久久久久久久久久免费av| 秋霞在线观看毛片| 亚洲一级一片aⅴ在线观看| 91精品三级在线观看| 久久久久精品性色| www.色视频.com| 伦精品一区二区三区| 女性被躁到高潮视频| 亚洲欧洲国产日韩| 中国美白少妇内射xxxbb| 在现免费观看毛片| 亚洲av成人精品一二三区| 国产不卡av网站在线观看| 免费观看的影片在线观看| 国产熟女午夜一区二区三区 | 国产成人免费无遮挡视频| 午夜日本视频在线| 又黄又爽又刺激的免费视频.| 国产av国产精品国产| 只有这里有精品99| 亚洲第一av免费看| 国产成人a∨麻豆精品| 欧美日韩成人在线一区二区| 国产免费又黄又爽又色| 午夜福利视频精品| 只有这里有精品99| 美女中出高潮动态图| 一区二区三区免费毛片| 国产 一区精品| 国产日韩欧美在线精品| 国产成人a∨麻豆精品| 国产精品偷伦视频观看了| 男女无遮挡免费网站观看| av国产久精品久网站免费入址| 国产一区二区三区av在线| 精品亚洲成国产av| 精品一区在线观看国产| 国产成人精品福利久久| 99久国产av精品国产电影| 2022亚洲国产成人精品| 欧美日韩亚洲高清精品| 亚洲精品色激情综合| 免费看不卡的av| 亚洲av电影在线观看一区二区三区| 欧美人与善性xxx| 久久午夜福利片| 视频在线观看一区二区三区|