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

    基于雷暴尺度集合預(yù)報(bào)的頻率匹配降水預(yù)報(bào)研究

    2022-02-06 03:20:38陳訓(xùn)來朱江山林良勛王德立陳元昭孔凡鈾
    氣象與環(huán)境科學(xué) 2022年6期
    關(guān)鍵詞:個(gè)例雷暴強(qiáng)降水

    陳訓(xùn)來, 朱江山, 林良勛, 王德立, 王 蕊, 陳元昭, 孔凡鈾

    (1.深圳市氣象局,深圳 518040; 2.深圳南方強(qiáng)天氣研究重點(diǎn)實(shí)驗(yàn)室,深圳 518040; 3.中國(guó)科學(xué)院大氣物理研究所,北京 100029; 4.廣東省氣象臺(tái),廣州 510080; 5.美國(guó)俄克拉荷馬大學(xué)風(fēng)暴分析和預(yù)測(cè)中心,諾曼 732072)

    引 言

    大氣運(yùn)動(dòng)本身的混沌特性、模式誤差、參數(shù)化方案誤差等因素,導(dǎo)致數(shù)值預(yù)報(bào)存在不可避免的不確定性[1]。為了估計(jì)預(yù)報(bào)的不確定性,集合預(yù)報(bào)方法被廣泛使用在數(shù)值預(yù)報(bào)的研究和業(yè)務(wù)應(yīng)用當(dāng)中。集合預(yù)報(bào)不僅能夠提供具有量化基礎(chǔ)的概率預(yù)報(bào),反映預(yù)報(bào)的不確定性,還能夠提供一個(gè)技巧較高的確定性預(yù)報(bào),如集合平均預(yù)報(bào)[2]。目前,已經(jīng)發(fā)展出了許多基于集合預(yù)報(bào)的概率預(yù)報(bào)方法,如邏輯回歸[3]、貝葉斯模式平均[4-6]、集合敷料法[7-8]、非均勻高斯回歸[9]等。代刊等[10]對(duì)各種集合預(yù)報(bào)后處理方法作了較全面的綜述。

    與此同時(shí),基于集合預(yù)報(bào)的確定性預(yù)報(bào)仍然被研究者和預(yù)報(bào)員廣泛研究和使用,例如集合平均預(yù)報(bào)[2,11]。集合平均預(yù)報(bào)通過集合成員之間的算術(shù)平均,濾除各集合成員的差異部分,從而減小預(yù)報(bào)的不確定性,提高預(yù)報(bào)準(zhǔn)確率,常被作為基于集合預(yù)報(bào)的最優(yōu)確定性預(yù)報(bào)。但是,對(duì)于時(shí)空分布變率較大的物理量(如降水)來說,集合平均預(yù)報(bào)方法存在較為明顯的不足。雖然,通過過濾各集合成員之間的不同,集合平均預(yù)報(bào)對(duì)降水位置的預(yù)報(bào)往往優(yōu)于大多數(shù)集合成員,并且這一優(yōu)勢(shì)在長(zhǎng)期統(tǒng)計(jì)中會(huì)更加明顯[12],但是,當(dāng)各集合成員的降水場(chǎng)的時(shí)空分布差異較大時(shí),集合成員之間的算術(shù)平均常常會(huì)過度平滑降水場(chǎng)的強(qiáng)度,造成非真實(shí)的降水強(qiáng)度分布,常體現(xiàn)為擴(kuò)大的降水范圍和減弱的降水強(qiáng)度。針對(duì)這一問題,Ebert[12]使用頻率匹配(Frequency Matching,FM)方法,按照等頻率映射關(guān)系,將集合平均預(yù)報(bào)的累計(jì)頻率分布函數(shù)(Cumulative Distribution Function,CDF)映射到所有集合成員的CDF上去,從而改善了集合平均預(yù)報(bào)的強(qiáng)度分布偏差。Zhu等[13]又進(jìn)一步使用頻率匹配方法,將全球集合預(yù)報(bào)各成員預(yù)報(bào)的CDF訂正到降水觀測(cè)資料的CDF上,較好地訂正了各成員預(yù)報(bào)的強(qiáng)度偏差,同時(shí)也改善集合平均預(yù)報(bào)的效果。李俊等[14]對(duì)華中地區(qū)AREM模式確定性預(yù)報(bào)進(jìn)行了頻率匹配訂正試驗(yàn),也得到相似的結(jié)果。李俊等[15]對(duì)集合預(yù)報(bào)進(jìn)行頻率匹配訂正試驗(yàn),指出直接對(duì)集合平均進(jìn)行頻率匹配訂正的效果優(yōu)于對(duì)各集合成員分別作頻率匹配訂正再求算術(shù)平均。李莉等[16]使用頻率匹配方法訂正了T213降水預(yù)報(bào),結(jié)果表明該方法對(duì)T213降水預(yù)報(bào)的偏差有明顯改善。周迪等[17]針對(duì)四川盆地暴雨,采用Gamma函數(shù)來擬合頻率分布曲線,進(jìn)而使用頻率匹配方法有效地訂正了T213暴雨集合預(yù)報(bào)系統(tǒng)性誤差。以上研究發(fā)現(xiàn),頻率匹配訂正法對(duì)大范圍降水過程訂正效果較好,尤其是對(duì)有無降水預(yù)報(bào)和50 mm以下量級(jí)降水預(yù)報(bào)的改善較為明顯。但由于頻率匹配降水預(yù)報(bào)訂正法在處理預(yù)報(bào)和觀測(cè)降水頻率時(shí),均是針對(duì)某一關(guān)鍵區(qū)域進(jìn)行統(tǒng)計(jì)分析和訂正,不同區(qū)域的累積概率分布函數(shù)具有明顯的差異,因此訂正效果受不同區(qū)域天氣氣候特征影響較大[16-17]。目前針對(duì)華南地區(qū)頻率匹配訂正方法的研究和業(yè)務(wù)應(yīng)用還比較少,有必要進(jìn)一步開展針對(duì)性的研究。

    由于對(duì)流尺度天氣系統(tǒng)時(shí)空尺度小、動(dòng)力與物理過程高度非線性,對(duì)流天氣的精細(xì)化預(yù)報(bào)面臨巨大的挑戰(zhàn),因此,發(fā)展高分辨率的雷暴尺度集合預(yù)報(bào)顯得尤為重要。近十幾年來,高性能計(jì)算機(jī)的發(fā)展,為雷暴尺度集合預(yù)報(bào)巨大計(jì)算量的問題提供了必要、可行的計(jì)算環(huán)境,相關(guān)研究和應(yīng)用日益受到國(guó)內(nèi)外學(xué)者和數(shù)值預(yù)報(bào)業(yè)務(wù)機(jī)構(gòu)重視,在雷暴尺度集合預(yù)報(bào)模式初值擾動(dòng)技術(shù)、集合同化技術(shù)等方面取得了豐碩的研究成果,為數(shù)不少的業(yè)務(wù)系統(tǒng)也投入運(yùn)行[18-19],但頻率匹配訂正方法在雷暴尺度集合降水預(yù)報(bào)中的應(yīng)用研究還不多。深圳市氣象局自2013年開始引進(jìn)和應(yīng)用美國(guó)俄克拉荷馬大學(xué)雷暴分析和預(yù)報(bào)中心(CAPS)技術(shù),發(fā)展了深圳雷暴尺度(4 km分辨率)集合預(yù)報(bào)系統(tǒng)(Shenzhen Storm-Scale Ensemble Forecast system,SZ-SSEF)),在華南暴雨、強(qiáng)對(duì)流等預(yù)報(bào)業(yè)務(wù)中發(fā)揮了重要的參考價(jià)值[20-21]。江崟等[21]針對(duì)深圳雷暴尺度集合預(yù)報(bào)系統(tǒng),分析了頻率匹配方法的降水偏差訂正技術(shù)在典型暴雨過程的應(yīng)用,但缺少更多降水個(gè)例的分析。

    華南地處低緯,瀕臨海洋,受到低緯度熱帶天氣系統(tǒng)與中高緯度天氣系統(tǒng)的交替影響,天氣氣候異常復(fù)雜多變,暴雨發(fā)生頻率高、雨強(qiáng)大、致災(zāi)性強(qiáng),皆居全國(guó)前列[22-23]。因此,本文選取了華南地區(qū)2018年汛期4-9月103個(gè)降水個(gè)例,對(duì)深圳市氣象局雷暴尺度集合預(yù)報(bào)系統(tǒng)進(jìn)行了批量試驗(yàn),以檢驗(yàn)頻率匹配方法在雷暴尺度集合降水預(yù)報(bào)中的優(yōu)點(diǎn)和不足。華南汛期持續(xù)時(shí)間長(zhǎng),針對(duì)華南地區(qū)汛期降水預(yù)報(bào)進(jìn)行頻率匹配訂正研究,可以更全面地評(píng)估頻率匹配方法在降水預(yù)報(bào)中的效果,同時(shí)對(duì)于防災(zāi)減災(zāi)救災(zāi)也具有較好的實(shí)際應(yīng)用意義。

    1 資料與方法

    1.1 SZ-SSEF系統(tǒng)

    SZ-SSEF采用單向兩重嵌套網(wǎng)格,其中外層區(qū)域12 km網(wǎng)格距中尺度集合預(yù)報(bào)(D01)的水平格點(diǎn)數(shù)為211×179,內(nèi)層區(qū)域4 km網(wǎng)格距雷暴尺度集合預(yù)報(bào)(D02)的水平格點(diǎn)數(shù)為423×363,垂直方向?yàn)?1層(圖1、表1)。外層區(qū)域采用歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts,ECMWF)的0.125 °× 0.125 °資料和美國(guó)環(huán)境預(yù)報(bào)中心(National Centers for Environmental Prediction,NCEP)全球集合預(yù)報(bào)系統(tǒng)(Global Ensemble Forecast System,GEFS)1°×1°資料作為初始場(chǎng)和邊界條件,內(nèi)層成員利用美國(guó)風(fēng)暴分析和預(yù)報(bào)中心開發(fā)的ARPS(The Advanced Regional Prediction System)模式的三維變分系統(tǒng)(Three-dimensional variational,3DVAR)和云分析模塊分別同化廣東省全省11部多普勒雷達(dá)的徑向風(fēng)和反射率資料,預(yù)報(bào)部分采用WRF (Weather Research and Forecast Model)模式V3.5.1[18]。

    表1 深圳雷暴尺度集合預(yù)報(bào)系統(tǒng) 4 km區(qū)域模式成員的配置

    圖1 深圳雷暴尺度集合預(yù)報(bào)系統(tǒng)(SZ-SSEF)模式區(qū)域設(shè)置

    Du等[24]研究認(rèn)為,8個(gè)集合成員就可以獲得集合平均90%的效果。高峰等[25]研究認(rèn)為,集合預(yù)報(bào)效果總體上預(yù)報(bào)技巧隨成員數(shù)增加而提高,但當(dāng)集合成員數(shù)達(dá)到5~13 時(shí),預(yù)報(bào)技巧呈飽和特征。考慮到計(jì)算能力及性價(jià)比,目前深圳雷暴尺度集合預(yù)報(bào)系統(tǒng)成員數(shù)為10個(gè),其中M00成員采用ECMWF作為初始場(chǎng)和邊界條件,其他集合成員采用GEFS作為初始場(chǎng)和邊界條件,不同成員選擇不同的微物理過程和邊界層參數(shù)化方案,最終形成多模式初值、多微物理過程和多邊界層過程的雷暴尺度集合預(yù)報(bào)系統(tǒng)。深圳雷暴尺度集合預(yù)報(bào)系統(tǒng)內(nèi)層網(wǎng)格區(qū)域每間隔3 h啟動(dòng)一次,每天啟動(dòng)8次,預(yù)報(bào)時(shí)效為36 h。目前,該系統(tǒng)已經(jīng)成為深圳市氣象局業(yè)務(wù)預(yù)報(bào)的一個(gè)主要支撐[20-21]。

    1.2 預(yù)報(bào)和檢驗(yàn)資料

    本文收集了SZ-SSEF系統(tǒng)在2017年和2018年汛期4-9月的266次華南地區(qū)降水歷史預(yù)報(bào)個(gè)例,其中,2017年163次個(gè)例,2018年103次個(gè)例。選取了SZ-SSEF內(nèi)層4 km網(wǎng)格,每日12 UTC起報(bào)(世界時(shí),下同),預(yù)報(bào)時(shí)間12-36 h時(shí)段內(nèi)的24 h累計(jì)降水的集合平均預(yù)報(bào)作為訂正對(duì)象。因?yàn)閿?shù)據(jù)傳輸、計(jì)算機(jī)系統(tǒng)故障等因素 ,兩年的歷史預(yù)報(bào)資料都存在數(shù)據(jù)缺失的情況。

    本文使用中國(guó)國(guó)家氣象中心發(fā)布的0.05°×0.05°網(wǎng)格化逐時(shí)實(shí)況降水融合分析數(shù)據(jù)作為檢驗(yàn)數(shù)據(jù)[26-27]。降水格點(diǎn)分析數(shù)據(jù)首先被插值到SZ-SSEF的4 km內(nèi)層網(wǎng)格上,然后與預(yù)報(bào)進(jìn)行對(duì)比。

    1.3 頻率匹配(FM)訂正方法

    本文采用的頻率匹配FM訂正方法包括以下三個(gè)步驟進(jìn)行:

    (1)假設(shè)被訂正預(yù)報(bào)日期為d,選取d±45天的時(shí)間窗口內(nèi),當(dāng)年和往年所有可用的歷史預(yù)報(bào)個(gè)例,作為訓(xùn)練訂正模型的歷史個(gè)例。為了有足夠的樣本數(shù)而又大致考慮到近期相似的總體天氣形勢(shì),本文采取前后45天滑動(dòng)平均的方法計(jì)算降水的平均頻率。為模擬業(yè)務(wù)運(yùn)行,滿足訂正開始時(shí),觀測(cè)數(shù)據(jù)已經(jīng)可以獲得的條件,當(dāng)年d-1天之后的個(gè)例不用作訓(xùn)練樣本。

    (2)在整個(gè)SZ-SSEF 4 km區(qū)域上,分別計(jì)算被訂正預(yù)報(bào)的CDF和訓(xùn)練集個(gè)例的降水分析資料CDF。計(jì)算訓(xùn)練集降水分析數(shù)據(jù)的CDF時(shí),先將訓(xùn)練集中所有個(gè)例的降水分析數(shù)據(jù)集合到一起,再計(jì)算CDF。CDF的計(jì)算采用非參數(shù)方法,對(duì)0 mm到500 mm,每1 mm間隔的閾值,分別計(jì)算其對(duì)應(yīng)的累計(jì)頻率。

    (3)設(shè)被訂正預(yù)報(bào)和降水分析資料的CDF分別為CDFfcst和CDFobs。對(duì)任意格點(diǎn)上的原始預(yù)報(bào)fraw,計(jì)算其對(duì)應(yīng)的累計(jì)頻率CDFfcst(fraw),再將其帶入降水分析數(shù)據(jù)CDFobs的反函數(shù),求出降水分析數(shù)據(jù)CDFobs中相同頻率對(duì)應(yīng)的降水值,將其作為訂正預(yù)報(bào)fcal。其具體方法如下:

    (1)

    頻率匹配訂正方法的詳細(xì)解釋可參考文獻(xiàn)[13-14]。

    1.4 檢驗(yàn)方法

    本文使用Gilbert Skill Score (GSS)評(píng)分和BIAS評(píng)分對(duì)訂正前后的集合平均降水預(yù)報(bào)進(jìn)行檢驗(yàn)。兩項(xiàng)評(píng)分都是針對(duì)某一特定閾值,計(jì)算預(yù)報(bào)和觀測(cè)在某一方面的相似性,其計(jì)算都基于二分類預(yù)報(bào)聯(lián)列表,即將預(yù)報(bào)-觀測(cè)關(guān)系分為四類情況(表2)。

    表2 二分類預(yù)報(bào)聯(lián)列表

    GSS評(píng)分(又稱為Equitable Threat Score,ETS評(píng)分)針對(duì)某一特定閾值,檢驗(yàn)預(yù)報(bào)與觀測(cè)之間的重疊情況,并剔除隨機(jī)命中。其具體計(jì)算方法如下:

    (2)

    (3)

    其中,h為命中數(shù)量,hr為隨機(jī)命中數(shù)量,m為漏報(bào)數(shù)量,f為錯(cuò)警數(shù)量,t為檢驗(yàn)樣本總數(shù)。

    BIAS評(píng)分針對(duì)某一特定閾值,比較預(yù)報(bào)與觀測(cè)之間的面積差異,當(dāng)評(píng)分等于1時(shí),預(yù)報(bào)面積與觀測(cè)面積相等,為完美預(yù)報(bào);大于1時(shí),表示預(yù)報(bào)面積偏大(濕偏差);小于1時(shí),表示預(yù)報(bào)面積偏小(干偏差)。其具體計(jì)算方法如下:

    (4)

    1.5 試驗(yàn)設(shè)計(jì)

    本文選取2018年汛期4-9月華南地區(qū)汛期103次降水個(gè)例過程,對(duì)SZ-SSEF系統(tǒng)12時(shí)(UTC)起報(bào),預(yù)報(bào)時(shí)間12-36 h時(shí)段的24 h累計(jì)降水量的集合平均預(yù)報(bào)進(jìn)行了頻率匹配訂正預(yù)報(bào)試驗(yàn)。2017年汛期4-9月出現(xiàn)的163次降水個(gè)例數(shù)據(jù)只用作訂正模型的訓(xùn)練樣本。

    2 結(jié)果與分析

    2.1 個(gè)例試驗(yàn)1:2018年超強(qiáng)臺(tái)風(fēng)“山竹”

    超強(qiáng)臺(tái)風(fēng)“山竹”于2018年9月16日在江門市臺(tái)山海宴鎮(zhèn)沿海地區(qū)登陸,是2018年登陸我國(guó)最強(qiáng)的臺(tái)風(fēng),給華南大部分地區(qū)帶來了暴雨到大暴雨局部特大暴雨的降水,廣東省出現(xiàn)了歷史罕見的風(fēng)雨浪潮影響,全省平均雨量93.1 mm,有167個(gè)氣象水文站錄得超過250 mm的降水,其中陽(yáng)江陽(yáng)東大八鎮(zhèn)記錄到過程最大雨量481.5 mm。“山竹”造成廣東省18個(gè)市1060個(gè)鄉(xiāng)鎮(zhèn)受災(zāi),因?yàn)?zāi)死亡5人,直接經(jīng)濟(jì)損失超過75億元。

    這次臺(tái)風(fēng)過程是本文所有個(gè)例中降水最強(qiáng)的一次。本文首先選取該臺(tái)風(fēng)個(gè)例,對(duì)頻率匹配訂正的效果進(jìn)行檢驗(yàn)和分析。圖2是SZ-SSEF系統(tǒng)2018年9月15日12 UTC起報(bào),預(yù)報(bào)時(shí)間為12-36 h時(shí)段的24 h累計(jì)降水預(yù)報(bào)的原始集合平均預(yù)報(bào)、FM訂正預(yù)報(bào)和相應(yīng)的降水融合分析資料。圖2(a)中可以看到,原始集合平均預(yù)報(bào)的降水落區(qū)偏大,強(qiáng)降水區(qū)域(50 mm和100 mm量級(jí))的面積也明顯大于降水分析資料的面積,存在較為明顯的濕偏差。頻率匹配訂正預(yù)報(bào)(圖2b)相比于原始集合平均預(yù)報(bào),各個(gè)量級(jí)的降水面積都有較為明顯的收縮,與融合分析資料降水的分布更加接近(圖2c)。但是,原始預(yù)報(bào)和訂正預(yù)報(bào)的降水極值區(qū)域面積(超過200 mm)差別不大,都與融合分析資料的較為接近。

    圖2 2018年9月15日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段的24 h累計(jì)降水量預(yù)報(bào)單位:mm;(a)原始集合平均預(yù)報(bào),(b)FM訂正預(yù)報(bào),(c)降水分析資料

    圖3是這次降水過程整個(gè)模式區(qū)域上計(jì)算的原始集合平均預(yù)報(bào)、頻率匹配訂正預(yù)報(bào)和降水分析資料的累積頻率。

    圖3 2018年9月15日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段的24 h累計(jì)降水預(yù)報(bào)和相應(yīng)降水分析資料在模式區(qū)域D02上的CDF藍(lán)色為原始集合平均預(yù)報(bào),紅色為FM訂正預(yù)報(bào),黑色為降水分析資料

    預(yù)報(bào)累積頻率和融合分析資料累積頻率之間的比較,不考慮降水位置和形態(tài)的差異,單純檢驗(yàn)預(yù)報(bào)和融合分析資料之間的強(qiáng)度分布的相似性,可以看到,對(duì)于小于160 mm的降水,原始集合平均預(yù)報(bào)與降水分析數(shù)據(jù)存在較為明顯的差異,原始集合平均預(yù)報(bào)的累積頻率小于融合分析的累積頻率,表明模式系統(tǒng)集合平均預(yù)報(bào)降水具有系統(tǒng)性偏大的特征,而通過頻率匹配訂正預(yù)報(bào)后,降水預(yù)報(bào)的累積頻率與融合分析的較為接近,明顯優(yōu)于原始集合平均預(yù)報(bào)結(jié)果;對(duì)大于160 mm的強(qiáng)降水,原始集合平均預(yù)報(bào)累積頻率與實(shí)況分析較為一致,頻率匹配訂正預(yù)報(bào)后改進(jìn)效果不明顯。

    圖4是針對(duì)0.1 mm,1 mm,10 mm,25 mm,50 mm和100 mm等不同閾值的降水,原始集合平均預(yù)報(bào)和頻率匹配訂正預(yù)報(bào)的GSS評(píng)分和BIAS評(píng)分。由BIAS評(píng)分可以看出,原始集合平均預(yù)報(bào)對(duì)各個(gè)閾值的降水事件有過報(bào)現(xiàn)象。而頻率匹配訂正則較好地修正了原始集合平均預(yù)報(bào)的偏差,各個(gè)閾值的BIAS評(píng)分都與1較為接近。從GSS評(píng)分來看,訂正預(yù)報(bào)在0.1 mm至10 mm閾值上有優(yōu)勢(shì),說明訂正預(yù)報(bào)的小雨量級(jí)的降水落區(qū)更加準(zhǔn)確,這也與圖2中兩組預(yù)報(bào)的情況相符。在25 mm至100 mm閾值的GSS評(píng)分中,原始預(yù)報(bào)和訂正預(yù)報(bào)各有好壞,難以分出絕對(duì)的優(yōu)劣。這可能是因?yàn)闃O端降水預(yù)報(bào)中,降水中心的位置偏差對(duì)GSS評(píng)分的影響較大。原始集合平均預(yù)報(bào)中更大的降水面積更容易與觀測(cè)重合,對(duì)GSS評(píng)分有正貢獻(xiàn)。而頻率匹配訂正縮小了降水的區(qū)域,更容易發(fā)生錯(cuò)位,引起“錯(cuò)警”和“漏報(bào)”的雙重懲罰,對(duì)GSS評(píng)分有負(fù)影響[28-29]。

    圖4 2018年9月15日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段的24 h累計(jì)降水預(yù)報(bào)的GSS和BIAS評(píng)分藍(lán)色為原始集合平均預(yù)報(bào),紅色為FM訂正預(yù)報(bào)

    2.2 個(gè)例試驗(yàn)2:2018年5月9日西南季風(fēng)強(qiáng)降水

    2018年5月9-10日受西南季風(fēng)和切變線共同影響,在廣東西部和廣西東部出現(xiàn)了暴雨到大暴雨的降水,強(qiáng)降水主要集中在粵西市縣、江門、肇慶和廣西玉林,其中粵西陽(yáng)江最大降雨量達(dá)227.3 mm (圖5c)。SZ-SSEF系統(tǒng)較好地預(yù)報(bào)出了此次強(qiáng)降水過程,從2018年5月8日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段24 h累計(jì)降水預(yù)報(bào)的原始集合平均雨量分布來看(圖5a),大雨以上的降水區(qū)域主要在廣東西部沿海地區(qū),以及廣西東部地區(qū),暴雨以上降水主要集中在廣東陽(yáng)江、江門及廣西玉林地區(qū),預(yù)報(bào)降水強(qiáng)度和落區(qū)與實(shí)況分析較為一致,但SZ-SSEF系統(tǒng)的集合平均雨量偏大,強(qiáng)降水落區(qū)范圍較實(shí)況區(qū)域偏大,降水空?qǐng)?bào)較為明顯。通過頻率匹配方法訂正后,降水量級(jí)預(yù)報(bào)偏大的現(xiàn)象得到訂正,強(qiáng)降水區(qū)域范圍也明顯減小(圖5b),與實(shí)況分析區(qū)域更加接近,降低了降水空?qǐng)?bào)的現(xiàn)象。另外,由于SZ-SSEF系統(tǒng)原始集合平均對(duì)福建和江西交界處的強(qiáng)降水中心預(yù)報(bào)位置偏北、降水范圍偏小,頻率匹配方法訂正后對(duì)降水落區(qū)和強(qiáng)度改進(jìn)效果不明顯,這是由于頻率匹配方法的頻率統(tǒng)計(jì)是在一個(gè)區(qū)域內(nèi)進(jìn)行的,而沒有考慮具體降雨位置、雨區(qū)的形狀或走向誤差,所以頻率匹配方法無法訂正雨區(qū)位置的偏差[14]。

    圖5 2018年5月8日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段的24 h累計(jì)降水量預(yù)報(bào)單位:mm;(a)原始集合平均預(yù)報(bào),(b)FM訂正預(yù)報(bào),(c)降水分析資料

    圖6為此次季風(fēng)強(qiáng)降水過程累積頻率曲線分布圖。從圖中可以看出,對(duì)小于60 mm的降水,頻率匹配訂正后,降水預(yù)報(bào)的累積頻率明顯優(yōu)于原始集合平均預(yù)報(bào)的,尤其對(duì)小于10 mm的降水訂正效果比較明顯,降水預(yù)報(bào)的累積頻率與實(shí)況分析基本吻合;而對(duì)大于60 mm的降水,由于原始集合平均預(yù)報(bào)累積頻率與實(shí)況分析較為一致,頻率匹配訂正預(yù)報(bào)的改進(jìn)有限。

    圖6 2018年5月8日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段的24 h累計(jì)降水預(yù)報(bào)和相應(yīng)降水分析資料在模式區(qū)域D02上的CDF藍(lán)色為原始集合平均預(yù)報(bào),紅色為FM訂正預(yù)報(bào),黑色為降水分析資料

    從頻率匹配方法訂正前后預(yù)報(bào)的GSS和BIAS評(píng)分來看(圖7),經(jīng)過頻率匹配訂正后,對(duì)不同閾值的降水預(yù)報(bào)GSS評(píng)分和BIAS偏差評(píng)分都得到提高,GSS評(píng)分都高于訂正前的結(jié)果,其中對(duì)晴雨(0.1 mm)和大暴雨以上(100 mm)的預(yù)報(bào)評(píng)分提高較為明顯。對(duì)0.1 mm至50 mm量級(jí)的降水預(yù)測(cè),原始集合平均預(yù)報(bào)BIAS評(píng)分大于1,且隨著量級(jí)的增大BIAS評(píng)分為增大的趨勢(shì),說明這次過程有明顯的濕偏差,但對(duì)100 mm量級(jí)以上的降水,BIAS評(píng)分小于1,存在干偏差的情況,有漏報(bào)的現(xiàn)象,這從圖5中也可以看出。頻率匹配方法訂正后,對(duì)各個(gè)閾值的降水預(yù)報(bào)面積都得到訂正,BIAS評(píng)分更接近1,說明訂正使原始集合平均的“有雨或無雨”的定性降水預(yù)報(bào),以及大暴雨量級(jí)預(yù)報(bào)偏小的現(xiàn)象有了改善。

    圖7 2018年5月8日12時(shí)(UTC)起報(bào)的12-36 h時(shí)段的24 h累計(jì)降水預(yù)報(bào)的GSS和BIAS評(píng)分藍(lán)色為原始集合平均預(yù)報(bào),紅色為FM訂正預(yù)報(bào)

    2.3 連續(xù)試驗(yàn)

    本文對(duì)2018年汛期4-9月103次降水個(gè)例進(jìn)行了連續(xù)試驗(yàn),并針對(duì)0.1 mm,1 mm,10 mm,25 mm,50 mm和100 mm等不同降水閾值,檢驗(yàn)了所有降水個(gè)例,以及前汛期(4-6月)和后汛期(7-9月)降水預(yù)報(bào)的原始集合平均預(yù)報(bào)和頻率匹配訂正預(yù)報(bào)的GSS評(píng)分和BIAS評(píng)分(圖8),其中前汛期降水過程為61次,后汛期降水過程為42次。從汛期所有降水個(gè)例的BIAS評(píng)分來看(圖8a),原始集合平均預(yù)報(bào)對(duì)50 mm以下閾值存在較明顯的濕偏差, BIAS評(píng)分超過1.5,對(duì)10 mm閾值的最大BIAS超過1.9,但對(duì)100 mm閾值以上的BIAS較小,為0.95,表明在2018年汛期SZ-SSEF系統(tǒng)的原始集合平均對(duì)強(qiáng)降水有較好的預(yù)測(cè)能力。頻率匹配方法較好地訂正了原始集合平均預(yù)報(bào)50 mm以下閾值的濕偏差,BIAS評(píng)分得到明顯改善,更接近于1。從所有降水過程的GSS評(píng)分(圖8a)來看,頻率匹配訂正對(duì)0.1 mm和1 mm閾值的GSS評(píng)分有較明顯的改善,這可能是因?yàn)槟J綄?duì)范圍較大的小雨量級(jí)降水落區(qū)的位置預(yù)報(bào)較好,但濕偏差造成了較高的錯(cuò)警率,頻率匹配訂正消除了偏差,減少了大片虛報(bào)的小雨量級(jí)區(qū)域,降低了錯(cuò)警率,進(jìn)而改善了GSS評(píng)分。頻率匹配訂正對(duì)于10 mm至100 mm閾值的GSS評(píng)分影響不顯著。其中,對(duì)10 mm閾值的GSS評(píng)分略有改善,對(duì)25 mm和100 mm閾值的GSS評(píng)分稍有下降,這可能是因?yàn)轭l率匹配訂正使強(qiáng)降水面積收縮,高分辨率雷暴尺度集合預(yù)報(bào)系統(tǒng)的降水預(yù)報(bào)更加難以與觀測(cè)重疊,“錯(cuò)警”和“漏報(bào)”的雙重懲罰,導(dǎo)致GSS評(píng)分下降。李俊等[14]也指出,在原理上頻率匹配訂正法不能用來訂正雨區(qū)位置的偏差,在降水落區(qū)較為準(zhǔn)確的情形下,訂正效果較好;如果降水落區(qū)預(yù)報(bào)出現(xiàn)明顯偏差時(shí),使用該方法并不能改進(jìn)預(yù)報(bào)結(jié)果。

    華南季風(fēng)氣候特征極其顯著,每年降水集中在兩個(gè)時(shí)期,前汛期以鋒面降水和暖區(qū)降水為主,后汛期以熱帶氣旋、東風(fēng)波等熱帶系統(tǒng)降水為主[30-32]。前汛期和后汛期環(huán)流形勢(shì)、影響系統(tǒng)不同,造成的降水落區(qū)和數(shù)值模式的可預(yù)報(bào)性也有不同[33-34]。因此,本文在對(duì)2018年整個(gè)汛期降水過程檢驗(yàn)的基礎(chǔ)上,分別對(duì)前汛期和后汛期的降水預(yù)報(bào)進(jìn)行檢驗(yàn)和分析。對(duì)比圖8(b)和圖8(c)可以看出,對(duì)不同閾值的降水預(yù)報(bào),SZ-SSEF系統(tǒng)原始集合平均對(duì)前汛期降水預(yù)報(bào)的GSS評(píng)分低于后汛期的,且前汛期BIAS評(píng)分高于后汛期的;對(duì)0.1 mm至10 mm量級(jí)的降水預(yù)報(bào),前汛期和后汛期GSS評(píng)分都比較高,但隨著降水量級(jí)的增大,前汛期的GSS評(píng)分迅速下降,尤其是對(duì)50 mm和100 mm閾值的降水預(yù)報(bào),前汛期GSS評(píng)分分別下降到0.05和0.02,明顯低于后汛期的0.17和0.21。這表明SZ-SSEF系統(tǒng)對(duì)后汛期熱帶氣旋等天氣系統(tǒng)的強(qiáng)降水預(yù)報(bào)效果較好,后汛期降水預(yù)報(bào)濕偏差的程度低于前汛期的偏差程度。華南前汛期以暖區(qū)暴雨為主,暖區(qū)暴雨產(chǎn)生的降雨量一般遠(yuǎn)大于鋒面系統(tǒng)形成的降雨量,而目前數(shù)值模式對(duì)暖區(qū)暴雨的預(yù)報(bào)能力十分有限[35],這可能是SZ-SSEF系統(tǒng)對(duì)前汛期降水預(yù)報(bào)評(píng)分偏低的原因。頻率匹配方法訂正后,對(duì)前汛期和后汛期的0.1 mm至1 mm閾值降水預(yù)測(cè)的GSS評(píng)分都有明顯的改善,對(duì)10 mm閾值降水預(yù)報(bào)的GSS評(píng)分略有提升,對(duì)50 mm暴雨以上的GSS評(píng)分改善不明顯;對(duì)不同閾值降水預(yù)報(bào)的偏差也明顯降低,訂正后BIAS評(píng)分與1更接近。從整體評(píng)分來看,頻率匹配方法訂正對(duì)后汛期降水預(yù)測(cè)改進(jìn)的效果更為明顯,GSS評(píng)分提高幅度高于前汛期的。

    圖8 2018年4-9月汛期103次降水個(gè)例試驗(yàn)GSS和BIAS評(píng)分統(tǒng)計(jì)結(jié)果(a)為所有個(gè)例評(píng)分,(b)為前汛期個(gè)例評(píng)分,(c)為后汛期個(gè)例評(píng)分;藍(lán)色為原始集合平均預(yù)報(bào),紅色為FM訂正預(yù)報(bào)

    上述分析結(jié)果表明,SZ-SSEF系統(tǒng)對(duì)小雨量級(jí)的預(yù)報(bào)能力較好,但隨著降水量級(jí)增大,SZ-SSEF系統(tǒng)預(yù)報(bào)能力逐漸降低;對(duì)于不同時(shí)期降水,預(yù)報(bào)能力有所區(qū)別,對(duì)前汛期降水預(yù)報(bào)的GSS評(píng)分低于后汛期的,且前汛期BIAS評(píng)分高于后汛期的,降水預(yù)報(bào)濕偏差更明顯。頻率匹配方法訂正后,減小降水預(yù)報(bào)偏差BIAS評(píng)分,提高對(duì)0.1 mm至1 mm閾值降水的GSS評(píng)分,但對(duì)達(dá)到50 mm暴雨以上的GSS評(píng)分改進(jìn)不明顯,未來還需通過更多的試驗(yàn)樣本進(jìn)行檢驗(yàn),以獲得更為可靠的統(tǒng)計(jì)結(jié)果。

    3 結(jié)論與討論

    本文采用頻率匹配方法,對(duì)深圳雷暴尺度集合預(yù)報(bào)系統(tǒng)SZ-SSEF在2018年4-9月汛期103次降水個(gè)例的集合平均降水預(yù)報(bào)進(jìn)行了訂正預(yù)報(bào)試驗(yàn),以檢驗(yàn)頻率匹配方法在雷暴尺度集合降水預(yù)報(bào)中的效果。總的來說,頻率匹配訂正法對(duì)改善雷暴尺度集合平均預(yù)報(bào)的降水強(qiáng)度分布有顯著效果,其具體表現(xiàn)如下:

    (1)深圳雷暴尺度集合預(yù)報(bào)系統(tǒng)對(duì)小雨量級(jí)的預(yù)報(bào)能力較好,但隨著降水量級(jí)增大,預(yù)報(bào)能力逐漸降低,對(duì)2018年前汛期降水預(yù)報(bào)的GSS評(píng)分低于后汛期的,且前汛期BIAS評(píng)分高于后汛期的,降水預(yù)報(bào)濕偏差更明顯。

    (2)頻率匹配訂正可以較好地訂正原始集合平均預(yù)報(bào)對(duì)降水量的過報(bào)情況,特別是對(duì)50 mm以下閾值降水預(yù)報(bào)的BIAS評(píng)分有明顯改善,較為接近1的完美評(píng)分。

    (3)對(duì)于低閾值降水預(yù)報(bào)來說,頻率匹配訂正可以將集合平均預(yù)報(bào)的降水區(qū)域范圍明顯減小,面積縮小到與融合分析資料相似的大小,減小了原始集合平均預(yù)報(bào)的錯(cuò)警率,進(jìn)而改善了GSS評(píng)分。

    (4)對(duì)于50 mm暴雨以上強(qiáng)降水預(yù)報(bào)來說,無論單個(gè)集合成員還是集合平均預(yù)報(bào),都很難精確捕捉到強(qiáng)降水中心的位置,做到格點(diǎn)對(duì)格點(diǎn)的重合。頻率匹配訂正對(duì)降水的位置偏差沒有訂正能力,同時(shí)會(huì)縮小強(qiáng)降水區(qū)域的面積。這可能會(huì)導(dǎo)致訂正預(yù)報(bào)的強(qiáng)降水區(qū)域更加難以與分析資料重合,出現(xiàn)“錯(cuò)警”和“漏報(bào)”,進(jìn)而降低GSS評(píng)分。這也說明,GSS等點(diǎn)對(duì)點(diǎn)評(píng)分方法對(duì)強(qiáng)降水預(yù)報(bào)的檢驗(yàn)存在不足,未來需要嘗試其他方法,將強(qiáng)降水的強(qiáng)度、位置誤差分開檢驗(yàn)。

    關(guān)于頻率匹配方法對(duì)強(qiáng)降水預(yù)報(bào)(如100 mm閾值)的訂正效果,本文所選取的2018年汛期103個(gè)降水個(gè)例中,強(qiáng)降水事件較為有限,其結(jié)果可能存在一定的偶然性。未來還需要進(jìn)一步擴(kuò)大樣本容量,或者單獨(dú)挑選強(qiáng)降水個(gè)例進(jìn)行檢驗(yàn),獲取更加全面、可靠的統(tǒng)計(jì)結(jié)果。

    猜你喜歡
    個(gè)例雷暴強(qiáng)降水
    新德里雷暴
    2020年江淮地區(qū)夏季持續(xù)性強(qiáng)降水過程分析
    一例膀胱鏡下留置尿管困難拔除個(gè)例分析
    一次東移型西南低渦引發(fā)的強(qiáng)降水診斷分析
    阜新地區(qū)雷暴活動(dòng)特點(diǎn)研究
    綠色科技(2018年22期)2019-01-15 05:36:08
    廣西富川縣雷暴氣候特征分析
    四川省春季一次飛機(jī)增雨作業(yè)個(gè)例分析
    四川盆地南部一次冬季飛機(jī)增雨作業(yè)個(gè)例分析
    一次森林滅火人工增雪個(gè)例分析
    西藏科技(2016年8期)2016-09-26 09:00:53
    四川盆地西南部短時(shí)強(qiáng)降水天氣特征分析
    午夜福利在线在线| 嫩草影院入口| av天堂中文字幕网| 欧美一级a爱片免费观看看| 欧美一级毛片孕妇| tocl精华| 一个人观看的视频www高清免费观看| 精品久久久久久,| 成人国产综合亚洲| 亚洲av中文字字幕乱码综合| 校园春色视频在线观看| 操出白浆在线播放| 丰满乱子伦码专区| 日韩 欧美 亚洲 中文字幕| 悠悠久久av| 一级黄色大片毛片| 国产激情欧美一区二区| 18禁黄网站禁片午夜丰满| 欧美激情在线99| 欧美日韩中文字幕国产精品一区二区三区| 国内揄拍国产精品人妻在线| 老汉色∧v一级毛片| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清三级在线| 免费在线观看亚洲国产| 一本久久中文字幕| 男女视频在线观看网站免费| 男女视频在线观看网站免费| 91麻豆精品激情在线观看国产| 久久人人精品亚洲av| 此物有八面人人有两片| 国产午夜福利久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久久免费视频| 国产精品久久久久久久久免 | 可以在线观看毛片的网站| 精品久久久久久成人av| 很黄的视频免费| 最后的刺客免费高清国语| 最近视频中文字幕2019在线8| 午夜a级毛片| 久久久成人免费电影| 国产不卡一卡二| 欧美av亚洲av综合av国产av| 精品国产亚洲在线| 日本成人三级电影网站| 日本一本二区三区精品| 国产成人欧美在线观看| 午夜老司机福利剧场| 亚洲av一区综合| 99久久精品热视频| 久久国产乱子伦精品免费另类| 久久精品91蜜桃| 在线天堂最新版资源| 老司机深夜福利视频在线观看| 国产在线精品亚洲第一网站| 热99re8久久精品国产| 日韩亚洲欧美综合| 黄片大片在线免费观看| 免费看a级黄色片| 午夜日韩欧美国产| 国产又黄又爽又无遮挡在线| 亚洲第一欧美日韩一区二区三区| 国产中年淑女户外野战色| bbb黄色大片| 亚洲成av人片免费观看| 亚洲精品国产精品久久久不卡| 久久这里只有精品中国| 亚洲欧美日韩无卡精品| 欧美+亚洲+日韩+国产| 国产午夜福利久久久久久| 成年版毛片免费区| 9191精品国产免费久久| 免费观看精品视频网站| 丰满人妻熟妇乱又伦精品不卡| 欧美3d第一页| 性色avwww在线观看| 舔av片在线| 亚洲精品国产精品久久久不卡| 色综合亚洲欧美另类图片| 少妇丰满av| 97人妻精品一区二区三区麻豆| 日本免费一区二区三区高清不卡| 我的老师免费观看完整版| 亚洲人成网站在线播放欧美日韩| 成人国产综合亚洲| 国产精品国产高清国产av| 欧美日韩国产亚洲二区| 91在线观看av| 亚洲欧美日韩东京热| 香蕉久久夜色| 国产亚洲精品一区二区www| 欧美+日韩+精品| 搡老熟女国产l中国老女人| 人人妻人人看人人澡| 国产亚洲精品综合一区在线观看| 久久国产精品人妻蜜桃| 亚洲国产精品久久男人天堂| 欧美3d第一页| 免费看日本二区| 国产麻豆成人av免费视频| 一进一出抽搐动态| 美女cb高潮喷水在线观看| 久9热在线精品视频| 青草久久国产| 国产成人aa在线观看| 午夜激情福利司机影院| 日本撒尿小便嘘嘘汇集6| 狂野欧美白嫩少妇大欣赏| tocl精华| 成人av在线播放网站| 黄色片一级片一级黄色片| 国产成人av教育| 真人一进一出gif抽搐免费| 最近在线观看免费完整版| 给我免费播放毛片高清在线观看| 国产欧美日韩一区二区精品| 成人av在线播放网站| 精品欧美国产一区二区三| 精品欧美国产一区二区三| 两人在一起打扑克的视频| 欧美zozozo另类| 国产高清有码在线观看视频| 女人十人毛片免费观看3o分钟| 两个人视频免费观看高清| 日本 欧美在线| 亚洲最大成人手机在线| 成人av在线播放网站| 国产单亲对白刺激| 亚洲avbb在线观看| 午夜久久久久精精品| 国产精品野战在线观看| 亚洲国产精品合色在线| 麻豆成人av在线观看| 又粗又爽又猛毛片免费看| 亚洲av成人不卡在线观看播放网| 美女免费视频网站| 他把我摸到了高潮在线观看| 岛国视频午夜一区免费看| 中文字幕人妻丝袜一区二区| 日韩欧美免费精品| 男女之事视频高清在线观看| 欧美成人a在线观看| www.www免费av| 一本一本综合久久| 女人被狂操c到高潮| 中文字幕人妻熟人妻熟丝袜美 | 国产精品影院久久| 神马国产精品三级电影在线观看| 啦啦啦免费观看视频1| 日韩欧美一区二区三区在线观看| 国产精品一区二区三区四区久久| 欧美成人a在线观看| 国产69精品久久久久777片| 天堂av国产一区二区熟女人妻| 国内精品一区二区在线观看| 成人特级黄色片久久久久久久| 日本黄大片高清| 久久久久久九九精品二区国产| 99久久久亚洲精品蜜臀av| 伊人久久大香线蕉亚洲五| 国产三级在线视频| 免费电影在线观看免费观看| 久久久国产精品麻豆| 午夜亚洲福利在线播放| 亚洲aⅴ乱码一区二区在线播放| 毛片女人毛片| 国产中年淑女户外野战色| 亚洲第一电影网av| 99久久99久久久精品蜜桃| 日韩国内少妇激情av| 一区二区三区激情视频| 精品一区二区三区视频在线观看免费| 精品人妻1区二区| 五月伊人婷婷丁香| 久久99热这里只有精品18| 国产一区二区亚洲精品在线观看| 国产伦人伦偷精品视频| 在线观看午夜福利视频| 亚洲精品久久国产高清桃花| 18禁黄网站禁片午夜丰满| 亚洲一区高清亚洲精品| 久久国产精品人妻蜜桃| 欧美日韩国产亚洲二区| www日本在线高清视频| 一级黄色大片毛片| 国产乱人视频| 久久久成人免费电影| 露出奶头的视频| 听说在线观看完整版免费高清| 亚洲欧美一区二区三区黑人| 最后的刺客免费高清国语| 日韩有码中文字幕| 欧美中文综合在线视频| 午夜免费男女啪啪视频观看 | 在线观看66精品国产| 日韩欧美精品v在线| 国产黄a三级三级三级人| 亚洲国产精品成人综合色| 十八禁人妻一区二区| 精品欧美国产一区二区三| 亚洲人成网站在线播| 国产亚洲精品综合一区在线观看| 日本黄色视频三级网站网址| 中文字幕av成人在线电影| 国产伦精品一区二区三区视频9 | 欧美成人a在线观看| 操出白浆在线播放| 亚洲av电影不卡..在线观看| 小蜜桃在线观看免费完整版高清| 欧美激情久久久久久爽电影| 精品人妻偷拍中文字幕| 波多野结衣高清作品| 国产精品久久久久久人妻精品电影| 好男人在线观看高清免费视频| 一本精品99久久精品77| 99久久九九国产精品国产免费| 最好的美女福利视频网| 长腿黑丝高跟| 俺也久久电影网| 亚洲精品456在线播放app | 有码 亚洲区| 亚洲精华国产精华精| 国产69精品久久久久777片| 人妻丰满熟妇av一区二区三区| 女同久久另类99精品国产91| 中亚洲国语对白在线视频| 午夜免费激情av| 夜夜夜夜夜久久久久| 人妻丰满熟妇av一区二区三区| 少妇熟女aⅴ在线视频| 波多野结衣高清无吗| 久久中文看片网| 男女视频在线观看网站免费| 美女高潮的动态| av专区在线播放| av天堂中文字幕网| 男女午夜视频在线观看| 国产精品爽爽va在线观看网站| 久久久色成人| 午夜久久久久精精品| 色老头精品视频在线观看| 亚洲 国产 在线| 日韩成人在线观看一区二区三区| 国产老妇女一区| 级片在线观看| 哪里可以看免费的av片| 99国产极品粉嫩在线观看| 免费看日本二区| 亚洲av熟女| 69av精品久久久久久| 欧美中文日本在线观看视频| 9191精品国产免费久久| 国产爱豆传媒在线观看| 久久中文看片网| 啪啪无遮挡十八禁网站| 欧美中文日本在线观看视频| 一区二区三区国产精品乱码| 国产高潮美女av| 亚洲精品粉嫩美女一区| 三级毛片av免费| 在线十欧美十亚洲十日本专区| 欧美av亚洲av综合av国产av| 综合色av麻豆| 亚洲人成电影免费在线| 午夜精品一区二区三区免费看| 亚洲电影在线观看av| 国产高清激情床上av| 欧美黑人巨大hd| 国产精品久久久久久久久免 | 一级黄片播放器| 国产综合懂色| 欧美日韩精品网址| 亚洲成人中文字幕在线播放| 嫩草影院精品99| 亚洲人成伊人成综合网2020| 可以在线观看毛片的网站| 在线a可以看的网站| 黄色女人牲交| 欧美成人a在线观看| 老司机午夜福利在线观看视频| 一个人看视频在线观看www免费 | 熟女人妻精品中文字幕| 亚洲av中文字字幕乱码综合| 男人舔奶头视频| 久久精品综合一区二区三区| 国产精华一区二区三区| 亚洲av电影不卡..在线观看| 老司机在亚洲福利影院| 国产真人三级小视频在线观看| 18禁在线播放成人免费| 非洲黑人性xxxx精品又粗又长| 国产激情偷乱视频一区二区| 色哟哟哟哟哟哟| 美女高潮喷水抽搐中文字幕| 国产欧美日韩一区二区精品| 最新美女视频免费是黄的| 亚洲色图av天堂| 少妇人妻一区二区三区视频| 免费大片18禁| www国产在线视频色| 欧美精品啪啪一区二区三区| 国产男靠女视频免费网站| 欧美日韩中文字幕国产精品一区二区三区| 精品欧美国产一区二区三| 精品一区二区三区av网在线观看| 麻豆成人av在线观看| 九九久久精品国产亚洲av麻豆| 亚洲黑人精品在线| 女警被强在线播放| 在线播放无遮挡| 人人妻人人澡欧美一区二区| 亚洲性夜色夜夜综合| 三级毛片av免费| 欧美乱色亚洲激情| 波多野结衣高清无吗| 国产精品一区二区三区四区久久| 亚洲精品久久国产高清桃花| 国产aⅴ精品一区二区三区波| 日韩中文字幕欧美一区二区| 88av欧美| 夜夜夜夜夜久久久久| 在线十欧美十亚洲十日本专区| 午夜影院日韩av| 美女黄网站色视频| 嫁个100分男人电影在线观看| 午夜久久久久精精品| 国产老妇女一区| 一进一出好大好爽视频| 一区二区三区免费毛片| 午夜久久久久精精品| 久久精品91蜜桃| 国产一区二区在线av高清观看| 精品一区二区三区视频在线 | 欧美bdsm另类| 一个人观看的视频www高清免费观看| 日本在线视频免费播放| h日本视频在线播放| 日韩欧美一区二区三区在线观看| 色精品久久人妻99蜜桃| 亚洲在线观看片| 1000部很黄的大片| 色综合站精品国产| 欧美丝袜亚洲另类 | 一夜夜www| 黄色丝袜av网址大全| 欧美+日韩+精品| 看免费av毛片| 波多野结衣高清无吗| 一个人免费在线观看电影| 特大巨黑吊av在线直播| 欧美日韩精品网址| 宅男免费午夜| 麻豆国产av国片精品| 欧美日本视频| avwww免费| 久久午夜亚洲精品久久| 午夜视频国产福利| 亚洲国产精品sss在线观看| 夜夜爽天天搞| 精品国内亚洲2022精品成人| 日韩精品青青久久久久久| 国内毛片毛片毛片毛片毛片| 亚洲久久久久久中文字幕| 熟女电影av网| 精品久久久久久久毛片微露脸| 亚洲欧美精品综合久久99| 99视频精品全部免费 在线| 免费高清视频大片| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看 | 亚洲人与动物交配视频| 一本综合久久免费| 久久婷婷人人爽人人干人人爱| 久久久国产成人免费| 国产精品99久久久久久久久| 国产精品,欧美在线| 免费在线观看成人毛片| 网址你懂的国产日韩在线| 国产成人福利小说| 老司机午夜福利在线观看视频| 久久精品91蜜桃| 舔av片在线| 欧美中文综合在线视频| 人妻夜夜爽99麻豆av| 亚洲av第一区精品v没综合| 中文字幕人妻熟人妻熟丝袜美 | 亚洲人成伊人成综合网2020| 黄色女人牲交| 麻豆国产av国片精品| 内射极品少妇av片p| 欧美午夜高清在线| 亚洲人与动物交配视频| 色老头精品视频在线观看| 亚洲熟妇中文字幕五十中出| 成年女人毛片免费观看观看9| 成人特级黄色片久久久久久久| 丰满乱子伦码专区| 午夜福利成人在线免费观看| 成年女人看的毛片在线观看| 老汉色av国产亚洲站长工具| 国产精华一区二区三区| 日韩欧美在线乱码| av欧美777| 成人性生交大片免费视频hd| 亚洲一区二区三区不卡视频| 国产精品香港三级国产av潘金莲| 色精品久久人妻99蜜桃| 国产av麻豆久久久久久久| 国产野战对白在线观看| 岛国在线免费视频观看| 久久国产精品人妻蜜桃| 欧美日本亚洲视频在线播放| 在线观看午夜福利视频| 一级毛片高清免费大全| 亚洲片人在线观看| 亚洲欧美日韩高清专用| 欧美乱码精品一区二区三区| 国产精品久久久久久精品电影| 在线观看日韩欧美| 欧美精品啪啪一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 色综合欧美亚洲国产小说| 一级黄色大片毛片| 噜噜噜噜噜久久久久久91| 黄色丝袜av网址大全| 每晚都被弄得嗷嗷叫到高潮| 香蕉丝袜av| 久久人妻av系列| 午夜影院日韩av| 美女cb高潮喷水在线观看| netflix在线观看网站| 久久精品国产清高在天天线| 老司机午夜十八禁免费视频| 女同久久另类99精品国产91| 俄罗斯特黄特色一大片| 亚洲av成人精品一区久久| 男女午夜视频在线观看| 美女cb高潮喷水在线观看| 特级一级黄色大片| 亚洲国产欧洲综合997久久,| 亚洲无线观看免费| 9191精品国产免费久久| 男女做爰动态图高潮gif福利片| 国产麻豆成人av免费视频| 中文字幕av在线有码专区| 精品国产三级普通话版| 久9热在线精品视频| 免费看光身美女| 变态另类丝袜制服| 人妻久久中文字幕网| 午夜福利高清视频| 午夜福利在线观看吧| 国产探花在线观看一区二区| 欧美bdsm另类| 欧美色视频一区免费| 无限看片的www在线观看| 精品乱码久久久久久99久播| 很黄的视频免费| 在线天堂最新版资源| 亚洲av一区综合| 欧美一级a爱片免费观看看| 91av网一区二区| 女警被强在线播放| 一级黄色大片毛片| 最近视频中文字幕2019在线8| 99久久成人亚洲精品观看| 日韩精品青青久久久久久| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 美女黄网站色视频| 美女高潮的动态| 午夜福利高清视频| 91久久精品国产一区二区成人 | 国产色婷婷99| 国产真实伦视频高清在线观看 | 免费在线观看影片大全网站| 日本a在线网址| 精品国内亚洲2022精品成人| 久久精品夜夜夜夜夜久久蜜豆| 51国产日韩欧美| 亚洲国产高清在线一区二区三| 日韩欧美免费精品| 国产成人av激情在线播放| 亚洲欧美日韩无卡精品| 成人性生交大片免费视频hd| 日韩精品中文字幕看吧| 国产亚洲欧美在线一区二区| svipshipincom国产片| 欧美av亚洲av综合av国产av| 在线观看日韩欧美| 久久国产精品影院| 欧美一区二区亚洲| 亚洲欧美日韩高清在线视频| 久久亚洲精品不卡| 欧美丝袜亚洲另类 | 久久久色成人| 国产黄片美女视频| 日韩欧美免费精品| 中文字幕高清在线视频| 精品乱码久久久久久99久播| 日本三级黄在线观看| 少妇人妻一区二区三区视频| 国产中年淑女户外野战色| 日本黄色视频三级网站网址| 天天躁日日操中文字幕| 99国产精品一区二区三区| 久久久精品欧美日韩精品| 99热这里只有是精品50| 久久这里只有精品中国| 色av中文字幕| 美女大奶头视频| 国内精品久久久久久久电影| 精品人妻偷拍中文字幕| 久9热在线精品视频| 国产99白浆流出| 久久久久国内视频| 成人国产综合亚洲| 日韩av在线大香蕉| 婷婷亚洲欧美| 69人妻影院| 久久精品国产亚洲av涩爱 | 国产一区二区三区视频了| 国产亚洲av嫩草精品影院| 亚洲中文字幕日韩| 99riav亚洲国产免费| 一级黄色大片毛片| 中国美女看黄片| 国产精品99久久99久久久不卡| 日韩av在线大香蕉| 欧美日韩精品网址| 国产精品久久视频播放| 波野结衣二区三区在线 | 成人18禁在线播放| 午夜福利欧美成人| 亚洲一区高清亚洲精品| 村上凉子中文字幕在线| 午夜精品久久久久久毛片777| 亚洲精品在线美女| 日韩欧美精品免费久久 | 国产视频一区二区在线看| 亚洲欧美日韩卡通动漫| 少妇熟女aⅴ在线视频| 国产黄色小视频在线观看| 香蕉久久夜色| 我要搜黄色片| 免费av观看视频| 中文字幕熟女人妻在线| 色综合婷婷激情| or卡值多少钱| 韩国av一区二区三区四区| 18禁在线播放成人免费| 欧美日韩国产亚洲二区| 欧美一区二区亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 偷拍熟女少妇极品色| 久久精品国产亚洲av涩爱 | 熟女电影av网| 久久天躁狠狠躁夜夜2o2o| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美免费精品| 欧美乱色亚洲激情| 国产高清有码在线观看视频| 午夜日韩欧美国产| 极品教师在线免费播放| 久久久久免费精品人妻一区二区| 亚洲av电影在线进入| 午夜免费男女啪啪视频观看 | 最后的刺客免费高清国语| 日韩精品中文字幕看吧| 免费人成视频x8x8入口观看| 亚洲av二区三区四区| 欧美xxxx黑人xx丫x性爽| 亚洲成人免费电影在线观看| 午夜精品在线福利| 免费在线观看影片大全网站| 可以在线观看的亚洲视频| 久久精品人妻少妇| 色视频www国产| 亚洲成人中文字幕在线播放| 蜜桃亚洲精品一区二区三区| 欧美在线黄色| 成人亚洲精品av一区二区| 日韩 欧美 亚洲 中文字幕| 欧美日韩亚洲国产一区二区在线观看| 在线观看一区二区三区| svipshipincom国产片| 宅男免费午夜| 男女那种视频在线观看| 久久精品影院6| 小蜜桃在线观看免费完整版高清| 极品教师在线免费播放| 久久精品国产亚洲av涩爱 | 免费在线观看影片大全网站| 亚洲人与动物交配视频| 午夜免费成人在线视频| 色综合欧美亚洲国产小说| 白带黄色成豆腐渣| 啪啪无遮挡十八禁网站| 婷婷亚洲欧美| 国产精品久久视频播放| 亚洲国产高清在线一区二区三| 麻豆国产av国片精品| 久久精品人妻少妇| 亚洲精品456在线播放app | 国产成人av激情在线播放| 亚洲国产欧美人成| 91九色精品人成在线观看| 搡老妇女老女人老熟妇| 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 午夜福利高清视频| 夜夜夜夜夜久久久久| 国产一区二区在线观看日韩 | 在线天堂最新版资源|