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

    不同方法對青海2020年強降水模式產(chǎn)品預(yù)報性能的檢驗對比

    2022-05-17 06:25:32沈曉燕申燕玲杜華禮顏玉倩
    干旱氣象 2022年2期

    沈曉燕,申燕玲,權(quán) 晨,杜華禮,顏玉倩

    (1.青海省氣象科學(xué)研究所,青海 西寧 810001;2.青海省防災(zāi)減災(zāi)重點實驗室,青海 西寧 810001)

    引 言

    隨著數(shù)值預(yù)報技術(shù)的不斷發(fā)展及模式分辨率的逐步提高,數(shù)值模式已成為天氣預(yù)報的重要手段。預(yù)報檢驗是數(shù)值模式發(fā)展及應(yīng)用的重要組成部分,檢驗?zāi)軌虮嫖霾煌J街g的客觀差異[1],不僅可以反饋給模式開發(fā)者以提高其預(yù)報性能,還可以為模式應(yīng)用者提供有價值的參考信息。

    以往模式預(yù)報檢驗分析大多采用點對點的傳統(tǒng)檢驗方法[2-8],如TS評分、空報率、漏報率以及HSS評分、ORSS評分、HK評分、SEDS評分[9-10]等。傳統(tǒng)檢驗方法定量且易于理解,但預(yù)報時間和空間上的細微偏差均會導(dǎo)致漏報或空報,從而引起總體評分較低,且缺乏對預(yù)報和觀測之間誤差來源的詳細分析,有時還與模式預(yù)報的主觀評價不一致。因此,空間檢驗逐漸備受關(guān)注[11]??臻g檢驗方法可分為4種類型:鄰域法、尺度分離法、特征法、場變形法[12-15]。不同方法各有優(yōu)缺點,鄰域法也稱“迷糊法”,通過對高分辨率模式預(yù)報結(jié)果濾波,獲取不同空間尺度評分來體現(xiàn)模式預(yù)報能力[16],它通過放寬傳統(tǒng)檢驗中觀測與預(yù)報一一對應(yīng)的要求,給予與觀測場空間接近的更高預(yù)報評分,避免了將此類預(yù)報判定為無效預(yù)報[17];尺度分離法是分離出與尺度相關(guān)的誤差,需要通過小波分析等空間濾波器進行檢驗評估;特征法是檢驗降水系統(tǒng)的宏觀尺度特征,而場變形法則是對預(yù)報場進行形變以便更好地與觀測場進行匹配,與特征法相比它還能夠識別出預(yù)報與實況嚴(yán)重不符的情況,但計算過程復(fù)雜。FSS(fraction skill score)是目前應(yīng)用較為普遍的空間鄰域檢驗方法,若檢驗?zāi)康氖菫榱嗽\斷預(yù)報誤差來源,則特征法更為合適,而面向?qū)ο蟮臋z驗方法(method of object-based diagnostic evaluation,MODE)是最具代表性的特征法,它通過分析預(yù)報與觀測對象的空間特性,有效診斷出對應(yīng)組合對象的位置、大小等屬性差異[18-20]。通過對比TS(threat score)評分、FSS評分、MODE方法在模式降水預(yù)報性能檢驗中的應(yīng)用,發(fā)現(xiàn)傳統(tǒng)檢驗方法雖然有一定的缺陷,但能夠得出模式預(yù)報的整體性能,故而仍然具有一定的適用性[21-26];在降水量級較小情況下FSS評分與TS評分的評估效果相差不大,但在降水量級較大時TS評分不穩(wěn)定,而FSS評分依舊可以區(qū)分不同條件下模式預(yù)報性能[27],且可以判別預(yù)報效果較好的空間尺度(鄰域半徑);MODE方法通過計算降水對象的多種位置、形狀屬性,從多視角展現(xiàn)模式預(yù)報能力[28]。

    青海位于青藏高原東北部,強降水次數(shù)雖然較少,但因局地性強,在地形復(fù)雜特殊、生態(tài)環(huán)境脆弱條件下極易引發(fā)山洪、泥石流等災(zāi)害,因此提升模式降水預(yù)報準(zhǔn)確率對當(dāng)?shù)胤罏?zāi)減災(zāi)、人類生命財產(chǎn)安全意義重大,而預(yù)報檢驗是模式改進的重要前提。目前,全球及區(qū)域模式對強降水預(yù)報能力較弱,傳統(tǒng)點對點檢驗效果較差,無法體現(xiàn)模式性能差異,而空間檢驗則是傳統(tǒng)點對點檢驗的有益補充。為此,本研究針對青海省2020年7—8月強降水天氣個例,選取多種模式降水預(yù)報產(chǎn)品,采用傳統(tǒng)檢驗方法、鄰域法及面向?qū)ο蠓椒?,對比檢驗各模式在青海地區(qū)強降水預(yù)報中的性能,以期為不同模式在該地區(qū)的應(yīng)用及改進提供一定參考。

    1 資料及處理

    按照24 h降水量達到大雨及以上量級的站點超過3個的原則,選出2020年7—8月青海省[包括海西蒙古族藏族自治州(簡稱“海西州”)、海北藏族自治州(簡稱“海北州”)、海南藏族自治州(簡稱“海南州”)、玉樹藏族自治州(簡稱“玉樹州”)、黃南藏族自治州(簡稱“黃南州”)、果洛藏族自治州(簡稱“果洛州”)及西寧市、海東市8個行政區(qū)域]10次強降水天氣個例(表1),強降水落區(qū)主要集中在祁連山區(qū)、東部農(nóng)區(qū)、玉樹州東部、果洛州東部及海南州、黃南州地區(qū)。

    所用資料包括:(1)歐洲中期天氣預(yù)報中心(European Center of Medium-range Weather Forecasts,ECMWF)(水平分辨率12.5 km)和美國國家環(huán)境預(yù)報中心(National Center for Environmental Prediction,NCEP)(水平分辨率25 km)全球模式、中國氣象局全球同化預(yù)報系統(tǒng)(China Meteorological Administration global assimilation forecast system,CMA-GFS)(水平分辨率25 km)及GRAPES(global and regional assimilation and prediction system)區(qū)域中尺度數(shù)值預(yù)報系統(tǒng)(GRAPES regional mesoscale numerical prediction system,GRAPES-Meso)(水平分辨率9 km)24 h降水量預(yù)報產(chǎn)品;(2)青海省232個氣象站24 h降水量觀測資料,其中52個國家級地面觀測站,180個國家地面天氣站,站點分布如圖1所示;(3)中國氣象局國家氣象信息中心制作的青海區(qū)域地面-衛(wèi)星-雷達三源融合的實時降水格點分析產(chǎn)品(CMA multi-source precipitation analysis system,CMPAS),水平分辨率為5 km。文中的青海省及各市(州)行政邊界是基于國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號為GS(2019)1696號的標(biāo)準(zhǔn)地圖制作,底圖無修改。

    表1 2020年青海7—8月強降水天氣個例信息Tab.1 Information of heavy precipitation weather cases in Qinghai from July to August in 2020

    圖1 青海省地形高度(彩色區(qū),單位:m)與氣象站點(黑點)分布Fig.1 Distribution of topography elevation (colour areas, Unit: m) and meteorological stations (black dots) in Qinghai Province

    傳統(tǒng)檢驗過程中采用最鄰近插值法將模式網(wǎng)格數(shù)據(jù)插值到232個站點上,而空間檢驗過程中則采用雙線性插值法將模式網(wǎng)格數(shù)據(jù)插值到CMPAS產(chǎn)品網(wǎng)格上,空間檢驗范圍為青海省區(qū)域(89.2°E—102.95°E、31.4°N—39.4°N)。另外,根據(jù)24 h累計降水量,將降水等級劃分為小雨(0.1~9.9 mm)、中雨(10.0~24.9 mm)、大雨(25.0~49.9 mm)和暴雨(≥50.0 mm)。

    2 檢驗方法

    2.1 傳統(tǒng)檢驗法

    傳統(tǒng)檢驗選用TS評分、空報率(false alarm ratio,F(xiàn)AR)、漏報率(missing ratio,MR)、預(yù)報偏差(Bias)作為檢驗指標(biāo),計算公式如下:

    (1)

    (2)

    (3)

    (4)

    式中:k為降水等級;NAk為k等級降水預(yù)報正確的站數(shù);NBk為k等級降水空報的站數(shù);NCk為k等級降水漏報的站數(shù)。

    2.2 空間檢驗法

    FSS評分是目前應(yīng)用較為普遍的空間鄰域檢驗方法,它是基于FBS評分(fraction Brier score)的變形,利用方差技巧評分構(gòu)建思路,用來檢驗預(yù)報降水概率[27],其計算公式如下:

    (5)

    (6)

    式中:Pf、Po分別為單個鄰域窗內(nèi)預(yù)報和觀測降水發(fā)生的概率(可類比于鄰域窗內(nèi)降水的概率);N為檢驗區(qū)域內(nèi)鄰域窗個數(shù)。

    MODE檢驗方法是將觀測場和預(yù)報場進行平滑,設(shè)定雨量閾值,選取關(guān)注的降雨面積最小值提取對象。然后,通過計算對象質(zhì)心距離、目標(biāo)面積、最小邊界距離等屬性,將相似的觀測場和預(yù)報場目標(biāo)進行匹配,如果觀測場和預(yù)報場目標(biāo)存在多對多的情況,則將多對多的目標(biāo)進行合并。最后,通過綜合計算觀測和預(yù)報對應(yīng)對象的屬性,得出預(yù)報對象的整體相似度,其公式[28]如下:

    (7)

    式中:ai、wi分別表示對象第i個屬性及其權(quán)重系數(shù);Ci表示對象第i個屬性的置信水平;I(ai)表示對象第i個屬性的收益函數(shù),是ai的函數(shù)。通常認(rèn)為,目標(biāo)相似度達到0.7時,模式對雨帶的預(yù)報較為準(zhǔn)確。

    在計算FSS評分時使用了MET V9.1.1(model evaluation tools)。MET是由美國數(shù)值預(yù)報發(fā)展試驗中心研發(fā),主要針對WRF(weather research and forecasting)中尺度模式(也可應(yīng)用于其他模式)提供多種檢驗方法,其中FSS評分包含于基于格點-格點的鄰域空間檢驗方法。MODE的計算采用基于Python的Meteva程序庫,該程序庫由中國氣象局國家氣象中心預(yù)報技術(shù)研發(fā)室檢驗科研發(fā),同樣提供了多種檢驗算法。

    3 結(jié)果與分析

    3.1 整體檢驗評分

    圖2是2020年7—8月青海10次強降水天氣個例的傳統(tǒng)點對點檢驗TS評分、漏報率、空報率及偏差。可以看出,各模式小雨及以上量級預(yù)報TS評分相差不大,平均值介于0.85~0.87之間,ECMWF略優(yōu),CMA-GFS略差;GRAPES-Meso漏報較多,空報較少,而其他全球模式空報較多。針對中雨及以上量級,NCEP和ECMWF預(yù)報較好,TS評分分別為0.36、0.35,GRAPES-Meso預(yù)報相對較差,TS評分為0.29。其中,NCEP空報率、漏報率均較低,GRAPES-Meso漏報較多,CMA-GFS空報較多。針對大雨及以上量級,GRAPES-Meso預(yù)報相對較好,但TS評分僅為0.06,所有模式空報和漏報均較高,其中NCEP和GRAPES-Meso相對較優(yōu)。從預(yù)報偏差來看,小雨及以上量級的預(yù)報偏差較小,介于1.03~1.14之間,而中雨及以上量級和大雨及以上量級的模式預(yù)報均偏多,且CMA-GFS預(yù)報偏差均最大(1.72、2.85),預(yù)報偏差最小的分別是GRAPES-Meso(1.18)、ECMWF(1.38)。

    圖2 2020年7—8月青海強降水天氣個例不同模式降水預(yù)報的平均TS評分(a)、漏報率(b)、空報率(c)及偏差(d)Fig.2 The average TS score (a), MR (b), FAR (c), Bias (d) of precipitation forecast by different models for heavy precipitation weather cases in Qinghai Province from July to August 2020

    鄰域檢驗法是計算檢驗區(qū)域內(nèi)所有空間鄰域集合的誤差矩陣。鄰域是以每個格點為中心的網(wǎng)格窗口,窗口可以是圓形或正方形,不同鄰域窗對檢驗結(jié)果影響不大[17]。為簡便計算,本文選擇正方形鄰域窗,鄰域半徑為3×3、5×5、…、21×21網(wǎng)格窗口。由于鄰域檢驗使用的CMPAS降水格點分析產(chǎn)品空間分辨率為5 km,因此計算的鄰域半徑相應(yīng)為15×15、25×25、…、105×105 km。從圖3看出,隨著鄰域空間尺度的增大,整體上各量級降水的FSS評分逐漸增加;不同量級降水,各模式預(yù)報的FSS評分優(yōu)劣存在差異,小雨及以上量級FSS評分自GRAPES-Meso、NCEP、ECMWF、CMA-GFS依次減小,中雨及以上量級自NCEP、ECMWF、CMA-GFS、GRAPES-Meso依次減小,其結(jié)果與TS評分一致,而大雨及以上量級則自ECMWF、CMA-GFS、GRAPES-Meso、NCEP依次降低。

    分別針對中雨和大雨及以上量級進行MODE檢驗,并對比目標(biāo)屬性(表2)??梢钥闯觯椭杏昙耙陨狭考夘A(yù)報而言,NCEP整體相似度相對較好,ECMWF次之,GRAPES-Meso最差,這與TS評分和FSS評分一致。從10次個例的相關(guān)屬性(圖略)來看,NCEP質(zhì)心距離僅次于GRAPES-Meso,總體平均面積比(觀測和預(yù)報中較小者與較大者之比)及交叉面積最大,ECMWF各目標(biāo)屬性僅次于NCEP,而CMA-GFS交叉面積小,質(zhì)心距離整體偏大;GRAPES-Meso預(yù)報的7月至8月初幾次過程整體相似度低,交叉面積小,質(zhì)心距離遠,而對8月中下旬海東地區(qū)過程表現(xiàn)較好。從預(yù)報與觀測目標(biāo)質(zhì)心距離及方位角[圖4(a)]來看,不同模式預(yù)報整體較觀測偏西,這與青海中南部中雨及以上雨量預(yù)報偏多有關(guān)系。

    圖3 2020年7—8月青海強降水個例不同模式預(yù)報的平均FSS評分Fig.3 The average FSS score of precipitation forecast by different models for heavy precipitation cases in Qinghai Province from July to August 2020

    表2 2020年7—8月青海強降水天氣個例MODE檢驗的預(yù)報與觀測對象平均屬性Tab.2 The average attributes of forecast objects corresponding to observation for MODE test of heavy precipitation weather cases in Qinghai Province from July to August 2020

    圖4 2020年7—8月青海強降水天氣個例MODE檢驗的中雨及以上量級(a)和大雨及以上量級(b)預(yù)報與觀測對象質(zhì)心距離(單位:km)及方位角(4位數(shù)字為強降水發(fā)生日期)Fig.4 The centroid distance (Unit: km) and azimuth of moderate rain and above (a) and heavy rain and above (b) forecast objects corresponding to observation for MODE test of heavy precipitation weather cases in Qinghai Province from July to August 2020(four digits for the occurrence date of heavy rainfall)

    對于大雨及以上量級預(yù)報(表2),ECMWF、CMA-GFS、GRAPES-Meso整體相似度相差不大,而NCEP整體相似度明顯偏低;NCEP預(yù)報與觀測目標(biāo)質(zhì)心距離最大、面積比和交叉面積最小,F(xiàn)SS評分最低;ECMWF質(zhì)心距離最小、交叉面積最大,CMA-GFS面積比最大,兩種模式FSS評分較高;GRAPES-Meso質(zhì)心距離、面積比、交叉面積均相對表現(xiàn)較好。另外,ECMWF、GRAPES-Meso的大雨及以上量級中心預(yù)報均較觀測偏北,而NCEP預(yù)報的2次大雨及以上量級中心較觀測偏南[圖4(b)]。

    整體來看,不同模式對小雨及以上量級預(yù)報評分均較高,模式之間差異不大,但不同檢驗方法下評分最高的模式略有不同;4種模式的中雨及以上量級預(yù)報在不同檢驗方法下優(yōu)劣程度表現(xiàn)一致,均自NCEP、ECMWF、CMA-GFS、GRAPES-Meso依次降低,其中MODE檢驗可以分析預(yù)報與觀測對應(yīng)對象的質(zhì)心距離、方位角、面積比、交叉面積等屬性,進而獲得預(yù)報與觀測間的誤差來源;對于大雨及以上量級,傳統(tǒng)點對點檢驗評分均較低,很難區(qū)別模式間預(yù)報能力,F(xiàn)SS評分雖有效提高了模式間預(yù)報差異性的評估能力,但MODE檢驗則提供了模式應(yīng)用參考價值。

    3.2 典型個例檢驗評分

    2020年7—8月青海省大到暴雨超過10站的個例有3次(表1),均出現(xiàn)在8月,分別為5日08:00至6日08:00(簡稱“case1”)、22日08:00至23日08:00(簡稱“case2”)、28日20:00至29日20:00(簡稱“case3”),這3次個例大雨站數(shù)較多且較為集中。當(dāng)大雨站數(shù)較多且較為集中時,MODE檢驗的目標(biāo)對象屬性更為明顯。下面詳細對比分析這3次個例不同檢驗方法結(jié)果。

    從圖5看出,CMPAS降水格點分析產(chǎn)品上,case1和case2大雨及以上量級站點主要出現(xiàn)在海東市、海南州、黃南州,232站中分別有11、18站達到大雨量級,2、1站達到暴雨量級,最大降水量分別出現(xiàn)在海南州南部同德秀麻鄉(xiāng)(60.8 mm)及海南州北部共和縣石乃亥鄉(xiāng)(94.0 mm);case3大雨及以上量級站點主要出現(xiàn)在祁連山區(qū)、海東市,35站達到大雨量級,9站達到暴雨量級,最大降水量出現(xiàn)在湟源(81.5 mm)。在case1中,針對中雨及以上量級,ECMWF對青海東部降水預(yù)報較好,而NCEP、CMA-GFS預(yù)報的青海東部降水中心位置偏北,GRAPES-Meso預(yù)報的青海東部降水范圍偏小;除ECMWF外,其余模式對玉樹州北部均有空報。在case2、case3中,不同模式對青海東部降水預(yù)報相對較好,但對青海南部均有空報。

    綜上可見,就中雨及以上量級而言,不同模式對case1預(yù)報效果存在明顯差異,NCEP、ECMWF、CMA-GFS、GRAPES-Meso預(yù)報的TS評分依次為0.27、0.26、0.09、0.02[圖6(a)],最高與最低評分相差0.25。同樣,各模式預(yù)報的空間檢驗結(jié)果也差異較明顯,F(xiàn)SS評分[圖7(a)]與MODE檢驗結(jié)果[圖8(a)]一致,自ECMWF、NCEP、CMA-GFS、GRAPES-Meso依次降低,MODE檢驗的整體相似度依次為0.91、0.89、0.73、0.54。另外,從MODE檢驗的預(yù)報與觀測對象屬性來看,各模式對case1預(yù)報對象的質(zhì)心方位角均較觀測偏西北[圖4(a)],其中ECMWF預(yù)報與觀測對象質(zhì)心距離僅次于GRAPES-Meso,面積比、交叉面積最大;GRAPES-Meso預(yù)報與觀測對象質(zhì)心距離、面積比、交叉面積均最小[圖8(a)],說明預(yù)報較實況范圍偏小。

    case2和case3的中雨及以上量級預(yù)報TS評分明顯優(yōu)于case1,但在case2和case3中不同模式預(yù)報結(jié)果差異不明顯,最高與最低TS評分分別相差0.11、0.05[圖6(a)],且不同檢驗方法下各模式表現(xiàn)略有差異。在case2中,ECMWF、NCEP、GRAPES-Meso、CMA-GFS預(yù)報的TS評分分別為0.61、0.59、0.54、0.50,F(xiàn)SS評分自CMA-GFS、ECMWF、NCEP、GRAPES-Meso依次降低[圖7(a)],MODE整體相似度GRAPES-Meso、NCEP、CMA-GFS、ECMWF依次為0.83、0.73、0.72、0.70[圖8(a)]。從MODE檢驗來看,GRAPES-Meso質(zhì)心距離最小、面積比和交叉面積最大,整體相似度高,而ECMWF面積比偏小,整體相似度最低。由于CMA-GFS降水空間分布形態(tài)更為接近實況,故而FSS評分較高。另外,各模式對青海南部空報較多,故而預(yù)報對象的質(zhì)心方位角普遍較觀測偏西[圖4(a)]。在case3中,傳統(tǒng)TS評分GRAPES-Meso(0.55)表現(xiàn)較好[圖6(a)],而FSS評分及MODE檢驗整體相似度NCEP表現(xiàn)略好。從預(yù)報相對觀測的方位角來看,各模式均偏西南,與青海南部地區(qū)空報有關(guān),質(zhì)心距離CMA-GFS最小,NCEP次之,但NCEP的面積比和交叉面積最大,故而MODE檢驗整體相似度較高[圖8(a)]。

    當(dāng)各模式大雨及以上量級預(yù)報的傳統(tǒng)檢驗TS評分均為0時,空間檢驗則是傳統(tǒng)檢驗的有益補充。在case1中,各模式大雨及以上量級預(yù)報的TS評分均為0,模式差異無法體現(xiàn)[圖6(b)]。但從模式預(yù)報空間分布(圖5)來看,case1中ECMWF預(yù)報較好,NCEP落區(qū)位于玉樹州北部,為空報,而CMA-GFS落區(qū)位置偏東北,但范圍與實況接近,GRAPES-Meso漏報。空間檢驗顯示,ECMWF與CMA-GFS預(yù)報結(jié)果較優(yōu),兩種模式預(yù)報對象均位于實況東北部;ECMWF質(zhì)心距離較小、交叉面積較大,但面積比遠小于CMA-GFS,故ECMWF的FSS評分高于CMA-GFS[圖7(b)],而整體相似度小于CMA-GFS[圖8(b)]。

    圖5 2020年8月青海省3次暴雨天氣個例的降水分析產(chǎn)品和模式預(yù)報產(chǎn)品(單位:mm)(從上至下依次是CMPAS和ECMWF、NCEP、CMA-GFS、GRAPES-Meso,從左至右依次為case1、case2、case3)Fig.5 The precipitation analysis product and forecast products of model for three rainstorm weather cases in Qinghai Province in August 2020 (Unit: mm)(The precipitation products are CMPAS, ECMWF, NCEP, CMA-GFS and GRAPES-Meso from top to bottom in turn; the cases are case1, case2 and case3 from left to right in turn)

    圖6 2020年8月青海省3次暴雨個例不同模式中雨及以上量級(a)和大雨及以上量級(b)預(yù)報的TS評分Fig.6 TS score of moderate rain and above (a) and heavy rain and above (b) forecast by different models for three rainstorm cases in Qinghai Province in August 2020

    從圖5看出,case2,不同模式均報出了青海東部大雨及以上量級降水中心;case3,除NCEP外,其余模式均報出了青海東部大雨及以上量級降水中心。當(dāng)模式預(yù)報較好時,不同檢驗方法下模式評分均較高。在case2中,CMA-GFS在不同檢驗方法下均表現(xiàn)最優(yōu),其TS評分為0.24,遠大于其余模式,且面積比、交叉面積最大,整體相似度高(0.66);ECMWF、GRAPES-Meso預(yù)報對象質(zhì)心均位于實況北部,而NCEP、CMA-GFS預(yù)報值均位于實況西部,且NCEP略偏南。具體來看,CMA-GFS對青海東部及東南部預(yù)報均較好,但質(zhì)心距離較大、方位偏西,對青海南部地區(qū)有空報;ECMWF、GRAPES-Meso對青海東部預(yù)報較好,但對東南部存在漏報,而NCEP預(yù)報位置偏南。當(dāng)模式未預(yù)報出降水中心或與實況偏離較遠時,不同檢驗方法下模式評分均較差。在case3中,除NCEP表現(xiàn)較差外,其余模式TS評分相對較高且相差不大,GRAPES-Meso、CMA-GFS、ECMWF分別為0.27、0.24、0.22[圖6(b)]。不同檢驗方法下3種較優(yōu)模式也略有差異,F(xiàn)SS評分自ECMWF、GRAPES-Meso、CMA-GFS依次降低[圖7(b)],整體相似度自CMA-GFS(0.83)、GRAPES-Meso(0.79)、ECMWF(0.65)依次減小[圖8(b)]。由于CMA-GFS預(yù)報與觀測對應(yīng)目標(biāo)面積較為接近、交叉面積大,整體相似度最高;ECMWF空間分布形態(tài)與實況最為接近,質(zhì)心距離最小,F(xiàn)SS評分最高,但在祁連山區(qū)預(yù)報面積小,整體相似度不高;NCEP三種評分結(jié)果均最差。

    結(jié)合環(huán)流形勢(圖略)來看,case1降水區(qū)位于副熱帶高壓邊緣西南氣流中,水汽條件好。不同模式傳統(tǒng)檢驗TS評分差異較大時,空間檢驗結(jié)果(MODE和FSS評分)同樣差異較大, CMA-GFS和GRAPES-Meso表現(xiàn)均較差,不同檢驗方法得出的結(jié)論較為一致。case2、case3,青海中東部存在高空槽,低層有切變線,青海東部有地面輻合線。不同模式中雨及以上量級預(yù)報效果均較好,模式之間差異不大,但不同檢驗方法下表現(xiàn)較優(yōu)的模式略有不同,觀測與預(yù)報對應(yīng)對象在面積接近、交叉面積大情況下,MODE整體相似度較高,而在質(zhì)心距離接近、空間分布形態(tài)相似時,F(xiàn)SS評分較高;CMA-GFS和GRAPES-Meso在大雨及以上量級預(yù)報中表現(xiàn)較好。

    4 結(jié) 論

    (1)2020年7—8月青海強降水天氣個例,ECMWF、NCEP、GRAPES-Meso、CMA-GFS模式對小雨及以上量級預(yù)報的TS評分均較高(超過0.84),預(yù)報偏差略偏高(介于1.03~1.14之間);對中雨和大雨及以上量級預(yù)報均偏多,中雨及以上量級CMA-GFS預(yù)報偏差最大為1.72,GRAPES-Meso偏差最小為1.18,而大雨及以上量級ECMWF偏差較小(1.38),其余模式偏多明顯(介于2.12~2.85之間)。

    (2)2020年7—8月青海強降水天氣個例,小雨及以上量級,不同模式間相差不大,不同檢驗方法下表現(xiàn)較優(yōu)的模式略有差異;中雨及以上量級,不同模式間差異較為明顯,不同檢驗方法下各模式優(yōu)劣表現(xiàn)一致;大雨及以上量級,各模式預(yù)報均較差,TS評分差異不大,F(xiàn)SS評分能夠有效提高模式差異性的表現(xiàn)能力,MODE方法則提供降水對象的屬性表現(xiàn)。另外,中雨及以上量級預(yù)報較觀測普遍偏西,質(zhì)心距離介于1.69~1.93 km之間,與青海南部空報較多有關(guān),而大雨及以上量級預(yù)報普遍偏北,質(zhì)心距離介于0.62~1.12 km之間。

    (3)當(dāng)模式預(yù)報效果較好時,不同檢驗方法下模式評分均較高;當(dāng)模式未預(yù)報出降水中心或與實況偏離較遠時,不同檢驗方法下模式評分均較差。

    總體來看,4個模式產(chǎn)品中ECMWF、CMA_GFS、GRAPES_Meso對強降水有較好的參考價值,ECMWF預(yù)報與實況質(zhì)心距離最為接近,CMA_GFS預(yù)報與實況面積比最大,ECMWF、GRAPES_Meso預(yù)報較觀測偏北。然而,由于青海省強降水天氣個例嚴(yán)重偏少,還不足以概括每一種模式產(chǎn)品在本省較穩(wěn)定的偏差。另外,MODE檢驗結(jié)果對參數(shù)的選擇較為敏感,如平滑半徑、濾波閾值、合并距離等,致使檢驗結(jié)果缺乏一致性,可能會引起一定誤導(dǎo)。然而,鄰域方法不需要選擇參數(shù),能夠體現(xiàn)不同鄰域空間尺度下評分變化,對檢驗一些“模糊的”不包含匹配良好對象的預(yù)報表現(xiàn)較好。隨著高分辨率數(shù)值模式的廣泛應(yīng)用,如何有效利用空間檢驗方法以提高模式釋用還有待更進一步研究。

    久久中文字幕一级| 岛国毛片在线播放| 亚洲精品国产精品久久久不卡| 一夜夜www| 中亚洲国语对白在线视频| 久久性视频一级片| 国产精品久久久久久人妻精品电影| 久久国产精品大桥未久av| 欧美日韩亚洲综合一区二区三区_| 精品免费久久久久久久清纯 | 中文字幕av电影在线播放| 50天的宝宝边吃奶边哭怎么回事| 久久国产精品人妻蜜桃| 午夜免费鲁丝| 窝窝影院91人妻| 91精品国产国语对白视频| 亚洲熟女毛片儿| 国产精品一区二区精品视频观看| 在线观看免费视频日本深夜| 国产91精品成人一区二区三区| 欧美在线一区亚洲| 色综合欧美亚洲国产小说| 91在线观看av| 搡老岳熟女国产| 叶爱在线成人免费视频播放| 女人被狂操c到高潮| 91精品国产国语对白视频| 啦啦啦 在线观看视频| 黄片大片在线免费观看| xxxhd国产人妻xxx| 男人舔女人的私密视频| 国产精品自产拍在线观看55亚洲 | 亚洲av成人不卡在线观看播放网| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产乱子伦一区二区三区| 下体分泌物呈黄色| 久久精品国产亚洲av香蕉五月 | 热re99久久精品国产66热6| 丰满的人妻完整版| 成人黄色视频免费在线看| xxx96com| ponron亚洲| 亚洲精品美女久久久久99蜜臀| 又大又爽又粗| 高潮久久久久久久久久久不卡| 精品国产亚洲在线| 99热只有精品国产| 男女免费视频国产| 成人黄色视频免费在线看| 久久草成人影院| 国产精品一区二区在线不卡| 国精品久久久久久国模美| 久久精品亚洲精品国产色婷小说| 水蜜桃什么品种好| 亚洲精品国产区一区二| 人人妻,人人澡人人爽秒播| 日韩精品免费视频一区二区三区| 一进一出好大好爽视频| 天堂动漫精品| av有码第一页| 精品免费久久久久久久清纯 | 亚洲三区欧美一区| 大型黄色视频在线免费观看| 一区福利在线观看| 女性生殖器流出的白浆| 亚洲精品自拍成人| 两人在一起打扑克的视频| 最近最新中文字幕大全电影3 | 久久精品91无色码中文字幕| 久久久精品国产亚洲av高清涩受| 免费av中文字幕在线| 国产99久久九九免费精品| 久久精品亚洲精品国产色婷小说| 一边摸一边做爽爽视频免费| 国产深夜福利视频在线观看| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av香蕉五月 | 精品高清国产在线一区| 久久久久国产精品人妻aⅴ院 | 欧美乱妇无乱码| 亚洲成a人片在线一区二区| 人妻久久中文字幕网| 亚洲 欧美一区二区三区| www日本在线高清视频| 美女视频免费永久观看网站| 天天影视国产精品| 亚洲av日韩精品久久久久久密| 亚洲精华国产精华精| 美女 人体艺术 gogo| 校园春色视频在线观看| 91老司机精品| 天天躁日日躁夜夜躁夜夜| 少妇 在线观看| 欧美+亚洲+日韩+国产| 80岁老熟妇乱子伦牲交| 一个人免费在线观看的高清视频| 母亲3免费完整高清在线观看| 无人区码免费观看不卡| 黄片播放在线免费| 欧美激情久久久久久爽电影 | 久9热在线精品视频| 成人国语在线视频| 黄色a级毛片大全视频| 老熟妇乱子伦视频在线观看| 久久婷婷成人综合色麻豆| 精品亚洲成国产av| 亚洲视频免费观看视频| 90打野战视频偷拍视频| 成年人黄色毛片网站| 最近最新中文字幕大全电影3 | 1024视频免费在线观看| 最近最新中文字幕大全电影3 | 国产成人影院久久av| 两个人看的免费小视频| 乱人伦中国视频| 91av网站免费观看| 国产日韩一区二区三区精品不卡| 国产深夜福利视频在线观看| 满18在线观看网站| 国产日韩一区二区三区精品不卡| 精品福利观看| 色播在线永久视频| 精品一区二区三区av网在线观看| 国精品久久久久久国模美| 好看av亚洲va欧美ⅴa在| 色老头精品视频在线观看| 精品视频人人做人人爽| 国产欧美日韩一区二区精品| 亚洲自偷自拍图片 自拍| 国产淫语在线视频| 日韩大码丰满熟妇| 欧美最黄视频在线播放免费 | 亚洲精品久久成人aⅴ小说| 一边摸一边抽搐一进一小说 | 黄色视频,在线免费观看| 可以免费在线观看a视频的电影网站| 欧美日韩av久久| 久久精品国产清高在天天线| 国产精品秋霞免费鲁丝片| 日韩免费高清中文字幕av| av在线播放免费不卡| 十八禁人妻一区二区| 精品国内亚洲2022精品成人 | 成人黄色视频免费在线看| 色综合婷婷激情| 国产精品久久久av美女十八| 成年人黄色毛片网站| 亚洲中文日韩欧美视频| 一区二区日韩欧美中文字幕| 后天国语完整版免费观看| 国产在线精品亚洲第一网站| 9热在线视频观看99| 三上悠亚av全集在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 天堂√8在线中文| 亚洲七黄色美女视频| 热re99久久国产66热| 色播在线永久视频| 亚洲成人手机| 夜夜夜夜夜久久久久| 国产精品 欧美亚洲| 国产在视频线精品| 亚洲专区中文字幕在线| 免费黄频网站在线观看国产| 久久中文看片网| 色94色欧美一区二区| 丝袜人妻中文字幕| 欧美一级毛片孕妇| 9191精品国产免费久久| 日韩免费av在线播放| 99在线人妻在线中文字幕 | 999精品在线视频| 国产男女超爽视频在线观看| 9色porny在线观看| av天堂在线播放| 日韩人妻精品一区2区三区| av国产精品久久久久影院| 天天躁狠狠躁夜夜躁狠狠躁| 精品熟女少妇八av免费久了| 亚洲精品一卡2卡三卡4卡5卡| 欧美日本中文国产一区发布| 免费久久久久久久精品成人欧美视频| av线在线观看网站| 欧美日韩成人在线一区二区| 老司机亚洲免费影院| 亚洲精品一二三| 99久久国产精品久久久| 麻豆乱淫一区二区| 一夜夜www| 精品一区二区三区av网在线观看| av欧美777| 在线观看日韩欧美| 亚洲专区国产一区二区| 国产有黄有色有爽视频| av国产精品久久久久影院| 在线观看一区二区三区激情| 免费在线观看黄色视频的| 在线十欧美十亚洲十日本专区| 亚洲,欧美精品.| 精品一区二区三区四区五区乱码| 精品卡一卡二卡四卡免费| 好男人电影高清在线观看| av一本久久久久| 国产精品永久免费网站| 18禁裸乳无遮挡动漫免费视频| 久久人人爽av亚洲精品天堂| 成年版毛片免费区| 巨乳人妻的诱惑在线观看| 色94色欧美一区二区| 波多野结衣一区麻豆| 两个人看的免费小视频| 妹子高潮喷水视频| 18禁国产床啪视频网站| 国产成人精品在线电影| 精品电影一区二区在线| 久久精品aⅴ一区二区三区四区| 国产伦人伦偷精品视频| 一区二区三区精品91| 亚洲国产中文字幕在线视频| 9色porny在线观看| 亚洲三区欧美一区| 亚洲,欧美精品.| 18在线观看网站| 国产精品一区二区在线不卡| 亚洲五月色婷婷综合| www.999成人在线观看| 成人国语在线视频| 纯流量卡能插随身wifi吗| 国产免费现黄频在线看| 中出人妻视频一区二区| 悠悠久久av| 精品福利永久在线观看| 老司机深夜福利视频在线观看| 亚洲男人天堂网一区| 国产单亲对白刺激| 丝瓜视频免费看黄片| 久久精品成人免费网站| 美女国产高潮福利片在线看| 国内毛片毛片毛片毛片毛片| 怎么达到女性高潮| 1024香蕉在线观看| 色综合欧美亚洲国产小说| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩一区二区精品| 亚洲成人免费电影在线观看| 久久久久久久久久久久大奶| 老司机亚洲免费影院| 美女高潮到喷水免费观看| 国内久久婷婷六月综合欲色啪| 十八禁高潮呻吟视频| 中文欧美无线码| 美女国产高潮福利片在线看| 18在线观看网站| 色播在线永久视频| 国产一区二区三区综合在线观看| 国产日韩一区二区三区精品不卡| 国产精品久久久人人做人人爽| 国产精品美女特级片免费视频播放器 | 国产精品一区二区免费欧美| 麻豆av在线久日| 精品国产一区二区三区久久久樱花| 精品一区二区三区四区五区乱码| 亚洲午夜理论影院| 男女午夜视频在线观看| 老司机影院毛片| 国产麻豆69| 好看av亚洲va欧美ⅴa在| 女警被强在线播放| 男女高潮啪啪啪动态图| 亚洲一区中文字幕在线| 一a级毛片在线观看| a级毛片在线看网站| 国产精品一区二区在线不卡| av不卡在线播放| 日韩制服丝袜自拍偷拍| av有码第一页| 日本wwww免费看| 热re99久久国产66热| 超碰97精品在线观看| 亚洲av电影在线进入| 欧美成人免费av一区二区三区 | 国产精品自产拍在线观看55亚洲 | 久久人妻熟女aⅴ| 国产1区2区3区精品| 久久中文看片网| 欧美成人免费av一区二区三区 | 午夜精品久久久久久毛片777| 变态另类成人亚洲欧美熟女 | 国产男女超爽视频在线观看| 精品高清国产在线一区| 欧美日韩中文字幕国产精品一区二区三区 | 一二三四社区在线视频社区8| 国内久久婷婷六月综合欲色啪| 欧美av亚洲av综合av国产av| 十分钟在线观看高清视频www| 亚洲专区中文字幕在线| 欧美日韩福利视频一区二区| 老司机靠b影院| 丝袜人妻中文字幕| 欧美日韩av久久| 免费女性裸体啪啪无遮挡网站| 美女福利国产在线| 久久香蕉国产精品| 日日摸夜夜添夜夜添小说| 99热只有精品国产| 国内久久婷婷六月综合欲色啪| 搡老岳熟女国产| 黄色片一级片一级黄色片| 久久精品亚洲熟妇少妇任你| 黄片播放在线免费| 曰老女人黄片| 亚洲男人天堂网一区| 成人国语在线视频| 精品国产亚洲在线| 亚洲av日韩精品久久久久久密| 欧美成人免费av一区二区三区 | 亚洲片人在线观看| 制服诱惑二区| 日日摸夜夜添夜夜添小说| 亚洲国产精品sss在线观看 | 国产欧美亚洲国产| 国产精品秋霞免费鲁丝片| 香蕉久久夜色| 欧美午夜高清在线| 在线观看日韩欧美| av不卡在线播放| 日韩欧美三级三区| 美女 人体艺术 gogo| 成年版毛片免费区| а√天堂www在线а√下载 | 久久精品亚洲av国产电影网| 成人特级黄色片久久久久久久| 一边摸一边抽搐一进一小说 | 国产精品美女特级片免费视频播放器 | 亚洲在线自拍视频| 久久中文看片网| 国产片内射在线| 亚洲片人在线观看| 亚洲精品国产一区二区精华液| 国产又爽黄色视频| 在线av久久热| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 婷婷成人精品国产| 亚洲欧美色中文字幕在线| 亚洲黑人精品在线| xxx96com| 色在线成人网| 18在线观看网站| 一级毛片精品| 国产xxxxx性猛交| 99热国产这里只有精品6| 无限看片的www在线观看| 国产不卡av网站在线观看| 18禁黄网站禁片午夜丰满| 99国产精品99久久久久| 可以免费在线观看a视频的电影网站| 9热在线视频观看99| 亚洲精品久久成人aⅴ小说| 午夜精品久久久久久毛片777| 亚洲专区字幕在线| 久久精品国产亚洲av香蕉五月 | 777久久人妻少妇嫩草av网站| 亚洲成av片中文字幕在线观看| 国产成人免费观看mmmm| 国产精品亚洲一级av第二区| 99精品欧美一区二区三区四区| 日韩有码中文字幕| 又黄又爽又免费观看的视频| 日本黄色日本黄色录像| 校园春色视频在线观看| 美女视频免费永久观看网站| 国产成人精品无人区| 亚洲 欧美一区二区三区| av欧美777| 黄色片一级片一级黄色片| 女同久久另类99精品国产91| 久久精品国产亚洲av高清一级| 成人精品一区二区免费| 天天影视国产精品| 少妇粗大呻吟视频| av天堂在线播放| av超薄肉色丝袜交足视频| 在线观看免费高清a一片| 热99久久久久精品小说推荐| 老司机午夜福利在线观看视频| 日韩免费av在线播放| 母亲3免费完整高清在线观看| 国产亚洲精品一区二区www | 国产熟女午夜一区二区三区| 757午夜福利合集在线观看| 91精品三级在线观看| 多毛熟女@视频| 亚洲精品在线观看二区| 桃红色精品国产亚洲av| 超碰成人久久| 亚洲av电影在线进入| 日日爽夜夜爽网站| av视频免费观看在线观看| 在线av久久热| 国产不卡一卡二| www.熟女人妻精品国产| 99国产精品免费福利视频| 久久青草综合色| 国产一区二区激情短视频| 在线国产一区二区在线| 51午夜福利影视在线观看| a级毛片黄视频| 国产精品 国内视频| 久久香蕉精品热| 久久性视频一级片| 天天躁夜夜躁狠狠躁躁| 狂野欧美激情性xxxx| 午夜福利乱码中文字幕| 淫妇啪啪啪对白视频| 99在线人妻在线中文字幕 | av网站免费在线观看视频| 亚洲欧美精品综合一区二区三区| 欧美精品一区二区免费开放| 中文字幕制服av| 极品人妻少妇av视频| 亚洲国产毛片av蜜桃av| 夫妻午夜视频| 久久久久精品人妻al黑| 搡老熟女国产l中国老女人| 18禁美女被吸乳视频| 欧美精品亚洲一区二区| 后天国语完整版免费观看| 好男人电影高清在线观看| 久久久久久亚洲精品国产蜜桃av| 一区福利在线观看| 国产在线精品亚洲第一网站| 人人妻人人澡人人爽人人夜夜| 午夜福利影视在线免费观看| av超薄肉色丝袜交足视频| 日本黄色视频三级网站网址 | 久久国产精品男人的天堂亚洲| svipshipincom国产片| 欧美黄色片欧美黄色片| 在线国产一区二区在线| 人成视频在线观看免费观看| 妹子高潮喷水视频| 日本vs欧美在线观看视频| 国产真人三级小视频在线观看| 男女高潮啪啪啪动态图| 欧美激情极品国产一区二区三区| 又黄又粗又硬又大视频| 国产精品av久久久久免费| 精品高清国产在线一区| 免费少妇av软件| 天堂√8在线中文| 亚洲精品国产一区二区精华液| 天天躁狠狠躁夜夜躁狠狠躁| 久久久精品国产亚洲av高清涩受| 成人av一区二区三区在线看| 久久精品亚洲av国产电影网| 人人妻人人澡人人爽人人夜夜| 99国产精品一区二区蜜桃av | 国产不卡av网站在线观看| 国产高清国产精品国产三级| av有码第一页| 色老头精品视频在线观看| 在线观看午夜福利视频| 大型黄色视频在线免费观看| 一级片'在线观看视频| 国产一卡二卡三卡精品| 黄色毛片三级朝国网站| 国产成人精品久久二区二区91| 女同久久另类99精品国产91| 99久久国产精品久久久| 国产精品香港三级国产av潘金莲| 高潮久久久久久久久久久不卡| 大陆偷拍与自拍| 欧美国产精品va在线观看不卡| 最新的欧美精品一区二区| 午夜精品国产一区二区电影| 成人三级做爰电影| 91麻豆精品激情在线观看国产 | 日韩欧美国产一区二区入口| 国产精品久久电影中文字幕 | 国产人伦9x9x在线观看| 亚洲av日韩精品久久久久久密| 欧美国产精品va在线观看不卡| 桃红色精品国产亚洲av| 女人高潮潮喷娇喘18禁视频| 国产高清videossex| 国产午夜精品久久久久久| 国产精品久久久久成人av| 18在线观看网站| 欧美+亚洲+日韩+国产| 国产亚洲欧美精品永久| 国产淫语在线视频| www.精华液| 999久久久精品免费观看国产| 亚洲第一av免费看| 国产熟女午夜一区二区三区| 一级毛片高清免费大全| 高清视频免费观看一区二区| 中文字幕色久视频| 日韩欧美在线二视频 | 19禁男女啪啪无遮挡网站| 国产精品永久免费网站| 美女高潮喷水抽搐中文字幕| 日日夜夜操网爽| 亚洲av美国av| 真人做人爱边吃奶动态| 亚洲精品美女久久av网站| 一级毛片女人18水好多| 视频区欧美日本亚洲| 搡老熟女国产l中国老女人| 免费高清在线观看日韩| 亚洲精品av麻豆狂野| 精品国产美女av久久久久小说| xxx96com| 老汉色av国产亚洲站长工具| 丝袜美足系列| 亚洲人成77777在线视频| 免费观看精品视频网站| 日韩三级视频一区二区三区| 精品国产乱子伦一区二区三区| 亚洲精品一二三| 91大片在线观看| 日韩欧美三级三区| 亚洲在线自拍视频| 亚洲av成人av| 亚洲精品av麻豆狂野| 国产精品.久久久| 王馨瑶露胸无遮挡在线观看| 国产免费现黄频在线看| 夜夜夜夜夜久久久久| 99久久99久久久精品蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 中文字幕av电影在线播放| 成年人免费黄色播放视频| 男男h啪啪无遮挡| 午夜福利在线免费观看网站| 欧美日韩国产mv在线观看视频| 欧美日韩乱码在线| 不卡av一区二区三区| 51午夜福利影视在线观看| 欧美日韩瑟瑟在线播放| 精品福利永久在线观看| 亚洲熟女精品中文字幕| 亚洲欧美日韩另类电影网站| 青草久久国产| 午夜视频精品福利| 久久久久久亚洲精品国产蜜桃av| 深夜精品福利| 国产精品乱码一区二三区的特点 | 一个人免费在线观看的高清视频| 一级黄色大片毛片| 国产区一区二久久| 亚洲男人天堂网一区| 超碰成人久久| 中文字幕人妻熟女乱码| 午夜福利乱码中文字幕| 成人18禁高潮啪啪吃奶动态图| 久久久久精品人妻al黑| 丝袜美足系列| 午夜视频精品福利| 俄罗斯特黄特色一大片| 中文字幕精品免费在线观看视频| а√天堂www在线а√下载 | 深夜精品福利| 一级毛片精品| 精品久久久久久久毛片微露脸| 精品久久久精品久久久| 男女免费视频国产| 变态另类成人亚洲欧美熟女 | 一a级毛片在线观看| 久久影院123| 亚洲精品国产精品久久久不卡| 黑人操中国人逼视频| 国产成人精品无人区| 精品一区二区三区视频在线观看免费 | tube8黄色片| 欧美中文综合在线视频| 欧美午夜高清在线| 亚洲精品久久午夜乱码| 一个人免费在线观看的高清视频| 亚洲自偷自拍图片 自拍| 天天操日日干夜夜撸| 亚洲熟女毛片儿| 国产极品粉嫩免费观看在线| 日韩制服丝袜自拍偷拍| 女人被狂操c到高潮| 日韩欧美一区视频在线观看| 久久人妻av系列| 女人久久www免费人成看片| 香蕉久久夜色| 最近最新免费中文字幕在线| 国产无遮挡羞羞视频在线观看| 美女 人体艺术 gogo| 美女午夜性视频免费| av天堂久久9| 久久中文字幕一级| 国产精品久久久久久精品古装| 亚洲avbb在线观看| 午夜福利免费观看在线| 丰满饥渴人妻一区二区三| 大码成人一级视频| 亚洲中文av在线| 成人特级黄色片久久久久久久| av网站在线播放免费| 国产男女超爽视频在线观看| 婷婷成人精品国产| www.熟女人妻精品国产| 久久人妻熟女aⅴ| 久久人人97超碰香蕉20202| 久久ye,这里只有精品| 男女免费视频国产| av线在线观看网站| 人人妻人人添人人爽欧美一区卜| 欧美精品人与动牲交sv欧美| 亚洲成人免费av在线播放|