章大全 陳麗娟, 2

?

基于DERF2.0的月平均溫度概率訂正預(yù)報

章大全1陳麗娟1, 2

1中國氣象局國家氣候中心，北京100081;2南京信息工程大學(xué)氣象災(zāi)害預(yù)報預(yù)警與評估協(xié)同創(chuàng)新中心，南京210044

國家氣候中心第二代月動力延伸模式回算資料的分析表明，二代模式月平均溫度預(yù)報與觀測實況仍然存在較大偏差，模式預(yù)報有較大改進(jìn)空間。本文采用非參數(shù)百分位映射法對模式月平均溫度預(yù)報進(jìn)行概率訂正，該方法基于模式集合平均給出的確定性預(yù)報，結(jié)合模式回算資料各集合成員計算得到的模式概率密度分布，給出確定性預(yù)報在模式概率密度分布中的百分位值，并將百分位值投影到觀測資料的概率密度分布中，得到模式預(yù)報的概率訂正值。對訂正前后模式預(yù)報的檢驗評估顯示，該訂正方案不僅有效降低了模式預(yù)報與實況的均方根誤差（RMSE），對月平均溫度距平分布的預(yù)報技巧也有所改善，不同超前時間模式預(yù)報的預(yù)測技巧評分（PS）和距平相關(guān)系數(shù)（ACC）均有提升，同時模式預(yù)報誤差的大小對訂正效果無明顯影響。從分月的訂正預(yù)報結(jié)果來看，對夏季各月的溫度預(yù)測技巧的提升整體高于冬季各月。

月動力延伸模式 檢驗 概率 訂正

1 引言

隨著資料同化技術(shù)以及模式性能的提高，動力延伸預(yù)報模式已經(jīng)成為月尺度氣候預(yù)測的主要工具（Miyakoda et al., 1983, 1986；李維京，2012；鄭志海，2013）?；跀?shù)值動力氣候模式，歐洲中期天氣預(yù)報中心（ECMWF）、美國氣候預(yù)測中心（CPC）和東京氣候中心（TCC）等國際先進(jìn)業(yè)務(wù)中心都已經(jīng)開展了月尺度氣候預(yù)測業(yè)務(wù)（Frédéric, 2004; Saha et al., 2010, 2014）。國家氣候中心在“九五”期間和中國科學(xué)院大氣物理研究所合作研制了中等分辨率的全球大氣環(huán)流譜模式BCC_ AGCM1.0，并基于該模式建立了國家氣候中心第一代月動力延伸模式業(yè)務(wù)系統(tǒng)DERF1.0。該系統(tǒng)于2005年應(yīng)用于預(yù)測業(yè)務(wù)，成為我國短期氣候預(yù)測業(yè)務(wù)的主要工具之一，并在全球氣候預(yù)測信息交換中發(fā)揮了重要作用（陳麗娟和李維京，1999；丁一匯等，2002；李維京等，2005；覃志年等，2010；賈小龍等，2013）。2005年起國家氣候中心基于美國國家大氣研究中心（NCAR）的大氣環(huán)流模式CAM3.0，發(fā)展了第二代全球大氣環(huán)流譜模式，并且經(jīng)改進(jìn)優(yōu)化，于2011年3月完成BCC_AGCM2.2版本。BCC_AGCM2.2模式分辨率較之CAM2.0有所提高，水平分辨率從原來的T42提高到T106，垂直方向為26層（吳統(tǒng)文等，2013）。在BCC_AGCM2.2基礎(chǔ)上建立的第二代月動力業(yè)務(wù)延伸模式業(yè)務(wù)系統(tǒng)（DERF2.0）已于2014年正式投入業(yè)務(wù)運行。

DERF2.0預(yù)報產(chǎn)品采用滯后平均方案進(jìn)行逐日滾動預(yù)報，每日根據(jù)不同起報時間（間隔6 h）得到4個樣本成員的預(yù)報場，預(yù)報集合樣本數(shù)為20，積分時間為50 d。二代模式歷史回算資料的檢驗評估表明，其對全球及區(qū)域溫度、降水和環(huán)流的氣候態(tài)及年際變率等多個要素的預(yù)測能力總體要高于DERF1.0的預(yù)測技巧。何慧根等（2014）利用平均方差技巧等指標(biāo)對DERF2.0月平均溫度回算資料進(jìn)行評估發(fā)現(xiàn)，DERF2.0在西南地區(qū)的確定性預(yù)測效果較差，模式仍然有很大的改進(jìn)空間。此外DERF2.0在實時業(yè)務(wù)預(yù)測和歷史回算檢驗中分別采用了T639和NCEP兩套不同初始場，引入了額外的誤差，使得模式直接輸出的預(yù)報技巧提升有限。由于現(xiàn)有氣候系統(tǒng)模式的不完善，直接用模式結(jié)果開展短期氣候預(yù)測必然存在著一定的系統(tǒng)性誤差（陳紅等，2008；頡衛(wèi)華等，2010；郭準(zhǔn)等，2011）。在科研及業(yè)務(wù)應(yīng)用中，亟需發(fā)展和建立適應(yīng)我國氣候特點的模式誤差訂正系統(tǒng)，以提高我國短期氣候預(yù)測的技巧。

針對模式溫度預(yù)報與實況之間的偏差，如何充分挖掘模式預(yù)報信息，利用歷史資料對模式輸出結(jié)果進(jìn)行后期訂正，進(jìn)一步提高模式預(yù)報準(zhǔn)確率，是亟待解決的科研及業(yè)務(wù)問題。本文采用非參數(shù)百分位映射法對模式月平均溫度預(yù)報進(jìn)行概率訂正，該方法是近年來國外在模式誤差訂正中應(yīng)用較為普遍的一種方法。非參數(shù)百分位映射法的原理最早由Panofsky and Brier（1968）提出，Déqué（2007）首先將其應(yīng)用于區(qū)域氣候模式預(yù)報誤差訂正中，隨后該方法被廣泛應(yīng)用于模式誤差訂正、統(tǒng)計降尺度以及衛(wèi)星資料反演等（Dobler and Ahrens, 2008; Piani et al., 2010; Dosio and Paruolo, 2011; Douglas, 2013）。該方法通過將模式集合預(yù)報轉(zhuǎn)化為概率預(yù)報，并結(jié)合歷史觀測資料進(jìn)行概率訂正，在一定程度上提升了月平均溫度的預(yù)報技巧。Raisanen and Raty（2013）對10種模式溫度誤差訂正方案進(jìn)行對比，交叉檢驗結(jié)果顯示百分位映射法訂正效果最為明顯。由于將每個模式確定性預(yù)報都轉(zhuǎn)化為基于模式預(yù)報概率分布的百分位值，訂正預(yù)報結(jié)果包含了相對模式氣候態(tài)的偏移和極端性兩方面的信息。為了檢驗訂正方案對模式預(yù)報的改進(jìn)效果，給出了模式訂正前后回報檢驗（1983～2012年）結(jié)果的對比以及應(yīng)用實例。

2 資料和方法

本文所用資料主要包括：（1）國家氣候中心提供的全國160個站點1951～2012年月平均溫度資料；（2）DERF2.0的1983～2014年氣溫回算資料，模式氣候態(tài)采用1983～2012年共30年平均值。將模式格點數(shù)據(jù)插值到站點時采用雙線性插值算法。在分析模式相對實況的偏差時，模式資料選取起報時間為每月1日的月平均溫度預(yù)報（超前時間為0天）。為了系統(tǒng)地評估模式的預(yù)測能力，本文選用短期氣候業(yè)務(wù)常用預(yù)測評分方法：預(yù)測技巧評分（PS）、距平相關(guān)系數(shù)（ACC）、距平同號率（ACR）、均方根誤差（RMSE）等4種指標(biāo)（陳桂英和趙振國，1998）對預(yù)報結(jié)果進(jìn)行定量評估。其中PS評分是在預(yù)報與實況的距平符號一致率基礎(chǔ)上加上異常級加權(quán)得分構(gòu)成，包含對預(yù)報結(jié)果在氣候趨勢和氣候異常兩方面的檢驗評估。當(dāng)預(yù)報與實況完全一致時，預(yù)測評分PS為100。

3 DERF2.0月平均溫度預(yù)報場誤差分析

為定量評估DERF2.0模式回算資料與實況之間偏差，將模式每月1日起報的月平均溫度減去觀測得到模式預(yù)報誤差場，并分別給出模式誤差場的平均值、均方根誤差以及模式和觀測的偏度系數(shù)分布。其中偏度系數(shù)是描述氣候變量概率分布特征的重要統(tǒng)計量，表征分布形態(tài)與平均值偏離的程度，并作為分布不對稱的測度。偏度系數(shù)的計算公式（魏鳳英等，2008a，2008b）如下：

圖1分別給出了模式1月和7月溫度回報誤差場的30 a平均值和均方根誤差的空間分布。由圖1可以看出，除東北北部、內(nèi)蒙古北部、華南南部以及新疆中部部分地區(qū)外，模式1月平均溫度預(yù)報較之實況整體以偏低為主，其中東北中東部、內(nèi)蒙古東部、華北西部、新疆西部及北部局部、西南地區(qū)北部以及西藏等地平均溫度偏低4°C以上（圖1a）。模式7月溫度與實況偏差呈東西反向分布，105°E以西地區(qū)整體以偏冷為主，且大部地區(qū)偏冷幅度超過4°C。東部地區(qū)除東南沿海部分地區(qū)以外，模式溫度整體以偏高為主，其中偏暖的中心位于華北南部、黃淮西部和西南地區(qū)東部（圖1b）。從模式均方根誤差分布來看，1月的模式均方根誤差明顯高于7月，即模式冬季預(yù)報的偏差高于夏季。1月份我國東北大部、內(nèi)蒙古東部、新疆北部、西藏中部以及西南地區(qū)東部局部的模式均方根誤差可達(dá)2°C以上（圖1c）。7月模式平均溫度均方根誤差較小，但在華北南部、西北地區(qū)東部、黃淮、江淮和江南西部等地也超過1°C（圖1d）。

圖1 模式1月和7月平均溫度回報誤差的均值與均方根誤差分布（單位：°C）：（a）1月誤差平均值；（b）7月誤差平均值；（c）1月均方根誤差；（d） 7月均方根誤差

圖2分別為模式和觀測的1月和7月平均溫度偏度系數(shù)分布。從模式和觀測的對比結(jié)果來看，1月模式偏度分布與實況差異較大。實況中（圖2b）東北北部、內(nèi)蒙古中東部、華北、黃淮、江漢和西藏東部等地月平均溫度以左偏為主，溫度最大概然值低于月平均溫度，說明該區(qū)域出現(xiàn)負(fù)溫度距平概率較高。模式數(shù)據(jù)的偏度分布（圖2a）在東北南部、內(nèi)蒙古大部、華北、黃淮、西北地區(qū)東部、新疆北部、西南地區(qū)西部、西藏南部和東南沿海等地以右偏為主，與實況存在較大差異。7月的模式（圖2c）與實況（圖2d）的偏度系數(shù)空間分布比較一致，模式及實況均顯示東北地區(qū)南部、內(nèi)蒙古中部、華北、黃淮、西北地區(qū)東部、西南地區(qū)大部以及西藏東部等地月平均溫度概率分布為左偏，其余大部地區(qū)以右偏為主，說明模式對夏季溫度具有較好的模擬能力。

圖2 模式和觀測的1月與7月平均溫度偏度系數(shù)分布：（a）模式1月偏度；（b）觀測1月偏度；（c）模式7月偏度；（d）觀測7月偏度

4 基于非參數(shù)百分位映射的概率誤差訂正

針對模式預(yù)報與溫度實況之間的偏差，如何充分挖掘模式預(yù)報信息，同時利用歷史觀測資料對模式預(yù)報進(jìn)行有效訂正？傳統(tǒng)的誤差訂正方法僅僅包含模式單一預(yù)報值（均值）的訂正，完備的概率誤差訂正應(yīng)考慮模式預(yù)報概率密度分布與觀測實況概率密度分布之間的偏差，并對之進(jìn)行訂正。本文采用基于非參數(shù)百分位映射的概率誤差訂正方法對模式預(yù)報進(jìn)行訂正，該方法的原理（Panofsky et al., 1968; Glahn and Lowry, 1972）是基于模式集合平均給出的確定性預(yù)報，結(jié)合模式回算數(shù)據(jù)集合成員計算得到的模式概率密度分布，給出確定性預(yù)報在模式概率密度分布中的百分位值，并將百分位值投影到觀測資料的概率密度分布中，得到模式確定性預(yù)報的概率訂正值。該方法的優(yōu)點在于將每個模式確定性預(yù)報都轉(zhuǎn)化為基于模式預(yù)報分布的百分位值，因此預(yù)報結(jié)果包含了相對模式氣候態(tài)的偏移（距平）和極端性（概率）兩方面的信息。由于訂正預(yù)報僅參考了模式預(yù)報的百分位信息，因此預(yù)報結(jié)果有效避免了模式較之實況的系統(tǒng)性偏差。

圖3為基于非參數(shù)百分位映射的概率誤差訂正方法示意圖。下標(biāo)o, c, f分別代表觀測數(shù)據(jù)、模式回算數(shù)據(jù)和實時預(yù)報數(shù)據(jù)。圖3中為利用 模式回算數(shù)據(jù)得到的變量的累計概率分布，對于變量在時刻的預(yù)報值，對應(yīng)的累積概 率為。在觀測數(shù)據(jù)累積概率分布中找到變量的對應(yīng)值，并作為模式訂正預(yù)報值。由于月平均溫度概率密度分布近似呈正態(tài)分布，在訂正方案中將模式及觀測資料的概率密度分布（PDF）轉(zhuǎn)換為累積概率分布（CDF），使得分布函數(shù)形式具有單調(diào)性，也保證了映射結(jié)果的唯一性。

圖3 基于非參數(shù)百分位映射的概率誤差訂正示意圖

模式概率誤差訂正方案的操作流程如圖4所示。分別基于模式回報集合成員和歷史觀測數(shù)據(jù)計算模式回報和觀測溫度的概率密度分布，然后利用三次樣條函數(shù)插值得到模式回報和觀測溫度的累計概率密度分布。在回報檢驗和實時預(yù)測中，首先得到模式集合平均的確定性預(yù)報，結(jié)合模式回算數(shù)據(jù)的累計概率分布得到模式的百分位值預(yù)報，在當(dāng)月觀測資料的累計概率分布中找出該百分位值對應(yīng)的溫度值，并作為模式概率訂正的預(yù)報值。

圖4 模式溫度概率訂正流程

5 回報檢驗

利用模式1983-～2012年回算資料及國家氣候中心160站月平均溫度資料，對該概率訂正方案進(jìn)行檢驗評估，分別給出模式直接輸出和訂正后的30年平均的回報PS評分、ACC和ACR。檢驗評估結(jié)果顯示（圖5），模式對月平均溫度的預(yù)報性能在各月存在較大差異，冬季溫度預(yù)報技巧明顯高于夏季。當(dāng)超前時間為10天和5天時，模式1月平均溫度預(yù)報與實況的距平相關(guān)系數(shù)可以達(dá)到0.3左右，而6、7月份低于0.1。隨著超前時間的減少，模式對月平均溫度的預(yù)測技巧整體呈上升趨勢，其中模式預(yù)報的ACC上升最為明顯。21日起報的各月模式預(yù)測PS評分、ACC、ACR的全年平均值分別為75.4、0.15和61.1%，1日起報的模式預(yù)測評分分別達(dá)到78.3、0.32和66.5%。

圖5 DERF2.0溫度概率訂正回報檢驗（1983～2012年）

檢驗結(jié)果還顯示，訂正方案對模式不同起報時間的預(yù)測效果均有所改善，且提升效果與起報時間無明顯相關(guān)。從逐月的結(jié)果來看，方案對夏季溫度預(yù)測技巧的提升要高于冬季，近30年6到8月預(yù)測PS評分和ACR平均提升了3.2和4.8%。經(jīng)過誤差訂正后，模式預(yù)報的均方根誤差大大降低，其中21日、26日和1日起報的模式預(yù)報均方根誤差分別由3.18降至1.07，由3.15降至1.05和由3.10降至1.01。

為進(jìn)一步評估訂正方案對在不同月平均溫度異常及模式偏差背景下對模式預(yù)報技巧的改進(jìn)效果，將訂正后的模式預(yù)報PS評分減去未經(jīng)訂正的預(yù)報評分，得到訂正方案的PS評分提高值，并分別給出PS評分提高值與全國平均溫度距平（圖6）及模式均方根誤差（圖7）的關(guān)系。由圖6可知，訂正方案對占總樣本量77%的模式預(yù)報技巧有改善，且不同起報時間模式預(yù)報的訂正效果差別不大，顯示了該訂正方案的穩(wěn)定性。從溫度距平與PS評分提高值的關(guān)系來看，當(dāng)全國平均溫度距平接近0°C時，訂正效果較好。預(yù)測技巧的提升存在溫度不對稱性，當(dāng)月平均溫度以負(fù)距平為主時，預(yù)測技巧提升更為明顯。當(dāng)全國平均溫度距平為正時，訂正方案仍然以改善預(yù)報技巧為主，但部分時次的模式預(yù)報評分有所下降。模式預(yù)報的均方根誤差與預(yù)報技巧提升的關(guān)系表明（圖7），模式誤差的大小對預(yù)測技巧的提升影響不大，在模式預(yù)報均方根誤差達(dá)到4°C時，該訂正方案仍然能夠?qū)δＪ筋A(yù)報進(jìn)行有效的訂正。

圖6 模式預(yù)報PS評分提高值與全國平均溫度距平（單位：°C）的關(guān)系：（a）模式超前10天；（b）模式超前5天；（c）模式超前0天

圖7 模式預(yù)報PS評分提高值與模式預(yù)報均方根誤差（單位：°C）的關(guān)系：（a）模式超前10天；（b)模式超前5天；（c）模式超前0天

6 訂正方案實例

2008年1月，我國大部地區(qū)氣溫顯著偏低，南方遭受了歷史罕見的低溫、雨雪、冰凍災(zāi)害（王凌等，2008；高輝等，2008；張慶云等，2008；譚桂容等，2010）。2013年7月我國南方多地出現(xiàn)高溫伏旱，江南、江淮、江漢以及重慶等地出現(xiàn)持續(xù)高溫天氣，全國100個氣象觀測站發(fā)生極端高溫事件（林玉成等，2013）。下面分別以這兩個月份為例，從模式冬季和夏季預(yù)報兩個角度對訂正流程及效果進(jìn)行說明。

監(jiān)測實況顯示（圖8a），2008年1月我國除西南及東部部分地區(qū)外，全國溫度顯著偏低，其中西北、華中、華南、西南部分地區(qū)月平均溫度偏低2°C以上。模式1月1日起報的結(jié)果較好的反映了全國大范圍氣溫偏低的氣候特征（圖8b），模式預(yù)報溫度距平的空間分布與實況基本一致。預(yù)報偏差主要出現(xiàn)在東北及華北地區(qū)，實況東北南部及華北氣溫以偏低為主，而模式預(yù)報東北南部、華北東部、黃淮等地溫度為正距平。此外模式對低溫異常量級估計不足，模式預(yù)報我國東部及南方大部分地區(qū)溫度距平在?0.5°C以內(nèi)，而實況在華中、華南及西南部分地區(qū)氣溫負(fù)異常超過了2°C。訂正后的模式預(yù)測（圖8c）在東北、華北的溫度正距平區(qū)范圍縮小，華北及黃淮地區(qū)氣溫為負(fù)距平，與實況基本一致。另外從模式預(yù)報的訂正量（圖8d）可知，訂正方案對東部及南方大范圍地區(qū)進(jìn)行了負(fù)距平訂正，氣溫負(fù)異常量級更加接近實況，模式預(yù)報的PS評分由84.1上升至92.4，RMSE則由3.08°C降低至1.20°C（表1）。

圖8 2008年1月平均溫度距平的訂正效果實例（模式起報時間1月1日，單位：°C）：（a）溫度距平實況；（b）模式預(yù)測；（c）模式概率訂正預(yù)測；（d）模式訂正量

表1 2008年1月DERF2.0模式溫度距平預(yù)報訂正前、后評分

2013年7月全國大部分地區(qū)氣溫較常年同期偏高，其中黃淮南部、江淮、江漢、江南大部、西南東部等地氣溫偏高1～2°C，部分地區(qū)偏高2°C以上（圖9a）。模式7月1日起報的結(jié)果（圖9b）與實況之間存在較大偏差。模式預(yù)報江南中西部、江漢、西南地區(qū)北部、華南中部等地氣溫偏低，與實況出現(xiàn)趨勢性偏差。此外模式預(yù)報東北、華北等地氣溫整體偏高，其中東北西部、內(nèi)蒙古中東部、華北北部等地溫度偏高1°C以上，而實況顯示該區(qū)域大部分地區(qū)氣溫正常偏高，東北北部和華北西部部分地區(qū)溫度為負(fù)距平。概率訂正結(jié)果（圖9c）減少了模式預(yù)報南方低溫區(qū)的范圍和強(qiáng)度，同時降低了東北和華北溫度正距平異常的量級，使之更加接近實況。

圖9 2013年7月平均溫度距平的訂正效果實例（模式起報時間7月1日，單位：°C）：（a）溫度距平實況；（b）模式預(yù)測；（c）模式概率訂正預(yù)測；（d）模式訂正量

從模式預(yù)報的訂正量可以看出，該訂正方案對我國江南、華南、江漢、西南地區(qū)北部等模式預(yù)報溫度偏低的區(qū)域進(jìn)行了正距平訂正，對東北、華北、內(nèi)蒙古等預(yù)報溫度偏高的區(qū)域做了負(fù)距平訂正，訂正方案基本把握住了模式預(yù)報與實況之間的趨勢性偏差（圖9d）。從預(yù)報評分來看（表2），模式預(yù)報的PS評分、ACC以及ACR分別上升了5.2、0.09和11.9%，而模式預(yù)報的均方根誤差（RMSE）由2.81降至1.19。

表2 2013年7月DERF2.0模式溫度距平預(yù)報訂正前、后評分

7 小結(jié)與討論

隨著國家氣候中心第二代月動力延伸模式投入業(yè)務(wù)運行，針對模式數(shù)據(jù)的回報檢驗和解釋應(yīng)用工作急需開展。本文針對模式月平均溫度預(yù)報中存在的偏差，利用模式回算數(shù)據(jù)和歷史觀測資料的概率分布進(jìn)行誤差訂正，并對訂正效果進(jìn)行檢驗評估，取得一些有價值的結(jié)果。

本文采用非參數(shù)百分位映射法對模式月平均溫度預(yù)報進(jìn)行概率訂正，該方法利用歷史回算資料給出的模式概率密度分布信息，結(jié)合歷史觀測資料對模式預(yù)報進(jìn)行訂正，剔除了模式較之實況的系統(tǒng)性偏差對預(yù)報結(jié)果的影響，改進(jìn)了模式預(yù)報技巧。對訂正前后模式預(yù)報的檢驗評估顯示，該訂正方案有效的降低了模式預(yù)報與實況之間的均方根誤差，對月平均溫度距平空間分布的預(yù)報效果也有所改善。不同超前時間模式預(yù)報的PS和ACC評分均有提升，且預(yù)報技巧的改善與起報時間無顯著相關(guān)。從分月的結(jié)果來看，模式夏季溫度預(yù)測技巧的提升整體要高于冬季。對訂正效果與溫度距平和模式誤差關(guān)系的分析表明，訂正方案對模式預(yù)測技巧的提升存在溫度不對稱性，當(dāng)月平均溫度以負(fù)距平為主時，預(yù)測技巧的提升更為明顯，同時該方案具有誤差穩(wěn)定性，即預(yù)報誤差的大小對模式預(yù)報訂正的效果無明顯影響。

在實際業(yè)務(wù)應(yīng)用中，不僅需要參考模式訂正預(yù)報，模式預(yù)報訂正量的分布也值得關(guān)注。模式預(yù)報訂正量的物理意義在于模式能夠完美模擬出觀測資料概率密度分布情景下的預(yù)報結(jié)果與模式直接輸出之間的差值，因此模式預(yù)報訂正量（正負(fù)距平）可作為對模式誤差趨勢的判斷指標(biāo)。

另外，DERF2.0的歷史回算資料僅從1983年開始，每次預(yù)報的集合成員數(shù)只有20個，觀測資料的長度也僅有60年，在計算模式月平均溫度概率密度分布時由于樣本量不足，影響了對概率密度分布估計的準(zhǔn)確性，也在一定程度上影響了模式預(yù)報訂正的效果。隨著樣本量的增加，該缺陷可得到部分改善。

參考文獻(xiàn)：（References）

陳桂英, 趙振國. 1998. 短期氣候預(yù)測評估方法和業(yè)務(wù)初估[J]. 應(yīng)用氣象學(xué)報, 9 (2): 178–185. Chen Guiying, Zhao Zhenguo. 1998. Assessment methods of short range climate prediction and their operational application [J]. J. Appl. Meteor. Sci. (in Chinese), 9 (2): 178–185.

陳紅, 郎咸梅, 周廣慶, 等. 2008. 2008年1月中國氣候異常的動力學(xué)預(yù)測及效果檢驗 [J]. 氣候與環(huán)境研究, 13 (4): 531–538. Chen Hong, Lang Xianmei, Zhou Guangqing, et al. 2008. Prediction of climate anomalies over China in January 2008 with IAP climate models and its verification [J]. Climatic Environ. Res. (in Chinese), 13 (4): 531–538, doi:10.3878/j.issn.1006-9585.2008.04.17.

陳麗娟, 李維京. 1999. 月動力延伸預(yù)報產(chǎn)品的評估和解釋應(yīng)用 [J]. 應(yīng)用氣象學(xué)報, 10 (4): 486–490. Chen Lijuan, Li Weijing. 1999. The score skill and interpretation of monthly dynamic extended range forecast [J]. J. Appl. Meteor. Sci. (in Chinese), 10 (4): 486–490, doi:10.3969/j.issn. 1001-7313.1999.04.013.

Déqué M. 2007. Frequency of precipitation and temperature extremes over France in an anthropogenic scenario: Model results and statistical correction according to observed values [J]. Global and Planetary Change, 57 (1–2): 16–26, doi:10.1016/j.gloplacha.2006.11.030.

丁一匯, 劉一鳴, 宋永加, 等. 2002. 我國短期氣候動力預(yù)測模式系統(tǒng)的研究及試驗 [J]. 氣候與環(huán)境研究, 7 (2): 236–246. Ding Yihui, Liu Yiming, Song Yongjia, et al. 2002. Research and experiments of the dynamical model system for short-term climate prediction [J]. Climatic Environ. Res. (in Chinese), 7 (2): 236–246, doi:10.3969/j.issn.1006-9585. 2002.02.011.

Dobler A, Ahrens B. 2008. Precipitation by a regional climate model and bias correction in Europe and South Asia [J]. Meteor. Zeitschrift, 17 (4): 499–509, doi:10.1127/0941-2948/2008/0306.

Dosio A, Paruolo P. 2011. Bias correction of the ENSEMBLES high- resolution climate change projections for use by impact models: Evaluation on the present climate [J]. J. Geophys. Res., 116 (D16): D16106, doi:10.1029/2011JD015934.

Douglas M. 2013. Bias correction, quantile mapping, and downscaling: Revisiting the inflation issue [J]. J. Climate, 26 (6): 2137–2143, doi:10.1175/JCLI-D-12-00821.1.

Frédéric V. 2004. Monthly forecasting at ECMWF [J]. Mon. Wea. Rev., 132 (12): 2761–2779, doi:10.1175/MWR2826.1.

高輝, 陳麗娟, 賈小龍, 等. 2008. 2008年1月我國大范圍低溫雨雪冰凍災(zāi)害分析Ⅱ. 成因分析 [J]. 氣象, 34 (4): 101–106. Gao Hui, Chen Lijuan, Jia Xiaolong, et al. 2008. Analysis of the severe cold surge, ice-snow and frozen disasters in South China during January 2008: Ⅱ. Possible climate causes [J]. Meteor. Mon. (in Chinese), 34 (4): 101–106,doi:10.7519/j.issn.1000-0526.2008.4.013.

Glahn H R, Lowry D A. 1972. The use of model output statistics (MOS) in objective weather forecasting [J]. J. Appl. Meteor., 11 (8): 1203–1211, doi:10.1175/1520-0450(1972)011<1203:TUOMOS>2.0.CO;2.

郭準(zhǔn), 吳春強(qiáng), 周天軍, 等. 2011. LASG/IAP和BCC大氣環(huán)流模式模擬的云輻射強(qiáng)迫之比較[J]. 大氣科學(xué), 35 (4): 739–752. Guo Zhun, Wu Chunqiang, Zhou Tianjun, et al. 2011. A comparison of cloud radiative forcings simulated by LASG/IAP and BCC atmospheric general circulation models [J]. Chinese J. Atmos. Sci. (in Chinese), 35 (4): 739–752, doi:10.3878/j.issn.1006-9895.2011.04.13.

何慧根, 李巧萍, 吳統(tǒng)文, 等. 2014. 月動力延伸預(yù)測模式業(yè)務(wù)系統(tǒng)DERF2.0對中國氣溫和降水的預(yù)測性能評估[J]. 大氣科學(xué), 38 (5): 950–964. He Huigen, Li Qiaoping, Wu Tongwen, et al. 2014. Temperature and precipitation evaluation of monthly dynamic extended range forecast operational system DERF2.0 in China [J]. Chinese J. Atmos. Sci. (in Chinese), 38 (5): 950–964, doi:10.3878/j.issn.1006-9895. 1401.13166.

賈小龍, 陳麗娟, 高輝, 等. 2013. 我國短期氣候預(yù)測技術(shù)進(jìn)展 [J]. 應(yīng)用氣象學(xué)報, 24 (6): 641–655. Jia Xiaolong, Chen Lijuan, Gao Hui, et al. 2013. Advances of the short-range climate prediction in China [J]. J. Appl. Meteor. Sci. (in Chinese), 24 (6): 641–655, doi:10.3969/j.issn. 1001-7313.2013.06.001.

頡衛(wèi)華, 吳統(tǒng)文. 2010. 全球大氣環(huán)流模式BCC_AGCM2.0.1對1998年夏季江淮流域降水過程的回報試驗研究 [J]. 大氣科學(xué), 34 (5): 962–978. Jie Weihua, Wu Tongwen. 2010. Hindcast for 1998 summer heavy precipitation in the Yangtze and Huaihe River valley using AGCM BCC_AGCM2.0.1 model [J]. Chinese J. Atmos. Sci. (in Chinese), 34 (5): 962–978, doi:10.3878/j.issn.1006-9895.2010.05.11.

李維京. 2012. 現(xiàn)代氣候業(yè)務(wù) [M]. 北京: 氣象出版社, 202–315. Li Weijing. 2012. Modern Climate Operation [M]. Beijing: China Meteorological Press, 202–315.

李維京, 張培群, 李清泉, 等. 2005. 動力氣候模式預(yù)測系統(tǒng)業(yè)務(wù)化及其應(yīng)用 [J]. 應(yīng)用氣象學(xué)報, 16（增刊）：1–11. Li Weijing, Zhang Peiqun, Li Qingquan, et al. 2005. Research and operational application of dynamical climate model prediction system [J]. J. Appl. Meteor. Sci. (in Chinese), 16 (Suppl.): 1–11, doi:10.3969/j.issn.1001-7313.2005.z1.001.

林玉成, 徐珺, 張芳華. 2013. 2013年7月大氣環(huán)流和天氣分析 [J]. 氣象, 39 (10): 1379–1384. Lin Yucheng, Xu Jun, Zhang Fanghua. 2013. Analysis of the July 2013 atmospheric circulation and weather [J]. Meteor. Mon. (in Chinese), 39 (10): 1379–1384, doi:10.7519/j.issn.1000-0526. 2013.10.018.

Miyakoda K, Gordon T, Caverly R, et al. 1983. Simulation of a blocking event in January 1977 [J]. Mon. Wea. Rev., 111 (4): 846–869, doi:10.1175/1520-0493(1983)111<0846:SOABEI>2.0.CO;2.

Miyakoda K, Sirutis J, Ploshay J. 1986. One-month forecast experiments—Without anomaly boundary forcings [J]. Mon. Wea. Rev., 114 (12): 2363– 2401, doi:10.1175/1520-0493(1986)114<2363:OMFEAB>2.0.CO;2.

Panofsky H A, Brier G W. 1968. Some Applications of Statistics to Meteorology [M]. University Park, Pa: The Pennsylvania State University Press, 224–225.

Piani C, Haerter J O, Coppola E. 2010. Statistical bias correction for daily precipitation in regional climate models over Europe [J]. Theor. Appl. Climatol., 99 (1–2): 187–192, doi:10.1007/s00704-009-0134-9.

覃志年, 陳麗娟, 唐紅玉, 等. 2010. 月尺度動力模式產(chǎn)品解釋應(yīng)用系統(tǒng)及預(yù)測技巧 [J]. 應(yīng)用氣象學(xué)報, 21 (5): 614–620. Qin Zhinian, Chen Lijuan, Tang Hongyu, et al. 2010. Downscaling methods and application system based on monthly-scale dynamical model outputs and forecast skill analysis [J]. J. Appl. Meteor. Sci. (in Chinese), 21 (5): 614–620, doi:10.3969/j.issn.1001-7313.2010.05.011.

Raisanen J, Raty O. 2013. Projections of daily mean temperature variability in the future: Cross-validation tests with ENSEMBLES regional climate simulations [J]. Clim. Dyn., 41 (5–6): 1553–1568, doi:10.1007/s00382- 012-1515-9.

Saha S, Moorthi S, Pan H L, et al. 2010. The NCEP climate forecast system reanalysis [J]. Bull. Amer. Meteor. Soc., 91 (8): 1015–1057, doi:10.1175/ 2010BAMS3001.1.

Saha S, Moorthi S, Wu X R, et al. 2014. The NCEP climate forecast system version 2 [J]. J. Climate, 27 (6): 2185–2208, doi:10.1175/JCLI-D-12- 00823.1.

譚桂容, 陳海山, 孫照渤, 等. 2010. 2008年1月中國低溫與北大西洋濤動和平流層異?；顒拥穆?lián)系 [J]. 大氣科學(xué), 34 (1): 175–183. Tan Guirong, Chen Haishan, Sun Zhaobo, et al. 2010. Linkage of the cold event in January 2008 over China to the North Atlantic Oscillation and stratospheric circulation anomalies [J]. Chinese J. Atmos. Sci. (in Chinese), 34 (1): 175–183, doi:10.3878/j.issn.1006-9895.2010.01.16.

王凌, 高歌, 張強(qiáng), 等. 2008. 2008年1月我國大范圍低溫雨雪冰凍災(zāi)害分析Ⅰ. 氣候特征與影響評估 [J]. 氣象, 34 (4): 95–100. Wang Ling, Gao Ge, Zhang Qiang, et al. 2008. Analysis of the severe cold surge, ice-snow and frozen disasters in South China during January 2008: Ⅰ. Climatic features and its impact [J]. Meteor. Mon. (in Chinese), 34 (4): 95–100, doi:10.7519/j.issn.1000-0526.2008.4.012.

魏鳳英. 2008a. 現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)（第二版）[M]. 北京: 氣象出版社, 19–20. Wei Fengying. 2008a. Modern Climate Statistical Diagnosis and Prediction Technology (Second edition) [M] (in Chinese). Beijing: China Meteorological Press, 19–20.

魏鳳英. 2008b. 中國冬半年最低氣溫概率分布特征 [J]. 氣候變化研究進(jìn)展, 4 (1): 8–11. Wei Fengying. 2008b. Probability distribution of minimum temperature in winter half years in China [J]. Advances in Climate Change Research (in Chinese), 4 (1): 8–11, doi:10.3969/j.issn. 1673-1719.2008.01.002.

吳統(tǒng)文, 宋連春, 劉向文, 等. 2013. 國家氣候中心短期氣候預(yù)測模式系統(tǒng)業(yè)務(wù)化進(jìn)展[J]. 應(yīng)用氣象學(xué)報, 24 (5): 533–543. Wu Tongwen, Song Lianchun, Liu Xiangwen, et al. 2013. Progress in developing the short-range operational climate prediction system of China National Climate Center [J]. J. Appl. Meteor. Sci. (in Chinese), 24 (5): 533–543, doi:10.3969/j.issn.1001-7313.2013.05.003.

張慶云, 宣守麗, 彭京備. 2008. La Ni?a年冬季亞洲中高緯環(huán)流與我國南方降雪異常關(guān)系 [J]. 氣候與環(huán)境研究, 13 (4): 385–394. Zhang Qingyun, Xuan Shouli, Peng Jingbei. 2008. Relationship between Asian circulation in the middle–high latitude and snowfall over South China during La Ni?a events [J]. Climatic Environ. Res. (in Chinese), 13 (4): 385–394, doi:10.3878/j.issn.1006-9585.2008.04.04.

鄭志海. 2013. 月動力延伸預(yù)報研究進(jìn)展回顧 [J]. 氣象科技進(jìn)展, 3 (1): 25–30. Zheng Zhihai. 2013. Review of the progress of dynamical extended-range forecasting studies [J]. Adv. Meteor. Sci. Technol. (in Chinese), 3 (1): 25–30, doi:10.3969/j.issn.2095-1973.2013.01.004.

Bias Correction in Monthly Means of Temperature Predictions of the Dynamic Extended Range Forecast Model

ZHANG Daquan1and CHEN Lijuan1, 2

1,,100081;2,,210044

Bias analysis of hindcast data from the monthly Dynamic Extended Range Forecast model (DERF2.0) indicates that the monthly mean temperature prediction performance is not good enough to be used operationally and needs to be improved. Quantile mapping of the non-parameter method is applied to correct the DERF2.0 model bias of monthly mean temperature. The first key step of this method is to utilize deterministic model output and calculate the cumulative density function (CDF) of hindcast data. Then, we can obtain the quantile mapping of the deterministic model output on the CDF of hindcast data. The second step is to calculate the CDF of observation data and map the quantile result of model output to the CDF of observational data. The model bias can be reduced to a certain extent after the above procedures. Hindcast verification shows that the method can significantly reduce the root mean square error (RMSE) of model output and improve the predictive skill of spatial distributions of monthly mean temperature anomalies (1983–2012). Prediction skill (PS) and anomaly correlation coefficient (ACC) scores between model output and observations with different lead times have been improved, and this improvement remains stable with different magnitudes of model bias. Comparison of model prediction skills of different months shows that the enhancement of prediction performance in summer is greater than in winter.

Monthly Dynamic Extended Range Forecast model, Verification, Probability, Bias correction

1006-9895(2016)05-1022-11

P456

A

10.3878/j.issn.1006-9895.1509.15228

2015-07-10；網(wǎng)絡(luò)預(yù)出版日期2015-09-23

章大全，男，1985年出生，博士，主要從事短期氣候預(yù)測研究。E-mail: zhangdq@cma.gov.cn

國家重點基礎(chǔ)研究發(fā)展計劃項目 2012CB955203，公益性行業(yè)（氣象）科研專項GYHY201306032、GYHY201406022，國家自然科學(xué)基金項目41205039，中國氣象局短期氣候預(yù)測創(chuàng)新團(tuán)隊項目

Funded by The National Basic Research Program of China (Grant 2012CB955203),Special Scientific Research Fund of Meteorological Public Welfare Profession of China (Grants GYHY201306032, GYHY201406022), National Natural Science Foundation of China (NSFC) (Grant 41205039), Program for Innovation Team in Short-term Climate Prediction of China Meteorological Administration