周永水,齊大鵬,吳古會(huì)
(貴州省氣象臺(tái),貴州 貴陽(yáng) 550002)
?
貴陽(yáng)14時(shí)氣溫?cái)?shù)值預(yù)報(bào)檢驗(yàn)及擬合訂正分析
周永水,齊大鵬,吳古會(huì)
(貴州省氣象臺(tái),貴州貴陽(yáng)550002)
利用2012—2014年T639、EC細(xì)網(wǎng)格24 h、48 h內(nèi)2 m溫度和站點(diǎn)觀測(cè)氣溫,在對(duì)比分析T639和EC細(xì)網(wǎng)格模式優(yōu)缺點(diǎn)的基礎(chǔ)上,選擇預(yù)報(bào)準(zhǔn)確率更高的EC細(xì)網(wǎng)格14時(shí)2 m溫度和14時(shí)觀測(cè)氣溫求解最優(yōu)回歸擬合方程,考察在不同天氣轉(zhuǎn)換過(guò)程中模式和擬合方程對(duì)14時(shí)氣溫預(yù)報(bào)的優(yōu)缺點(diǎn)可知:24 h內(nèi)回歸檢驗(yàn)擬合方程比模式預(yù)報(bào)準(zhǔn)確率高4.92%,對(duì)雨天和陰/多云天氣比模式預(yù)報(bào)有顯著提升,預(yù)報(bào)實(shí)時(shí)檢驗(yàn)也比模式高4.25%,對(duì)提高最高氣溫預(yù)報(bào)具有一定的參考價(jià)值。
14時(shí)氣溫;擬合;天氣轉(zhuǎn)換;準(zhǔn)確率
氣溫預(yù)報(bào)一直是天氣預(yù)報(bào)的重要組成部分,近年來(lái)由于全球氣候變化的影響,嚴(yán)寒、高溫?zé)崂?、短時(shí)強(qiáng)降水等極端天氣現(xiàn)象頻發(fā)[1][2],民眾對(duì)于天氣預(yù)報(bào)的要求和關(guān)注度也越來(lái)越高。隨著社會(huì)的發(fā)展以及人們生活水平的不斷提高,天氣預(yù)報(bào)已深入社會(huì)行業(yè)的各個(gè)方面。氣溫預(yù)報(bào)既是干旱、霜凍、凍雨等災(zāi)害天氣預(yù)報(bào)的基礎(chǔ),也是農(nóng)林牧漁部門(mén)的作業(yè)和規(guī)劃的重要參考,同時(shí)對(duì)于群眾的生活、工作以及出行都有重要的指示意義。
氣溫預(yù)報(bào)最近幾十年得到了快速的發(fā)展,目前已有臺(tái)站開(kāi)展和嘗試氣溫預(yù)報(bào)的精細(xì)化研究以及對(duì)數(shù)值預(yù)報(bào)的氣溫檢驗(yàn)并建立本地預(yù)報(bào)模型[3-6],但是對(duì)單個(gè)站點(diǎn)的氣溫預(yù)報(bào),主要關(guān)注的是最高氣溫和最低氣溫的預(yù)報(bào)。在2012—2014年這3 a全國(guó)城鎮(zhèn)天氣預(yù)報(bào)檢驗(yàn)評(píng)估報(bào)告中貴州地區(qū)24 h最高氣溫準(zhǔn)確率分別為64.23%、66.98%、68.95%,最低氣溫準(zhǔn)確率分別為81.97%、81.86%、86.97%,因此,在整體改進(jìn)要素預(yù)報(bào)的同時(shí),對(duì)最高氣溫的改進(jìn)和訂正顯得尤為重要和緊迫。在氣溫預(yù)報(bào)中,一般是在客觀的數(shù)值預(yù)報(bào)的基礎(chǔ)上進(jìn)行主觀訂正后做出最高氣溫和最低氣溫的預(yù)報(bào),因此,客觀的數(shù)值預(yù)報(bào)是氣溫預(yù)報(bào)準(zhǔn)確率的基礎(chǔ)和保證,而直接對(duì)氣溫這一要素做出預(yù)報(bào)并在貴州大范圍使用的數(shù)值模式主要是T639和EC細(xì)網(wǎng)格的2 m溫度預(yù)報(bào),這兩個(gè)模式的空間分辨率分別為0.125°×0.125°、1°×1°,時(shí)間分辨率均為3 h,這使得它們的正點(diǎn)預(yù)報(bào)時(shí)間段內(nèi)均會(huì)出現(xiàn)14時(shí)這一時(shí)間節(jié)點(diǎn),通過(guò)對(duì)2001—2010年10 a貴陽(yáng)歷史氣溫?cái)?shù)據(jù)做相關(guān)分析發(fā)現(xiàn),14時(shí)的氣溫和當(dāng)天最高氣溫相關(guān)系數(shù)為0.668 7,且通過(guò)0.01的置信度檢驗(yàn),對(duì)當(dāng)天最高氣溫有良好的反映。由此,14時(shí)的氣溫預(yù)報(bào)可成為最高氣溫預(yù)報(bào)的重要參考因子。
本文主要利用2012—2014年T639、EC細(xì)網(wǎng)格24 h、48 h內(nèi)2 m溫度在14時(shí)預(yù)報(bào)值插值到各個(gè)觀測(cè)站點(diǎn)上,再結(jié)合每個(gè)站點(diǎn)14時(shí)歷史氣溫進(jìn)行對(duì)比分析,找出預(yù)報(bào)準(zhǔn)確率更高的預(yù)報(bào)模式,并在考慮天氣現(xiàn)象轉(zhuǎn)換的情況下對(duì)該模式進(jìn)行擬合優(yōu)化,以期提高最高氣溫的預(yù)報(bào)準(zhǔn)確率。由于貴州觀測(cè)站點(diǎn)較多,方程的擬合以單個(gè)站點(diǎn)為基礎(chǔ),本文主要介紹貴陽(yáng)站檢驗(yàn)及擬合情況。
本文使用的溫度觀測(cè)和預(yù)報(bào)資料包括:2012—2014年貴陽(yáng)觀測(cè)站逐日最高氣溫、最低氣溫、云量、天空狀況、海平面氣壓以及2012—2014年T639、EC細(xì)網(wǎng)格預(yù)報(bào)資料, 若實(shí)時(shí)資料或模式資料的14時(shí)數(shù)據(jù)缺失,則該日實(shí)況和模式資料對(duì)比無(wú)意義,舍去該日數(shù)據(jù),因此在表1的對(duì)比中EC細(xì)網(wǎng)格和T639樣本數(shù)量有所差異。本文使用的方法包括:常規(guī)統(tǒng)計(jì)分析、多項(xiàng)式逼近擬合法。經(jīng)統(tǒng)計(jì)分析得知:EC細(xì)網(wǎng)格前一日20時(shí)和當(dāng)日08時(shí)起報(bào)場(chǎng)對(duì)當(dāng)日14時(shí)預(yù)報(bào)(+18 h、+6 h)準(zhǔn)確率分別為66.03%、64.83%;前兩日20和前一日08時(shí)起報(bào)場(chǎng)對(duì)當(dāng)日14時(shí)預(yù)報(bào)(+42 h、+30 h)準(zhǔn)確率分別為64.29、62.75%;T639的該形式統(tǒng)計(jì)準(zhǔn)確率24 h為54.37%、48.21%,48 h準(zhǔn)確率為50.39%、53.49%,在預(yù)報(bào)時(shí)長(zhǎng)比08時(shí)多12 h的情況下,EC和T639的20時(shí)起報(bào)場(chǎng)準(zhǔn)確率24 h內(nèi)均高于08時(shí)起報(bào)場(chǎng),48 h內(nèi)EC20時(shí)起報(bào)準(zhǔn)確率仍高于08時(shí)起報(bào),因此主要關(guān)注20時(shí)起報(bào)場(chǎng)的檢驗(yàn)和訂正。
在分析過(guò)程中,將天空狀況分為晴天、陰/多云、雨天3種天氣類(lèi)型,以便分析和比較在各種天氣轉(zhuǎn)換中模式預(yù)報(bào)準(zhǔn)確率。
對(duì)2012—2014年兩種模式資料和實(shí)況對(duì)比可知(表1),EC細(xì)網(wǎng)格24 h和48 h對(duì)14時(shí)的氣溫預(yù)報(bào)準(zhǔn)確率為66.03%、64.29%,T639的這一數(shù)據(jù)分別是54.37%、50.39%,EC細(xì)網(wǎng)格準(zhǔn)確率比T639高出10%。再對(duì)天空狀況進(jìn)行分類(lèi),分為晴、陰/多云、雨3種天氣類(lèi)型,以比較不同天氣類(lèi)型轉(zhuǎn)換過(guò)程中的模式準(zhǔn)確率,由此可知:在24 h的氣溫預(yù)報(bào)中,EC細(xì)網(wǎng)格準(zhǔn)確率普遍高T639 10%以上,只有在雨轉(zhuǎn)晴、陰/多云轉(zhuǎn)雨和雨轉(zhuǎn)多云的天氣轉(zhuǎn)換中,T639預(yù)報(bào)準(zhǔn)確率略高于EC細(xì)網(wǎng)格或準(zhǔn)確率相當(dāng);這可能與貴陽(yáng)的陰雨天氣以及陰天天氣中一部分由靜止鋒影響造成,而對(duì)于靜止鋒在貴陽(yáng)附近的擺動(dòng)預(yù)報(bào)準(zhǔn)確率T639高于EC細(xì)網(wǎng)格,因此在這兩類(lèi)天氣變化過(guò)程中對(duì)貴陽(yáng)氣溫的預(yù)報(bào)準(zhǔn)確率T639高于EC細(xì)網(wǎng)格。同樣,在48 h的氣溫預(yù)報(bào)中,除雨轉(zhuǎn)陰/多云天氣T639預(yù)報(bào)準(zhǔn)確率高于EC細(xì)網(wǎng)格外,其余天氣轉(zhuǎn)換下EC細(xì)網(wǎng)格均比T639高10%以上。因此在方程擬合中只擬合EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)的數(shù)據(jù)。
表1 T639/EC細(xì)網(wǎng)格24 h、48 h不同天氣下14時(shí)氣溫預(yù)報(bào)準(zhǔn)確率(單位:%)對(duì)比
根據(jù)EC細(xì)網(wǎng)格20時(shí)起報(bào)場(chǎng)18 h和42 h預(yù)報(bào)2 m溫度數(shù)據(jù),對(duì)14時(shí)觀測(cè)氣溫進(jìn)行多項(xiàng)式形式逼近,求出24 h內(nèi)氣溫最優(yōu)化擬合方程為:
y=0.006 439×x2+0.842 5×x+0.158 5
其中y為擬合后14時(shí)氣溫預(yù)報(bào),x為前一日20時(shí)起報(bào)場(chǎng)預(yù)報(bào)的當(dāng)日14時(shí)2 m溫度, 其擬合曲線如圖1所示,其中橫坐標(biāo)為EC預(yù)報(bào)14時(shí)2 m溫度,縱坐標(biāo)為與預(yù)報(bào)對(duì)應(yīng)的14時(shí)實(shí)況氣溫,坐標(biāo)內(nèi)實(shí)心點(diǎn)由EC預(yù)報(bào)2 m溫度和14時(shí)氣溫共同決定,兩者比值越接近1表明預(yù)報(bào)和實(shí)況氣溫越吻合,擬合曲線代表著大多數(shù)實(shí)心點(diǎn)的分布趨勢(shì),曲線傾角越接近45°,表明實(shí)心點(diǎn)的橫坐標(biāo)和縱坐標(biāo)的比值越接近1,預(yù)報(bào)越接近實(shí)況,反之則表示預(yù)報(bào)越偏離實(shí)況。利用該方程對(duì)2012—2014年歷史資料進(jìn)行24 h和48 h回歸檢驗(yàn),可知擬合方程24 h預(yù)報(bào)準(zhǔn)確率為70.95%,高出EC預(yù)報(bào)場(chǎng)約5%,48 h預(yù)報(bào)準(zhǔn)確率擬合方程預(yù)報(bào)準(zhǔn)確率為67.11%,比EC預(yù)報(bào)場(chǎng)高2.82%。
將擬合方程和EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)準(zhǔn)確率具體到不同的天氣類(lèi)型,由表2和表3可知:24 h內(nèi)擬合方程對(duì)晴轉(zhuǎn)雨、晴轉(zhuǎn)多云、晴轉(zhuǎn)晴的擬合14時(shí)氣溫準(zhǔn)確率分別為70.83%、76.47%、76.15%,這與EC細(xì)網(wǎng)格14時(shí)氣溫預(yù)報(bào)準(zhǔn)確率大致相同,但是在其他幾種天氣轉(zhuǎn)換中,擬合方程預(yù)報(bào)準(zhǔn)確率均有不同程度的上升,特別是在陰/多云轉(zhuǎn)陰/多云和雨轉(zhuǎn)晴的天氣中,擬合方程準(zhǔn)確率比EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)高10%以上;剩下幾種天氣轉(zhuǎn)換擬合方程準(zhǔn)確率均高于EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng),從高到低分別為雨轉(zhuǎn)陰/多云高出9.43%、雨轉(zhuǎn)雨高出7.69%、陰/多云轉(zhuǎn)雨高出7.14%,陰/多云轉(zhuǎn)晴高出4.3%,因此,擬合方程24 h內(nèi)對(duì)晴天的準(zhǔn)確率無(wú)顯著提高,對(duì)雨天和陰天/多云的天氣預(yù)報(bào)準(zhǔn)確率有較大提高。
48 h擬合方程和EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)對(duì)比可知:擬合方程對(duì)晴天天氣轉(zhuǎn)換以及陰/多云轉(zhuǎn)晴的預(yù)報(bào)準(zhǔn)確率和EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)維持在同一水準(zhǔn),對(duì)雨天和陰/多云天氣轉(zhuǎn)換比較于EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)有更高的準(zhǔn)確率,對(duì)其余幾種天氣轉(zhuǎn)換擬合方程預(yù)報(bào)準(zhǔn)確率均有不同程度提高,按照比EC細(xì)網(wǎng)格預(yù)報(bào)場(chǎng)提高幅度從低到高排序分別是:雨轉(zhuǎn)雨高20%、雨轉(zhuǎn)晴高17.64%、雨轉(zhuǎn)陰/多云高10%、陰/多云轉(zhuǎn)陰/多云高7.19%、陰/多云轉(zhuǎn)雨高2.44%。該方程對(duì)雨天天氣轉(zhuǎn)換的準(zhǔn)確率有較大程度提高。
圖1 24 hEC細(xì)網(wǎng)格及實(shí)況14時(shí)氣溫(單位:℃)擬合曲線Fig.1 The fitted curve of 24-hour EC fine grid and real-time 14∶00 temperature (unit: ℃)
天氣類(lèi)型轉(zhuǎn)換樣本數(shù)量/個(gè)準(zhǔn)確率/%所有天氣形勢(shì)63066.03/70.95晴轉(zhuǎn)雨2470.83/70.83晴轉(zhuǎn)陰/多云8577.65/76.47晴轉(zhuǎn)晴13076.15/76.15陰/多云轉(zhuǎn)雨4252.38/59.52陰/多云轉(zhuǎn)陰/多云15961.64/71.70陰/多云轉(zhuǎn)晴9370.97/75.27雨轉(zhuǎn)晴1838.89/50.00雨轉(zhuǎn)陰/多云5347.17/56.60雨轉(zhuǎn)雨2661.54/69.23
表3 48 h不同天氣轉(zhuǎn)換下EC細(xì)網(wǎng)格及擬合方程回歸檢驗(yàn)準(zhǔn)確率(%)對(duì)比
在實(shí)時(shí)檢驗(yàn)過(guò)程中,檢驗(yàn)數(shù)據(jù)為2015年1—4月EC細(xì)網(wǎng)格資料、實(shí)況資料,因EC細(xì)網(wǎng)格資料缺失,總共可檢驗(yàn)天數(shù)約90 d,若細(xì)分到上文規(guī)定的9種天氣類(lèi)型轉(zhuǎn)換下則樣本過(guò)少,不具有較好的代表性,因此實(shí)時(shí)檢驗(yàn)把天氣轉(zhuǎn)換只分為3種天氣類(lèi)型,即按當(dāng)日天氣分為晴、陰/多云、雨,分析表4可知:24 h內(nèi)預(yù)報(bào)準(zhǔn)確率擬合方程比EC細(xì)網(wǎng)格高5.25%,高出的部分由陰/多云(高2.5%)、雨天(高9.09%)天氣貢獻(xiàn),晴天的預(yù)報(bào)準(zhǔn)確率擬合方程低于EC細(xì)網(wǎng)格;這與歷史擬合數(shù)據(jù)具有較好的一致性,表明擬合方程的對(duì)EC細(xì)網(wǎng)格預(yù)報(bào)擬合的延續(xù)性和有效性,是提高最高氣溫預(yù)報(bào)的一個(gè)有效的方法。
48 h內(nèi)預(yù)報(bào)準(zhǔn)確率擬合方程和EC細(xì)網(wǎng)格比較沒(méi)有明顯提升,且48 h細(xì)網(wǎng)格預(yù)報(bào)準(zhǔn)確率遠(yuǎn)遠(yuǎn)低于歷史預(yù)報(bào)準(zhǔn)確率,這可能與實(shí)時(shí)檢驗(yàn)所采集資料時(shí)段主要處于春季有關(guān),春季天氣變化快,預(yù)報(bào)流場(chǎng)調(diào)整大,48 h的響應(yīng)時(shí)間會(huì)造成較大誤差。
表4 擬合方程及EC細(xì)網(wǎng)格實(shí)時(shí)檢驗(yàn)準(zhǔn)確率(單位:%)對(duì)比
①T639和EC細(xì)網(wǎng)格資料對(duì)2 m氣溫的預(yù)報(bào)準(zhǔn)確率,在預(yù)報(bào)時(shí)效更長(zhǎng)的情況下,24 h和48 h的20時(shí)起報(bào)場(chǎng)均好于08時(shí)起報(bào)場(chǎng);EC細(xì)網(wǎng)格資料2 m氣溫預(yù)報(bào)24 h、48 h預(yù)報(bào)準(zhǔn)確率高于T639。
②擬合方程在陰/多云、雨天的天氣下對(duì)2 m溫度預(yù)報(bào)有更高的準(zhǔn)確率。
③24 h內(nèi)擬合方程比EC細(xì)網(wǎng)格實(shí)時(shí)預(yù)報(bào)準(zhǔn)確率高,這和歷史擬合數(shù)據(jù)具有較好的一致性,表明擬合方程的對(duì)EC細(xì)網(wǎng)格預(yù)報(bào)擬合的延續(xù)性和有效性,對(duì)提高最高氣溫預(yù)報(bào)有一定的參考價(jià)值。
[1] 羅亞麗.極端天氣和氣候事件的變化[J].氣候變化研究進(jìn)展,2012.8(2):90-98.
[2] 杜小玲,彭芳,武文輝.貴州凍雨頻發(fā)地帶分布特征及成因分析[J].氣象,2010,36(5):92-97.
[3] 邱學(xué)興,王東勇,朱紅芳.鄉(xiāng)鎮(zhèn)精細(xì)化最高最低氣溫預(yù)報(bào)方法研究[J].氣象與環(huán)境學(xué)報(bào),2013,29(3):92-96.
[4] 魏秀蘭,侯艷麗,孔凡中.菏澤市氣溫精細(xì)化預(yù)報(bào)[J].山東氣象,2004,24(96):45.
[5] 羅貴東,吳海周,肖鵬.基于MOS方法建立達(dá)州溫度預(yù)報(bào)模型[J].貴州氣象,2015,39(4):16-21.
[6] 符鳳平,吳哲紅,高如玉.精細(xì)化預(yù)報(bào)的自動(dòng)檢驗(yàn)與技巧評(píng)分系統(tǒng)[J].貴州氣象,2015,39(5):14-18.
Verification for Numerical Prediction and Correction with Curve Fitting on the Air Temperature of Guiyang at 14∶00
ZHOU Yongshui,QI Dapeng,WU Guhui
(Guizhou Meteorological Observatory,Guiyang 550002,China)
Based on the comparative analysis of the advantages on air temperature prediction of T639 and ECMWF's numerical forecast produce within 48 hours, the optimal regression equation was obtained by fitting the observed temperature with temperature prediction on 14∶00 at 2 meters of ECMWF, whose forecast produce's accuracy rate is better than T639 model. After analyzing the advantages of the equation on the temperature prediction on 14∶00 in different weather, we know that: the forecast accuracy rate of the equation is 4.92% higher than that of ECMWF in regression test. There is a significant improvement about the forecast accuracy rate of the equation on rainy and cloudy days. And equation's forecast accuracy rate of the verification for air temperature Real-Time forecasts is 4.25% higher than that of ECMWF. It has a certain value for improving the forecast of the highest temperature.
air temperature at 14∶00; fitting; changes in the weather; accuracy rate
1003-6598(2016)02-0045-04
2015-08-12
周永水(1982-),男,工程師,主要從事天氣預(yù)報(bào)工作,E-mail:joever2000@163.com。
貴州省氣象局2015年創(chuàng)新團(tuán)隊(duì)“貴州春夏主要強(qiáng)對(duì)流天氣預(yù)報(bào)方法研究”、貴州省氣象局業(yè)務(wù)技術(shù)攻關(guān)小組“貴州省日最高氣溫預(yù)報(bào)方法研究”共同資助。
P457.3
B