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

    ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品在烏魯木齊市溫度預(yù)報(bào)中的檢驗(yàn)

    2015-01-04 09:52:40張超李娜賈健
    沙漠與綠洲氣象 2015年5期
    關(guān)鍵詞:格點(diǎn)最低氣溫準(zhǔn)確率

    張超,李娜,賈健

    (1.新疆氣象臺(tái),新疆烏魯木齊 830002;2.烏魯木齊市氣象局,新疆烏魯木齊 830002)

    ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品在烏魯木齊市溫度預(yù)報(bào)中的檢驗(yàn)

    張超1,李娜1,賈健2

    (1.新疆氣象臺(tái),新疆烏魯木齊 830002;2.烏魯木齊市氣象局,新疆烏魯木齊 830002)

    利用2011年7月—2014年6月的ECMWF細(xì)網(wǎng)格模式2m溫度預(yù)報(bào)產(chǎn)品,采用格點(diǎn)映射站點(diǎn)和雙線性插值法挑選最優(yōu)預(yù)報(bào)參考格點(diǎn),并用天氣學(xué)檢驗(yàn)方法,對(duì)該模式在不同季節(jié)、不同天氣時(shí)的烏魯木齊市溫度預(yù)報(bào)能力進(jìn)行檢驗(yàn)。檢驗(yàn)表明:該產(chǎn)品對(duì)制作烏魯木齊市未來24h和48h溫度預(yù)報(bào)具有很好的指導(dǎo)意義,未來24h全年最高、最低氣溫預(yù)報(bào)準(zhǔn)確率達(dá)74%~75%,未來48h最高、最低氣溫預(yù)報(bào)準(zhǔn)確率較未來24h偏低10%左右;模式對(duì)降雨、降雪時(shí)的溫度預(yù)報(bào)較好,大霧時(shí)最差,并且會(huì)使最低氣溫預(yù)報(bào)嚴(yán)重低于實(shí)況,對(duì)東南風(fēng)時(shí)最高氣溫預(yù)報(bào)好于最低氣溫,預(yù)報(bào)有逆溫時(shí),最高氣溫預(yù)報(bào)不可信,但最低氣溫預(yù)報(bào)可作參考。該產(chǎn)品整體在準(zhǔn)確率和穩(wěn)定性上均優(yōu)于國家氣象中心指導(dǎo)預(yù)報(bào),預(yù)報(bào)時(shí)效越長,優(yōu)勢(shì)越明顯。

    ECMWF細(xì)網(wǎng)格;溫度檢驗(yàn);2m溫度;烏魯木齊

    近年來,越來越多的氣象工作者通過數(shù)值預(yù)報(bào)模式產(chǎn)品檢驗(yàn)釋用提高預(yù)報(bào)準(zhǔn)確率和精細(xì)化水平,取得了較好的效果[1-9]。新疆接收ECMWF細(xì)網(wǎng)格模式預(yù)報(bào)產(chǎn)品已近4a,其在很多方面都優(yōu)于ECMWF粗網(wǎng)格模式。國內(nèi)一些氣象工作者利用該模式產(chǎn)品針對(duì)大霧、降水、溫度預(yù)報(bào)等方面做了檢驗(yàn)和釋用[10-13]。在新疆的ECMWF細(xì)網(wǎng)格模式預(yù)報(bào)產(chǎn)品的釋用檢驗(yàn)工作中,張俊蘭等[15]檢驗(yàn)了ECMWF細(xì)網(wǎng)格模式的LSP大尺度降水預(yù)報(bào)產(chǎn)品在2012年前冬3場(chǎng)暴雪天氣中的預(yù)報(bào)效果,萬瑜等[16]利用ECMWF細(xì)網(wǎng)格模式產(chǎn)品對(duì)一次烏魯木齊東南大風(fēng)做了釋用,都顯示出了ECMWF細(xì)網(wǎng)格模式產(chǎn)品出色的預(yù)報(bào)性能,大大提高了各類天氣預(yù)報(bào)的精細(xì)化程度,但這些僅是針對(duì)某次天氣過程的研究,缺乏對(duì)長時(shí)間序列預(yù)報(bào)效果的檢驗(yàn)。

    ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品的空間分辨率為25km,為烏魯木齊市單站預(yù)報(bào)提供更多的參考格點(diǎn),具有未來72h內(nèi)時(shí)間間隔3h的溫度預(yù)報(bào)產(chǎn)品,未來12~36h的最高、最低值可選取作為該模式次日未來24h的最高、最低氣溫預(yù)報(bào),未來36~60h的最高、最低值可選取作為該模式次日未來48h的最高、最低氣溫預(yù)報(bào),但此產(chǎn)品對(duì)烏魯木齊市溫度的預(yù)報(bào)效果如何,最佳參考格點(diǎn)如何選取并未進(jìn)行過客觀的檢驗(yàn)分析。由于烏魯木齊地形復(fù)雜,周邊模式格點(diǎn)所在海拔高度、下墊面性質(zhì)差異大,模式各格點(diǎn)預(yù)報(bào)數(shù)值差異也較大,需獲取能代表烏魯木齊溫度的最優(yōu)預(yù)報(bào)值。本文對(duì)2011年7月—2014年6月ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品資料未來24h及48h溫度預(yù)報(bào)進(jìn)行檢驗(yàn)評(píng)估,結(jié)合格點(diǎn)映射站點(diǎn)和雙線性插值2種方法獲取烏魯木齊溫度預(yù)報(bào)值,選出代表烏魯木齊溫度預(yù)報(bào)的最優(yōu)值,統(tǒng)計(jì)和歸納不同季節(jié)、不同天氣現(xiàn)象下ECMWF細(xì)網(wǎng)格模式對(duì)溫度預(yù)報(bào)的性能。此工作將為制作準(zhǔn)確、精細(xì)的溫度預(yù)報(bào)提供參考,并對(duì)做好復(fù)雜天氣背景下溫度預(yù)報(bào)具有重要意義。

    1 資料與方法

    1.1 資料

    實(shí)況溫度資料選用2011年7月—2014年6月共3a烏魯木齊市地面自動(dòng)觀測(cè)站逐小時(shí)溫度資料;ECMWF細(xì)網(wǎng)格模式資料采用2011年7月—2014年6月20h起報(bào)2m溫度預(yù)報(bào)資料,空間分辨率為0.25°×0.25°,時(shí)間分辨率為3h,數(shù)據(jù)缺失日不檢驗(yàn),共有1047d的有效檢驗(yàn)樣本;為同國家氣象中心指導(dǎo)預(yù)報(bào)進(jìn)行對(duì)比檢驗(yàn),還選用了2013年1—12月國家氣象中心下發(fā)的最高、最低氣溫指導(dǎo)預(yù)報(bào)。

    1.2 方法

    提取08—08時(shí)內(nèi)地面自動(dòng)觀測(cè)站逐小時(shí)溫度最高、最低值作為實(shí)況的最高、最低氣溫,采用格點(diǎn)映射站點(diǎn)和雙線性插值法提取預(yù)報(bào)檢驗(yàn)值,并將模式預(yù)報(bào)未來12~36h和36~60h的溫度最高、最低值,作為ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品未來24h和48h的最高、最低氣溫預(yù)報(bào)。

    格點(diǎn)映射站點(diǎn)方法:ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品為格點(diǎn)數(shù)據(jù),烏魯木齊站不在格點(diǎn)上,周邊存在4個(gè)預(yù)報(bào)值,分別位于烏魯木齊左上、左下、右上、右下4個(gè)方位,用以上4個(gè)方位格點(diǎn)數(shù)據(jù)映射到烏魯木齊站點(diǎn)上,用T左上、T左下、T右上、T右下格點(diǎn)映射數(shù)據(jù)代表ECMWF細(xì)網(wǎng)格模式2m溫度資料對(duì)烏魯木齊的預(yù)報(bào)值。

    雙線性插值法:又稱為雙線性內(nèi)插,在數(shù)學(xué)上,雙線性插值是有兩個(gè)變量的插值函數(shù)的線性插值擴(kuò)展,其核心思想是在兩個(gè)方向分別進(jìn)行一次線性插值。若未知函數(shù)f在點(diǎn)P=(x,y)的值,假設(shè)我們已知函數(shù)f在Q11=(x1,y1)、Q12=(x1,y2),Q21=(x2,y1)以及Q22=(x2,y2)四個(gè)點(diǎn)的值。首先在x方向進(jìn)行線性插值,然后在y方向進(jìn)行線性插值。本文將模式的格點(diǎn)映射值通過雙線性插值法插值到烏魯木齊站點(diǎn)上,得到T差值。

    檢驗(yàn)方法:采用預(yù)報(bào)準(zhǔn)確率、平均誤差、平均絕對(duì)誤差等統(tǒng)計(jì)方法檢驗(yàn)(公式如下),對(duì)烏魯木齊左上(T左上)、左下(T左下)、右上(T右上)、右下(T右下)及插值(T差值)后的ECMWF細(xì)網(wǎng)格模式的2m溫度預(yù)報(bào)值按月、季、天氣類型分別進(jìn)行檢驗(yàn)。

    其中,F(xiàn)i為第i站(次)預(yù)報(bào)溫度,Oi為第i站(次)實(shí)況溫度,K為1、2,分別代表|Fi-Oi|≤1℃、|Fi-Oi|≤2℃,MrK為預(yù)報(bào)正確的站(次)數(shù),MfK為預(yù)報(bào)的總站(次)數(shù)。

    評(píng)分技巧是指ECMWF細(xì)網(wǎng)格模式2m溫度預(yù)報(bào)產(chǎn)品的最高、最低氣溫預(yù)報(bào)平均絕對(duì)誤差相對(duì)于國家氣象中心指導(dǎo)預(yù)報(bào)的提高率,TMAEN為國家氣象中心指導(dǎo)預(yù)報(bào)最高、最低氣溫預(yù)報(bào)平均絕對(duì)誤差,TMAEF為ECMWF細(xì)網(wǎng)格2m溫度預(yù)報(bào)產(chǎn)品的最高、最低氣溫預(yù)報(bào)平均絕對(duì)誤差。

    2 預(yù)報(bào)檢驗(yàn)結(jié)果

    2.1 高、低溫預(yù)報(bào)逐月檢驗(yàn)

    根據(jù)中國氣象局對(duì)各省、地(市)氣象臺(tái)的考核標(biāo)準(zhǔn),用誤差為2℃以內(nèi)作為溫度預(yù)報(bào)準(zhǔn)確率的閾值。圖1為2011年7月—2014年6月ECMWF細(xì)網(wǎng)格模式2m溫度的未來24h最高、最低溫度預(yù)報(bào)準(zhǔn)確率。圖1a中,最高氣溫預(yù)報(bào)準(zhǔn)確率從高到低排列分別為T左下、T差值、T右上、T右下、T左上,左下格點(diǎn)T左下夏半年的最高氣溫預(yù)報(bào)準(zhǔn)確率最高,尤其是5—8月最高氣溫預(yù)報(bào)準(zhǔn)確率均超過了85%,7月準(zhǔn)確率最高,達(dá)到92.3%;右上格點(diǎn)溫度預(yù)報(bào)T右上的準(zhǔn)確率呈現(xiàn)出“雙峰型”,在春季和秋季準(zhǔn)確率明顯高于其他時(shí)期,4月、10月超過了T左下的準(zhǔn)確率,達(dá)到了84%;左上格點(diǎn)T左上的準(zhǔn)確率最低,5—8月不到50%,6、7月最低,不到20%;雙線性插值預(yù)報(bào)的準(zhǔn)確率均低于T左下。可見不同格點(diǎn)的預(yù)報(bào)能力存在明顯差異,插值的溫度預(yù)報(bào)T差值并沒有預(yù)期的優(yōu)越性,在春、夏、秋季均低于T左下的最高氣溫預(yù)報(bào)準(zhǔn)確率。冬季各格點(diǎn)預(yù)報(bào)能力普遍較差,預(yù)報(bào)準(zhǔn)確率都偏低,在30%~60%之間。

    由圖1b可知各參考點(diǎn)的預(yù)報(bào)能力具有更為明顯的季節(jié)性變化,T右上的最低氣溫預(yù)報(bào)準(zhǔn)確率同樣呈現(xiàn)出“雙峰型”,在3—5月和9—1月有著明顯的優(yōu)勢(shì),4月、5月、9月、10月的預(yù)報(bào)準(zhǔn)確率都超過了80%,其中10月的最低氣溫預(yù)報(bào)準(zhǔn)確率達(dá)到了91.8%,6月開始,T右上的最低氣溫預(yù)報(bào)準(zhǔn)確率開始明顯下降,7月降到了54.9%;而T差值預(yù)報(bào)在6月、7月表現(xiàn)較好,準(zhǔn)確率分別達(dá)到了69.7%、74.7%,對(duì)T右上夏季的不足有了較好的彌補(bǔ);在所有參考點(diǎn)中,T右下的最低氣溫預(yù)報(bào)能力最差,對(duì)烏市單站的最低氣溫預(yù)報(bào)缺乏參考意義;在12月—翌年2月,各格點(diǎn)最低氣溫預(yù)報(bào)準(zhǔn)確率也較其他月份差,準(zhǔn)確率在20%~60%,2月,T左上的最低氣溫預(yù)報(bào)準(zhǔn)確率最高,可達(dá)到59.2%,12月和1月的最低氣溫預(yù)報(bào)準(zhǔn)確率仍然是T右上最高,準(zhǔn)確率分別達(dá)到71.4%、55.2%。

    圖1 2011年7月—2014年6月ECMWF細(xì)網(wǎng)格2m溫度對(duì)烏魯木齊溫度未來24h預(yù)報(bào)準(zhǔn)確率/%

    圖2為2011年7月—2014年6月ECMWF細(xì)網(wǎng)格模式2m溫度未來24h的最高、最低氣溫的誤差分析,檢驗(yàn)?zāi)J筋A(yù)報(bào)的平均誤差和平均絕對(duì)誤差,其中平均誤差能夠量化預(yù)報(bào)和實(shí)況之間的平均偏差,也是模式預(yù)報(bào)的系統(tǒng)誤差,平均絕對(duì)誤差能夠量化預(yù)報(bào)與實(shí)況的總體偏差程度。由圖2a可以看出,4個(gè)格點(diǎn)最高氣溫預(yù)報(bào)平均誤差隨季節(jié)也有明顯變化,T右上和T左上在2—10月都為預(yù)報(bào)正誤差,即預(yù)報(bào)偏高誤差,在5—8月平均誤差都接近3℃,11—1月逐漸轉(zhuǎn)為負(fù)誤差;而T左下、T右下及T差值平均誤差月分布正好相反,在3—11月為負(fù)誤差,其中T左下的平均誤差最小,12—翌年2月逐漸轉(zhuǎn)為正誤差,開始轉(zhuǎn)換的時(shí)期正是冬春和秋冬的過渡時(shí)期。從圖2c來分析,T右下和T左上的預(yù)報(bào)值較實(shí)況偏離程度明顯高于其余參考點(diǎn),偏離值為2~3.5℃;T左下和T差值的預(yù)報(bào)值在3—11月較實(shí)況偏離程度較小,均低于1.6℃,但在11—2月偏離程度迅速加大,偏離值達(dá)2.5~3.5℃,T右上的平均絕對(duì)誤差呈現(xiàn)出了“雙谷型”,4月、10月為兩個(gè)低值,偏離值在1.5℃以下,其他時(shí)段在2~2.5℃之間。由圖2b可知,各格點(diǎn)最低氣溫的平均誤差季節(jié)分布規(guī)律與高溫類似,在夏半年T左下、T右下及T差值的溫度預(yù)報(bào)較實(shí)況略偏低,而在冬半年迅速轉(zhuǎn)為偏高,1月偏高幅度達(dá)到最強(qiáng),為3~4℃,轉(zhuǎn)換的時(shí)期也正是冬春和秋冬的過渡時(shí)期;T右上和T左上在5—8月溫度預(yù)報(bào)較實(shí)況有所偏高,其他月偏低,總體來看平均誤差在-2~2℃之間。由圖2d可以看出,在各月4個(gè)格點(diǎn)較實(shí)況偏離范圍在1~4.5℃之間,T右上的偏離程度最小,在1~2.5℃之間,冬季普遍高于其他季節(jié)。

    通過2011年7月—2014年6月ECMWF細(xì)網(wǎng)格模式2m溫度的未來48h最高、最低氣溫預(yù)報(bào)準(zhǔn)確率和預(yù)報(bào)誤差分析。發(fā)現(xiàn),各格點(diǎn)預(yù)報(bào)水平有同未來24h相似的逐月變化規(guī)律,依然具有很好的指導(dǎo)意義。全年最高、最低氣溫預(yù)報(bào)準(zhǔn)確率為65%左右,但較未來24h最高、最低氣溫預(yù)報(bào)準(zhǔn)確率偏低10%左右,其中春季偏低幅度最大,為10%~28%,夏秋季降低幅度最小,僅為1%~5%,這同春季模式調(diào)整頻繁有關(guān)。

    綜上所述,ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品對(duì)制作烏魯木齊未來24h和48h溫度都具有很好的指導(dǎo)意義,隨著季節(jié)的變化,不同參考點(diǎn)對(duì)烏魯木齊市高低溫預(yù)報(bào)能力也相應(yīng)變化,在冬春和秋冬的過渡時(shí)期,各格點(diǎn)的系統(tǒng)誤差會(huì)發(fā)生明顯變化,多數(shù)時(shí)間格點(diǎn)映射站點(diǎn)的方法要優(yōu)于雙線性插值法。經(jīng)過分析,可以把T左下的日最大值代表2月、3月、5月、6月、7月、8月和9月的最高氣溫預(yù)報(bào),T右上的日最大值代表1月、4月、10月、11月和12月的最高氣溫預(yù)報(bào),T右上的日最小值代表1月、3月、4月、5月、9月、10月、11月、12月的最低氣溫預(yù)報(bào),T差值的日最小值代表6月、7月的最低氣溫預(yù)報(bào)。T左上的日最小值代表2月、8月的最低氣溫預(yù)報(bào),未來24h全年最高、最低氣溫準(zhǔn)確率可以達(dá)到74%~75%,未來48h最高、最低氣溫預(yù)報(bào)準(zhǔn)確率較24h偏低10%左右,其中春季降低幅度最大,為10%~28%,夏秋季較小,為1%~5%。

    2.2 高、低溫預(yù)報(bào)季度檢驗(yàn)

    圖2 2011年7月—2014年6月ECMWF細(xì)網(wǎng)格2m溫度的未來24h預(yù)報(bào)誤差分析

    為了便于統(tǒng)計(jì)、分析,采用3—5月為春季、6—8月為夏季、9—11月為秋季、12月—翌年2月為冬季的季節(jié)劃分。重點(diǎn)分析ECMWF細(xì)網(wǎng)格模式2m溫度未來24h預(yù)報(bào)產(chǎn)品逐季節(jié)預(yù)報(bào)準(zhǔn)確率,發(fā)現(xiàn)有明顯的季節(jié)變化,夏半年的預(yù)報(bào)準(zhǔn)確率高出冬半10%~25%。最高氣溫準(zhǔn)確率在夏季最高,可達(dá)到88.7%,冬季最低,僅有53.3%,夏季高出冬季35%;而最低氣溫預(yù)報(bào)準(zhǔn)確率,秋季最高,為83.1%,冬季最低,準(zhǔn)確率為63.2%,秋季高出冬季20%,且最低氣溫預(yù)報(bào)準(zhǔn)確率在春、秋、冬季均高出同季節(jié)最高氣溫預(yù)報(bào)3%~9%,夏季,最高氣溫預(yù)報(bào)準(zhǔn)確率高出最低氣溫預(yù)報(bào)準(zhǔn)確率14%。最高氣溫預(yù)報(bào)在夏季偏離實(shí)況的程度最小,平均絕對(duì)誤差在1℃附近,冬季最大,平均絕對(duì)誤差接近2.5℃;最低氣溫預(yù)報(bào)在秋季偏離實(shí)況程度最小,冬季最大,并且最低氣溫預(yù)報(bào)的平均絕對(duì)誤差在春、秋、冬季均低于同季節(jié)最高氣溫預(yù)報(bào)的平均絕對(duì)誤差,在夏季高于最高氣溫預(yù)報(bào)的平均誤差。冬季最高、最低氣溫預(yù)報(bào)偏離實(shí)況程度最大,同時(shí)也是降低全年溫度預(yù)報(bào)準(zhǔn)確率的主要季節(jié),這與冬季逆溫影響有關(guān)。

    2.3 重要天氣過程的溫度預(yù)報(bào)檢驗(yàn)結(jié)果

    為檢驗(yàn)ECMWF細(xì)網(wǎng)格模式2m溫度預(yù)報(bào)未來24h產(chǎn)品在不同天氣過程中的預(yù)報(bào)能力,選取2011年7月—2014年6月影響天氣出現(xiàn)日,共304d有效樣本,依照不同天氣性質(zhì)特點(diǎn)分成降雨日、降雪日、東南大風(fēng)日、大霧日四類。由于雨雪春秋季多為雨雪轉(zhuǎn)換時(shí)期,不好區(qū)分降雨降雪狀態(tài),為方便研究,我們把4—10月08—08時(shí)出現(xiàn)0.1mm以上的降水定義為降雨日,共108d;11月—翌年3月08—08時(shí)出現(xiàn)0.1mm以上的降水定義為降雪日,共92d;全年08—08時(shí)出現(xiàn)極大風(fēng)速超過13m·s-1(5級(jí))的東南風(fēng)定義為東南風(fēng)日,共39d;全年08—08時(shí)出現(xiàn)能見度小于1km的霧,定義為大霧日,共65d,分別分析各類天氣發(fā)生時(shí)期,模式的溫度預(yù)報(bào)效果(表1)。

    表1 2011年7月—2014年6月ECMWF細(xì)網(wǎng)格2m溫度在不同天氣情況下最高、最低氣溫預(yù)報(bào)效果

    降雨日時(shí)最高、最低氣溫預(yù)報(bào)準(zhǔn)確率相對(duì)較高,分別為77.8%、78.7%,均超過了全年的平均準(zhǔn)確率,且穩(wěn)定性較好,平均絕對(duì)誤差都在1.5℃以下,從平均誤差看,模式預(yù)報(bào)普遍較實(shí)況偏高,最低氣溫預(yù)報(bào)表現(xiàn)更為明顯,正誤差達(dá)到了1℃,最高、最低氣溫預(yù)報(bào)的最大誤差分別為5.3℃和5℃,但出現(xiàn)大幅偏離實(shí)況的次數(shù)不多,108d中,最高、最低氣溫預(yù)報(bào)都只有3d偏離幅度超過4℃,且均為正誤差。

    降雪日最高、最低氣溫預(yù)報(bào)準(zhǔn)確率低于降雨日,分別為68.5%和65.2%,這由于降雪日時(shí)逆溫的破壞與維持都會(huì)對(duì)溫度預(yù)報(bào)有影響。最高氣溫預(yù)報(bào)的穩(wěn)定性略好于最低氣溫。最高、最低氣溫預(yù)報(bào)平均絕對(duì)誤差分別為1.5℃和1.7℃,平均誤差都較小,無法找出明顯的系統(tǒng)誤差,最大誤差分別為5.5℃和-5.4℃。在92d中,最高氣溫預(yù)報(bào)偏離幅度超過4℃的情況出現(xiàn)了4次,有2次正誤差,2次負(fù)誤差;最低氣溫預(yù)報(bào)偏離幅度超過4℃的情況出現(xiàn)了8次,有3次正誤差,5次負(fù)誤差。

    烏魯木齊市東南大風(fēng)出現(xiàn)時(shí),由于空氣干絕熱下沉運(yùn)動(dòng),會(huì)使溫度異常偏高,增加了對(duì)溫度預(yù)報(bào)的難度,ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品在東南風(fēng)日的溫度預(yù)報(bào)中最高氣溫預(yù)報(bào)較好,準(zhǔn)確率可以達(dá)到71.8%,最低氣溫預(yù)報(bào)較差,準(zhǔn)確率僅為59%。最高氣溫預(yù)報(bào)相對(duì)穩(wěn)定,平均絕對(duì)誤差為1.65℃,最大誤差為-4.6℃,偏離幅度超過4℃的情況僅出現(xiàn)了1次;最低氣溫預(yù)報(bào)穩(wěn)定性較差,平均絕對(duì)誤差為2.3℃,最大誤差到達(dá)-5.8℃,偏離幅度超過4℃的情況達(dá)到了7次,且均為負(fù)誤差,出現(xiàn)概率超過17%。東南風(fēng)對(duì)最低氣溫的影響遠(yuǎn)大于最高氣溫,東南風(fēng)的起風(fēng)時(shí)間,持續(xù)時(shí)間以及強(qiáng)度的不同會(huì)對(duì)模式的最低氣溫預(yù)報(bào)造成不同影響,東南大風(fēng)在夜間的持續(xù)時(shí)間越長,對(duì)低溫的影響越大,7次較大的低溫預(yù)報(bào)誤差中,有4次都出現(xiàn)在持續(xù)2d以上的強(qiáng)東南風(fēng)過程,2次出現(xiàn)在隆冬季節(jié)短暫的東南大風(fēng)過程。

    在各類天氣中,ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品受大霧影響最大,最高、最低氣溫預(yù)報(bào)準(zhǔn)確率為53.8%和52.3%,準(zhǔn)確率是幾類天氣過程中最低的,且預(yù)報(bào)的穩(wěn)定性均較差,平均絕對(duì)誤差都超過了2℃。最高氣溫預(yù)報(bào)的最大誤差達(dá)12℃,偏離幅度超過4℃的情況出現(xiàn)15次,出現(xiàn)概率超過了23%,其中6次為正誤差,9次為負(fù)誤差;最低氣溫預(yù)報(bào)的最大誤差為-5.8℃,偏離幅度超過4℃的情況出現(xiàn)10次,出現(xiàn)概率超過了15%,其中2次為正誤差,8次為負(fù)誤差??梢?,大霧天氣下多會(huì)使模式最低氣溫預(yù)報(bào)嚴(yán)重低于實(shí)況。逆溫是造成大霧天氣的重要條件,大霧時(shí)溫度預(yù)報(bào)準(zhǔn)確率低與模式無法很好處理逆溫時(shí)的溫度預(yù)報(bào)有很大關(guān)系。

    總體來看,模式對(duì)降雨、降雪時(shí)的溫度預(yù)報(bào)較好;大霧時(shí)最差,并且會(huì)使最低氣溫預(yù)報(bào)嚴(yán)重低于實(shí)況;對(duì)東南風(fēng)溫度預(yù)報(bào)時(shí)最高氣溫預(yù)報(bào)好于最低氣溫,東南風(fēng)的起風(fēng)時(shí)間,持續(xù)時(shí)間以及強(qiáng)度的不同會(huì)對(duì)模式的最低氣溫預(yù)報(bào)造成不同影響,東南大風(fēng)在夜間的持續(xù)時(shí)間越長,對(duì)低溫預(yù)報(bào)的影響越大。

    2.4 ECMWF細(xì)網(wǎng)格預(yù)報(bào)逆溫時(shí)的檢驗(yàn)結(jié)果

    將2011年7月—2014年6月ECMWF細(xì)網(wǎng)格上格點(diǎn)850hPa溫度與925hPa的溫度差值,作為逆溫預(yù)報(bào)值,將未來24h最強(qiáng)逆溫大于8℃作為一個(gè)逆溫預(yù)報(bào)日,共63d。檢驗(yàn)表明,最高、最低氣溫預(yù)報(bào)準(zhǔn)確率為31.74%和63.5%,最高溫度的預(yù)報(bào)明顯偏低。最高氣溫預(yù)報(bào)的最大誤差達(dá)12℃,偏離幅度超過4℃的情況出現(xiàn)21次,出現(xiàn)概率超過33%,其中13次為正誤差,8次為負(fù)誤差。最低氣溫預(yù)報(bào)的最大誤差為-5.6℃,偏離幅度超過4℃的情況出現(xiàn)8次,出現(xiàn)概率為13%,均為負(fù)誤差??梢?,當(dāng)ECMWF細(xì)網(wǎng)格預(yù)報(bào)明顯逆溫時(shí),最高溫度預(yù)報(bào)不可信,但最低溫度可作參考,并且會(huì)使最低氣溫預(yù)報(bào)低于實(shí)況,這同大霧天氣的溫度預(yù)報(bào)情況相似。

    2.5 ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品與國家氣象中心指導(dǎo)預(yù)報(bào)對(duì)比檢驗(yàn)

    為更好地衡量ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品在烏魯木齊市的預(yù)報(bào)能力,選取其2013年整年資料,同國家氣象中心指導(dǎo)預(yù)報(bào)溫度產(chǎn)品進(jìn)行對(duì)比檢驗(yàn),發(fā)現(xiàn)ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品整體在準(zhǔn)確率和穩(wěn)定性上均優(yōu)于國家氣象中心指導(dǎo)預(yù)報(bào),預(yù)報(bào)時(shí)效越長,優(yōu)勢(shì)越明顯,尤其在春季溫度預(yù)報(bào)中有明顯優(yōu)勢(shì)。從圖3a可以看出ECMWF細(xì)網(wǎng)格模式的未來24h預(yù)報(bào)產(chǎn)品在3—5月,最高氣溫預(yù)報(bào)準(zhǔn)確率高出國家氣象中心指導(dǎo)預(yù)報(bào)30%~40%;最低氣溫預(yù)報(bào)準(zhǔn)確率高出國家氣象中心預(yù)報(bào)準(zhǔn)確率20%~35%;從圖3b來看,最高氣溫預(yù)報(bào)較指導(dǎo)預(yù)報(bào)在3—6月、10月有30%~40%的正技巧,1月、9月、11月為負(fù)技巧,其中11月的負(fù)技巧達(dá)67%;對(duì)于最低氣溫預(yù)報(bào)而言,2—4月、10月有30%左右的正技巧,1月為64%的負(fù)技巧,其他月份預(yù)報(bào)水平相當(dāng)。圖3c可以看出,ECMWF細(xì)網(wǎng)格模式的未來48h最高氣溫預(yù)報(bào)準(zhǔn)確率僅在7月、11月低于國家氣象中心指導(dǎo)預(yù)報(bào),其余月份都高于指導(dǎo)預(yù)報(bào),最低氣溫預(yù)報(bào)準(zhǔn)確率在2月、5月、8月、9月低于國家氣象中心指導(dǎo)預(yù)報(bào),其余月份都高于指導(dǎo)預(yù)報(bào);圖3d可以看出,僅在11月ECMWF細(xì)網(wǎng)格模式的未來48h最高氣溫預(yù)報(bào)較指導(dǎo)預(yù)報(bào)有30%左右負(fù)技巧,其余時(shí)間都接近或?yàn)槊黠@正技巧。雖然ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品整體在準(zhǔn)確率和穩(wěn)定性上均優(yōu)于國家氣象中心指導(dǎo)預(yù)報(bào),但在秋冬季的一些時(shí)間也會(huì)劣于國家氣象中心指導(dǎo)預(yù)報(bào),兩種模式產(chǎn)品存在一定互補(bǔ)性。

    圖3 2013年ECMWF細(xì)網(wǎng)格2m溫度產(chǎn)品同國家氣象中心指導(dǎo)預(yù)報(bào)對(duì)比檢驗(yàn)結(jié)果

    3 結(jié)論

    (1)ECMWF細(xì)網(wǎng)格模式2m溫度預(yù)報(bào)產(chǎn)品對(duì)制作烏魯木齊未來24h及48h溫度預(yù)報(bào)具有很好的指導(dǎo)意義。隨著季節(jié)的變化,不同參考點(diǎn)對(duì)烏魯木齊市最高、最低氣溫預(yù)報(bào)能力也相應(yīng)變化。多數(shù)時(shí)間格點(diǎn)映射站點(diǎn)的方法要優(yōu)于雙線性插值法,經(jīng)過分析,可以把T左下的日最大值代表2月、3月、5月、6月、7月、8月和9月的最高氣溫預(yù)報(bào),T右上的日最大值代表1月、4月、10月、11月和12月的最高氣溫預(yù)報(bào),T右上的日最小值代表1月、3月、4月、5月、9月、10月、11月、12月的最低氣溫預(yù)報(bào),T差值的日最小值代表6月、7月的最低氣溫預(yù)報(bào)。T左上的日最小值代表代表2月、8月的最低氣溫預(yù)報(bào)。未來24h全年最高、最低氣溫預(yù)報(bào)準(zhǔn)確率可以達(dá)到74%~75%,未來48h最高、最低氣溫預(yù)報(bào)準(zhǔn)確率較24h偏低10%左右,其中春季降低幅度最大,為10%~ 28%,夏秋季較小,為1%~5%。

    (2)模式季節(jié)變化明顯,夏半年的預(yù)報(bào)準(zhǔn)確率高出冬半年10%~25%,冬季最高、最低氣溫預(yù)報(bào)偏離實(shí)況程度最大,同時(shí)也是降低全年溫度預(yù)報(bào)準(zhǔn)確率的主要季節(jié)。

    (3)模式對(duì)降雨、降雪時(shí)的溫度預(yù)報(bào)較好,大霧時(shí)最差,并且會(huì)使最低氣溫預(yù)報(bào)嚴(yán)重低于實(shí)況;對(duì)東南風(fēng)時(shí)最高氣溫預(yù)報(bào)好于最低氣溫,東南風(fēng)的起風(fēng)時(shí)間,持續(xù)時(shí)間以及強(qiáng)度的不同會(huì)對(duì)模式的最低氣溫預(yù)報(bào)造成不同影響,東南大風(fēng)在夜間的持續(xù)時(shí)間越長,對(duì)低溫的影響越大;預(yù)報(bào)逆溫時(shí),最高氣溫預(yù)報(bào)不可信,但最低氣溫可作參考,并且最低氣溫預(yù)報(bào)低于實(shí)況。

    (4)該產(chǎn)品整體在準(zhǔn)確率和穩(wěn)定性上均優(yōu)于國家氣象中心指導(dǎo)預(yù)報(bào),預(yù)報(bào)時(shí)效越長,優(yōu)勢(shì)越明顯,尤其在春季溫度預(yù)報(bào)中有明顯優(yōu)勢(shì)。但在秋冬季的一些時(shí)間也會(huì)劣于國家氣象中心指導(dǎo)預(yù)報(bào),兩種模式產(chǎn)品存在一定互補(bǔ)性。

    [1]賈麗紅,盧新玉,崔彩霞,等.新疆15地州首府城市最高、最低氣溫預(yù)報(bào)檢驗(yàn)[J].沙漠與綠洲氣象,2013,7(1):16-20.

    [2]肖明靜,隋明,范蘇丹,等.3種數(shù)值模式溫度預(yù)報(bào)產(chǎn)品在山東應(yīng)用和誤差分析與訂正[J].干旱氣象,2012,30(3):472-477.

    [3]康嵐,馮漢中,屠妮妮,等.Grapes模式預(yù)報(bào)西南地區(qū)夏季2m溫度檢驗(yàn)評(píng)估[J].高原山地氣象研究,2009,29(2):26-32.

    [4]歐徽寧,鄒哲馨,孫小龍.中央臺(tái)精細(xì)化溫度指導(dǎo)預(yù)報(bào)產(chǎn)品在賀州市的檢驗(yàn)[J].氣象研究與應(yīng)用,2012,33(2suppl1):150-152.

    [5]賈麗紅,張慧嵐,肖開提·多萊特,等.T639數(shù)值產(chǎn)品對(duì)影響新疆主要系統(tǒng)的預(yù)報(bào)檢驗(yàn)評(píng)估[J].氣象,2012,38(8):971-976.

    [6]佟華,姚明明,王雨,等.T213L31全球中期數(shù)值天氣預(yù)報(bào)系統(tǒng)2m溫度誤差源分析[J].氣象,2006,32(2):52-57.

    [7]張冰,魏建蘇,裴海瑛.2006年T213模式在江蘇的降水和溫度檢驗(yàn)評(píng)估[J].氣象科學(xué),2008,28(4):468-472.

    [8]李國翠,連志鸞,趙彥廠,等.石家莊溫度預(yù)報(bào)檢驗(yàn)及影響因子分析[J].氣象與環(huán)境學(xué)報(bào),2009,25(1):15-18.

    [9]吳愛敏,路亞奇,李祥科.利用ECMWF產(chǎn)品對(duì)慶陽極端氣溫釋用效果分析[J].干旱氣象,2009,27(3):288-292.

    [10]湯浩,賈麗紅.背景場(chǎng)和物理過程對(duì)WRF模式在新疆區(qū)域預(yù)報(bào)性能的影響[J].沙漠綠洲與氣象,2014,8(3):27-33.

    [11]陳曉紅,朱佳寧,周揚(yáng)帆.EC細(xì)網(wǎng)格數(shù)值預(yù)報(bào)產(chǎn)品在一次大霧預(yù)報(bào)中的釋用[J].大氣科學(xué)應(yīng)用與研究,2012,43(2):94-100.

    [12]陳海鳳,黃世芹,金建德,等.EC細(xì)網(wǎng)格數(shù)值預(yù)報(bào)在貴陽地區(qū)的釋用效果分析[J].貴州氣象,2014,38(4):22-25.

    [13]陳海鳳,黃世芹,金建德,等.基于ECMWF細(xì)網(wǎng)格模式產(chǎn)品的濕螺旋度在四川盆地強(qiáng)降水預(yù)報(bào)中的應(yīng)用檢驗(yàn)[J].氣象2014,40(9):1088-1096.

    [14]陳海鳳,黃世芹,金建德,等.T639、ECMWF細(xì)網(wǎng)格模式對(duì)2012年5—8月四川盆地降水預(yù)報(bào)的天氣學(xué)檢驗(yàn)[J].高原山地氣象,2013,38(1):22-25.

    [15]張俊蘭,李圓圓,張超.ECMWF細(xì)網(wǎng)格模式降水產(chǎn)品在北疆暴雪中的應(yīng)用檢驗(yàn)[J].沙漠與綠洲氣象,2013,7 (4):7-13.

    [16]萬瑜,曹興,竇新英,等.ECMWF細(xì)網(wǎng)格數(shù)值預(yù)報(bào)產(chǎn)品在烏魯木齊東南大風(fēng)中的釋用[J].沙漠與綠洲氣象,2014,8(1):32-38.

    Test of ECMWF Refined Net Numerical 2m Temperature Forecast Data in Temperature Forecasting in Urumqi

    ZHANG Chao1,LI Na1,JIA Jian2

    (1.Xinjiang meteorological Observatory,Urumqi 830002,China;2.Urumqi meteorological Bureau,Urumqi 830002,China)

    Using ECMWF refined net numerical 2m temperature forecast data from July 2011 to June 2014,grid mapping site and bilinear interpolation method are used to select the optimized reference points,finally apply the synoptic testing method to test its temperature forecasting ability indifferent seasons and different types of weather in Urumqi.The result shows that it is helpful to making nest 24 hours and 48hours temperature forecast,The next 24hours prediction accuracy of annual maximum and minimum temperature is 74%~75%,the next 48hours forecast accuracy is 10%,lower than the next 24h;forecast accuracy of temperature is better in rainy days and snowy days,worse in foggy days;the actual minimum temperature is significantly below the forecast temperature;the forecast accuracy of maximum temperature is better than the forecast accuracy when southeastly wind happens;when the temperature inversion is forecast,the maximum temperature forecast cannot be trusted,but the minimum temperature forecast is trusty.The accuracy and stability of ECMWF refined net numerical 2m temperature forecast data are better than guidance forecast data from the National Meteorological Center of CMA.The longer the forecast period is,the more obvious advantages there are.

    ECMWF refined net;temperature test;2m temperature;Urumqi

    P457.3

    A

    1002-0799(2015)05-0062-07

    張超,李娜,賈健.ECMWF細(xì)網(wǎng)格模式2m溫度產(chǎn)品在烏魯木齊市溫度預(yù)報(bào)中的檢驗(yàn)[J].沙漠與綠洲氣象,2015,9(5):62-68.

    10.3969/j.issn.1002-0799.2015.05.010

    2014-05-07;

    2014-07-21

    中國沙漠氣象科學(xué)研究基金(Sqj2014009)和中國氣象局預(yù)報(bào)員專項(xiàng)(CMAYBY2014-078)共同資助。

    張超(1987-),男,工程師,現(xiàn)從事短期天氣預(yù)報(bào)分析與預(yù)報(bào)服務(wù)工作。E-mail:64418696@qq.com

    猜你喜歡
    格點(diǎn)最低氣溫準(zhǔn)確率
    帶有超二次位勢(shì)無限格點(diǎn)上的基態(tài)行波解
    一種電離層TEC格點(diǎn)預(yù)測(cè)模型
    乳腺超聲檢查診斷乳腺腫瘤的特異度及準(zhǔn)確率分析
    健康之家(2021年19期)2021-05-23 11:17:39
    不同序列磁共振成像診斷脊柱損傷的臨床準(zhǔn)確率比較探討
    2015—2017 年寧夏各天氣預(yù)報(bào)參考產(chǎn)品質(zhì)量檢驗(yàn)分析
    北辰地區(qū)日最低氣溫特征及影響因素研究
    天津科技(2019年3期)2019-03-30 07:17:14
    67年來濟(jì)南最高和最低氣溫變化特征
    帶可加噪聲的非自治隨機(jī)Boussinesq格點(diǎn)方程的隨機(jī)吸引子
    高速公路車牌識(shí)別標(biāo)識(shí)站準(zhǔn)確率驗(yàn)證法
    格點(diǎn)和面積
    国产成人精品福利久久| 亚洲精品色激情综合| 国产淫片久久久久久久久| 丰满乱子伦码专区| 日韩大片免费观看网站| 国产极品粉嫩免费观看在线 | 亚洲伊人久久精品综合| 日韩中文字幕视频在线看片| 久久久久精品久久久久真实原创| 插逼视频在线观看| 日韩一区二区视频免费看| 久久精品夜色国产| 婷婷色综合www| 我要看日韩黄色一级片| 国产日韩一区二区三区精品不卡 | 午夜福利影视在线免费观看| 春色校园在线视频观看| 久久青草综合色| 伊人久久精品亚洲午夜| 伊人亚洲综合成人网| 亚洲av成人精品一二三区| 国产精品女同一区二区软件| 色视频www国产| 免费人成在线观看视频色| 高清午夜精品一区二区三区| 日韩中字成人| 在线观看三级黄色| 国产亚洲欧美精品永久| 一级片'在线观看视频| 日韩精品有码人妻一区| 永久免费av网站大全| 成人免费观看视频高清| 十八禁网站网址无遮挡 | 久久热精品热| 成人毛片60女人毛片免费| 男人舔奶头视频| 搡女人真爽免费视频火全软件| 久久精品久久精品一区二区三区| 国产欧美日韩一区二区三区在线 | 国产真实伦视频高清在线观看| 国产视频首页在线观看| 热re99久久精品国产66热6| 久久久欧美国产精品| 男的添女的下面高潮视频| 99久久中文字幕三级久久日本| 青青草视频在线视频观看| 亚洲精品成人av观看孕妇| 成人午夜精彩视频在线观看| 婷婷色麻豆天堂久久| 亚洲欧洲精品一区二区精品久久久 | www.色视频.com| 丝瓜视频免费看黄片| 国产中年淑女户外野战色| 亚洲精品成人av观看孕妇| 建设人人有责人人尽责人人享有的| 在线观看www视频免费| 三上悠亚av全集在线观看 | 观看美女的网站| 99热这里只有是精品在线观看| 国产精品国产av在线观看| 久久人妻熟女aⅴ| 亚洲精品日韩av片在线观看| freevideosex欧美| 免费av不卡在线播放| 十八禁高潮呻吟视频 | 日韩熟女老妇一区二区性免费视频| 国产免费又黄又爽又色| 精品一品国产午夜福利视频| 欧美激情极品国产一区二区三区 | 国产成人午夜福利电影在线观看| 日韩视频在线欧美| 老女人水多毛片| 国产精品久久久久久av不卡| 国产高清国产精品国产三级| 午夜免费鲁丝| 国产精品一区二区在线不卡| 黄色怎么调成土黄色| 我要看黄色一级片免费的| 黑人巨大精品欧美一区二区蜜桃 | 国产日韩欧美在线精品| 少妇熟女欧美另类| 国产色婷婷99| 久热这里只有精品99| 国产日韩欧美亚洲二区| 一个人免费看片子| 观看免费一级毛片| 日本91视频免费播放| 人妻 亚洲 视频| av福利片在线观看| 免费不卡的大黄色大毛片视频在线观看| 美女cb高潮喷水在线观看| 国产成人精品一,二区| 人妻夜夜爽99麻豆av| 一本—道久久a久久精品蜜桃钙片| 青春草亚洲视频在线观看| 国产成人a∨麻豆精品| 激情五月婷婷亚洲| 精品一区二区免费观看| 69精品国产乱码久久久| 亚洲国产精品专区欧美| 日本wwww免费看| 亚洲精品乱久久久久久| 免费久久久久久久精品成人欧美视频 | 久久6这里有精品| 国产深夜福利视频在线观看| 成人免费观看视频高清| 久久免费观看电影| 哪个播放器可以免费观看大片| 男女边摸边吃奶| 久久久久久久久久久免费av| 99国产精品免费福利视频| 大香蕉久久网| 欧美日韩一区二区视频在线观看视频在线| 人妻夜夜爽99麻豆av| 久久精品久久久久久噜噜老黄| 毛片一级片免费看久久久久| 老女人水多毛片| 欧美精品一区二区大全| 精品少妇久久久久久888优播| 在线亚洲精品国产二区图片欧美 | 日韩中字成人| 国产精品久久久久成人av| 精品一区二区三卡| 欧美日韩在线观看h| 日日摸夜夜添夜夜爱| 国产乱人偷精品视频| 免费观看在线日韩| 欧美少妇被猛烈插入视频| 妹子高潮喷水视频| 人体艺术视频欧美日本| 成人特级av手机在线观看| 午夜福利影视在线免费观看| 久久人人爽人人片av| 久久精品国产亚洲网站| 日日摸夜夜添夜夜爱| 久久国产亚洲av麻豆专区| 最后的刺客免费高清国语| 国产精品99久久99久久久不卡 | 婷婷色综合www| 老司机影院成人| 国产 精品1| 国产视频内射| 能在线免费看毛片的网站| 美女xxoo啪啪120秒动态图| 精品少妇内射三级| 亚洲欧美清纯卡通| 人妻系列 视频| 亚洲国产毛片av蜜桃av| av网站免费在线观看视频| 亚洲情色 制服丝袜| 在线观看www视频免费| 午夜影院在线不卡| 乱码一卡2卡4卡精品| 熟女人妻精品中文字幕| 美女福利国产在线| 丝袜脚勾引网站| 精品亚洲乱码少妇综合久久| 国产一区有黄有色的免费视频| 在线观看免费视频网站a站| av福利片在线观看| 亚洲av福利一区| 在线 av 中文字幕| 97超碰精品成人国产| 深夜a级毛片| 大片免费播放器 马上看| 男女无遮挡免费网站观看| 日本与韩国留学比较| 美女内射精品一级片tv| 国产精品不卡视频一区二区| 99久久中文字幕三级久久日本| 亚洲欧美一区二区三区国产| 日韩,欧美,国产一区二区三区| 亚洲av男天堂| 国产精品一二三区在线看| 国产精品.久久久| 日韩大片免费观看网站| 国产日韩欧美视频二区| 中文字幕人妻丝袜制服| 久久久久精品久久久久真实原创| 亚洲国产精品专区欧美| 国产一区二区在线观看日韩| 在线观看免费日韩欧美大片 | 在线精品无人区一区二区三| 日本黄色日本黄色录像| 国产永久视频网站| av在线播放精品| 又黄又爽又刺激的免费视频.| 一级,二级,三级黄色视频| 美女xxoo啪啪120秒动态图| 国产精品久久久久久精品电影小说| 亚洲av在线观看美女高潮| 在线观看人妻少妇| 亚洲性久久影院| 精品人妻熟女毛片av久久网站| 国产欧美亚洲国产| 中文天堂在线官网| 日本vs欧美在线观看视频 | 欧美 日韩 精品 国产| 国产精品国产av在线观看| 最新中文字幕久久久久| 亚洲精品第二区| 看十八女毛片水多多多| 久久影院123| 国产日韩欧美亚洲二区| av免费观看日本| 日本欧美国产在线视频| 三级经典国产精品| 免费大片18禁| 欧美97在线视频| 哪个播放器可以免费观看大片| 男人和女人高潮做爰伦理| 国产av一区二区精品久久| 日韩欧美 国产精品| 成人国产麻豆网| 国产视频首页在线观看| 9色porny在线观看| 久久久久久久精品精品| av卡一久久| 国产日韩欧美视频二区| 亚洲精品国产av蜜桃| 国产欧美日韩精品一区二区| 国产熟女欧美一区二区| 精品久久久久久电影网| 久久久久网色| 亚洲国产日韩一区二区| 丰满乱子伦码专区| 在线观看www视频免费| 成人综合一区亚洲| 女性生殖器流出的白浆| 欧美日韩综合久久久久久| 22中文网久久字幕| 七月丁香在线播放| 久久综合国产亚洲精品| 老司机亚洲免费影院| 老熟女久久久| 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产日韩一区二区| 青春草亚洲视频在线观看| av在线老鸭窝| 国内少妇人妻偷人精品xxx网站| 五月玫瑰六月丁香| 日日摸夜夜添夜夜爱| 51国产日韩欧美| 香蕉精品网在线| 91精品国产九色| 国产伦理片在线播放av一区| 91午夜精品亚洲一区二区三区| 在线观看三级黄色| 99热这里只有是精品在线观看| 人人妻人人澡人人爽人人夜夜| 精品国产一区二区久久| 一区二区三区免费毛片| 久久精品久久久久久久性| 热re99久久精品国产66热6| 少妇精品久久久久久久| 大陆偷拍与自拍| 亚洲av.av天堂| 免费观看无遮挡的男女| 在线亚洲精品国产二区图片欧美 | 精华霜和精华液先用哪个| 国产精品一区www在线观看| 老司机影院毛片| 久久综合国产亚洲精品| 久热这里只有精品99| 亚洲av男天堂| 亚洲精品一区蜜桃| 这个男人来自地球电影免费观看 | 九九爱精品视频在线观看| 老女人水多毛片| 九九在线视频观看精品| 一级毛片 在线播放| 王馨瑶露胸无遮挡在线观看| 十八禁高潮呻吟视频 | 午夜激情福利司机影院| 人体艺术视频欧美日本| 好男人视频免费观看在线| 国产精品99久久久久久久久| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产欧美日韩在线播放 | 亚洲精品乱码久久久v下载方式| 亚洲真实伦在线观看| 国产精品熟女久久久久浪| 欧美激情极品国产一区二区三区 | 一级a做视频免费观看| 十八禁高潮呻吟视频 | 欧美一级a爱片免费观看看| 亚洲精品日韩av片在线观看| 秋霞在线观看毛片| 亚洲激情五月婷婷啪啪| 在线 av 中文字幕| 寂寞人妻少妇视频99o| 精品卡一卡二卡四卡免费| 免费人成在线观看视频色| 欧美亚洲 丝袜 人妻 在线| 黄色怎么调成土黄色| 国产黄色视频一区二区在线观看| 在线观看www视频免费| 另类亚洲欧美激情| 亚洲国产精品一区二区三区在线| 亚洲自偷自拍三级| 女性被躁到高潮视频| 少妇精品久久久久久久| 欧美精品一区二区大全| 99九九线精品视频在线观看视频| 久久久久久久国产电影| 亚洲天堂av无毛| 免费在线观看成人毛片| 一级,二级,三级黄色视频| 日韩亚洲欧美综合| 免费观看av网站的网址| 国产免费一区二区三区四区乱码| av福利片在线观看| 一本—道久久a久久精品蜜桃钙片| 免费观看的影片在线观看| 日本91视频免费播放| 日韩成人伦理影院| 99热6这里只有精品| 色5月婷婷丁香| 国产 一区精品| 久久 成人 亚洲| 91成人精品电影| 亚洲av日韩在线播放| 久久人人爽人人爽人人片va| 国产乱来视频区| 看十八女毛片水多多多| 美女福利国产在线| 国产黄片美女视频| 精品国产国语对白av| 大香蕉久久网| 嫩草影院入口| 波野结衣二区三区在线| 一区二区三区免费毛片| 91精品一卡2卡3卡4卡| 在线观看免费日韩欧美大片 | 久久热精品热| 九九久久精品国产亚洲av麻豆| 日韩成人av中文字幕在线观看| 亚洲不卡免费看| 欧美最新免费一区二区三区| 中国三级夫妇交换| 国产 精品1| 亚洲av二区三区四区| 一级毛片电影观看| 亚洲丝袜综合中文字幕| 亚洲精品色激情综合| 熟妇人妻不卡中文字幕| 国产熟女欧美一区二区| 国产亚洲午夜精品一区二区久久| 亚洲欧美清纯卡通| 男人添女人高潮全过程视频| 在线观看免费日韩欧美大片 | 我要看黄色一级片免费的| av天堂中文字幕网| 午夜91福利影院| 亚洲经典国产精华液单| 午夜老司机福利剧场| av福利片在线观看| 亚洲欧美清纯卡通| 亚洲欧美精品自产自拍| 下体分泌物呈黄色| 国产黄色免费在线视频| 人妻少妇偷人精品九色| 欧美最新免费一区二区三区| 国国产精品蜜臀av免费| 国产精品伦人一区二区| 国产精品人妻久久久久久| 国产免费一级a男人的天堂| 国产成人91sexporn| 王馨瑶露胸无遮挡在线观看| 久久久精品免费免费高清| 国内揄拍国产精品人妻在线| 国产精品熟女久久久久浪| 国产又色又爽无遮挡免| 国产真实伦视频高清在线观看| 中文字幕精品免费在线观看视频 | 久久久a久久爽久久v久久| 女人精品久久久久毛片| videos熟女内射| 丰满人妻一区二区三区视频av| 一个人免费看片子| 免费久久久久久久精品成人欧美视频 | 亚洲av在线观看美女高潮| 尾随美女入室| 午夜av观看不卡| 91在线精品国自产拍蜜月| 久久精品国产鲁丝片午夜精品| 久久热精品热| 久久午夜福利片| 久久这里有精品视频免费| 成人国产麻豆网| 欧美精品亚洲一区二区| 大香蕉久久网| 国产av国产精品国产| 亚洲中文av在线| 国产av国产精品国产| 成人免费观看视频高清| 9色porny在线观看| 免费人成在线观看视频色| 日本wwww免费看| 99热6这里只有精品| 伊人久久精品亚洲午夜| 蜜臀久久99精品久久宅男| 噜噜噜噜噜久久久久久91| 蜜桃久久精品国产亚洲av| 精品人妻偷拍中文字幕| 99久久综合免费| 精品视频人人做人人爽| 五月伊人婷婷丁香| 人妻少妇偷人精品九色| av在线观看视频网站免费| 久久人人爽人人片av| 男女免费视频国产| 国产高清有码在线观看视频| 91成人精品电影| 亚洲精品久久久久久婷婷小说| av天堂久久9| 老熟女久久久| 国产免费福利视频在线观看| 在线观看美女被高潮喷水网站| 久久久国产欧美日韩av| 男女啪啪激烈高潮av片| 精品少妇内射三级| 国产又色又爽无遮挡免| 色视频在线一区二区三区| 国产精品99久久久久久久久| 中文字幕av电影在线播放| 久久国产精品大桥未久av | 国产有黄有色有爽视频| 亚洲国产精品成人久久小说| 国产精品福利在线免费观看| 丝瓜视频免费看黄片| 又大又黄又爽视频免费| 免费观看性生交大片5| 国产欧美另类精品又又久久亚洲欧美| 在线免费观看不下载黄p国产| 国产在线免费精品| 少妇的逼好多水| 中文字幕人妻熟人妻熟丝袜美| 亚洲真实伦在线观看| 久久久久网色| 国产高清不卡午夜福利| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品一区二区三区在线| 韩国av在线不卡| 哪个播放器可以免费观看大片| 老熟女久久久| 亚洲内射少妇av| 久久毛片免费看一区二区三区| 黄色怎么调成土黄色| 波野结衣二区三区在线| 国产精品久久久久久精品电影小说| 22中文网久久字幕| 日韩亚洲欧美综合| 国产有黄有色有爽视频| 啦啦啦视频在线资源免费观看| 在线播放无遮挡| 日日撸夜夜添| 国产精品久久久久久精品古装| 国产黄片视频在线免费观看| 老女人水多毛片| 美女主播在线视频| 午夜激情福利司机影院| 纵有疾风起免费观看全集完整版| 91久久精品国产一区二区成人| a级片在线免费高清观看视频| 2018国产大陆天天弄谢| 夜夜爽夜夜爽视频| 国产日韩一区二区三区精品不卡 | 女性生殖器流出的白浆| 在现免费观看毛片| 欧美精品一区二区免费开放| 精品亚洲成a人片在线观看| 少妇的逼好多水| 99国产精品免费福利视频| 亚洲国产最新在线播放| 亚洲精品自拍成人| 青春草国产在线视频| 午夜日本视频在线| 亚洲欧美一区二区三区国产| 亚洲av日韩在线播放| 丰满少妇做爰视频| 成年人免费黄色播放视频 | 国内精品宾馆在线| 国产精品一区二区三区四区免费观看| 18禁在线播放成人免费| 精品一品国产午夜福利视频| 午夜激情久久久久久久| 高清在线视频一区二区三区| 最黄视频免费看| 91精品国产九色| 国产 一区精品| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 国产爽快片一区二区三区| 日韩中字成人| 啦啦啦啦在线视频资源| 高清黄色对白视频在线免费看 | 人体艺术视频欧美日本| 久久人人爽人人片av| 日韩不卡一区二区三区视频在线| h视频一区二区三区| 国产黄片视频在线免费观看| 男女免费视频国产| 国产精品久久久久久久久免| 在线观看三级黄色| 久久精品国产亚洲av涩爱| 99久久精品国产国产毛片| 在线观看av片永久免费下载| 熟妇人妻不卡中文字幕| 夫妻午夜视频| 色94色欧美一区二区| 水蜜桃什么品种好| 久久久久精品久久久久真实原创| 亚洲丝袜综合中文字幕| 久久国产亚洲av麻豆专区| 天天操日日干夜夜撸| 18禁动态无遮挡网站| 亚州av有码| 麻豆成人av视频| 中文字幕精品免费在线观看视频 | 精品人妻偷拍中文字幕| 自线自在国产av| 啦啦啦在线观看免费高清www| 人妻 亚洲 视频| 国产精品女同一区二区软件| 内射极品少妇av片p| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 久久午夜综合久久蜜桃| 欧美高清成人免费视频www| 有码 亚洲区| 久久久a久久爽久久v久久| 搡女人真爽免费视频火全软件| 欧美成人精品欧美一级黄| 91aial.com中文字幕在线观看| 国产日韩欧美视频二区| 亚洲激情五月婷婷啪啪| 99热全是精品| 日本与韩国留学比较| 性色av一级| 国产一区二区在线观看av| 国产成人午夜福利电影在线观看| 国模一区二区三区四区视频| 国产精品偷伦视频观看了| 国产高清三级在线| 国产黄片视频在线免费观看| 亚洲精品aⅴ在线观看| 狠狠精品人妻久久久久久综合| 亚洲欧洲日产国产| 免费av不卡在线播放| 国产日韩欧美视频二区| 18禁裸乳无遮挡动漫免费视频| 久久久久久久久久成人| 精品一区二区三区视频在线| 久久精品久久久久久噜噜老黄| 成人午夜精彩视频在线观看| av国产久精品久网站免费入址| 汤姆久久久久久久影院中文字幕| 精品亚洲成国产av| 嫩草影院新地址| 边亲边吃奶的免费视频| 99热这里只有精品一区| 国产视频首页在线观看| 你懂的网址亚洲精品在线观看| 丰满饥渴人妻一区二区三| 爱豆传媒免费全集在线观看| 少妇被粗大的猛进出69影院 | 在线亚洲精品国产二区图片欧美 | 亚洲欧美精品专区久久| 欧美丝袜亚洲另类| 久久国产乱子免费精品| 亚洲成人一二三区av| 日本色播在线视频| 如何舔出高潮| 国产成人91sexporn| 精品久久久噜噜| 欧美少妇被猛烈插入视频| 高清视频免费观看一区二区| 久久国产亚洲av麻豆专区| 夜夜看夜夜爽夜夜摸| 熟女人妻精品中文字幕| 久久国产乱子免费精品| 国产欧美日韩精品一区二区| 美女xxoo啪啪120秒动态图| 国产在线免费精品| 亚洲欧美清纯卡通| 22中文网久久字幕| 水蜜桃什么品种好| 免费人妻精品一区二区三区视频| 亚洲av综合色区一区| 日本av免费视频播放| 久久久久久久大尺度免费视频| 国产欧美日韩综合在线一区二区 | 欧美高清成人免费视频www| 偷拍熟女少妇极品色| av有码第一页| 大香蕉97超碰在线| a级毛片在线看网站| 国产在线一区二区三区精| 欧美性感艳星| 亚洲精品aⅴ在线观看| 丰满饥渴人妻一区二区三| 亚洲欧洲国产日韩| 亚洲av成人精品一区久久| 国产成人精品婷婷| 免费看不卡的av| 国产黄片视频在线免费观看| 国产精品一区二区三区四区免费观看| 大片免费播放器 马上看| av天堂中文字幕网| 热99国产精品久久久久久7| 亚洲精品视频女| 麻豆精品久久久久久蜜桃|