新增積雪深度客觀預(yù)報技術(shù)研究及其應(yīng)用*

符嬌蘭 陳博宇 陳 雙 董 全 曹 勇 楊舒楠

國家氣象中心，北京 100081

提 要：利用地面歷史觀測數(shù)據(jù)對我國冬季典型降雪過程的雪水比進行了統(tǒng)計分析，在此基礎(chǔ)上，利用改進后的Cobb算法、ECMWF IFS模式預(yù)報及中央氣象臺降水和相態(tài)網(wǎng)格預(yù)報構(gòu)建了雪水比和新增積雪深度物理預(yù)報模型。結(jié)果表明：我國降雪過程的雪水比分布較廣，最小為0.3∶1,最大為100∶1,僅有4%左右的雪水比為10∶1；雪水比具有明顯的區(qū)域和季節(jié)特征，新疆、西北地區(qū)東部、華北以及東北地區(qū)雪水比大于其他區(qū)域，冬季較春、秋季雪水比大；改進后的Cobb算法建立了隨地形高度變化的云識別方案，利用云內(nèi)溫度與雪水比統(tǒng)計關(guān)系及垂直速度構(gòu)建了雪水比和新增積雪深度預(yù)報模型，其預(yù)報的雪水比和新增積雪深度空間分布特征較原Cobb算法的預(yù)報更接近實況；結(jié)合雪水比預(yù)報模型和中央氣象臺降雪網(wǎng)格技術(shù)的新增積雪深度預(yù)報較ECMWF IFS模式預(yù)報有明顯的改進，相對模式新增積雪深度預(yù)報TS評分提高率基本在10%以上，尤其對20 cm以上新增積雪深度預(yù)報能力提升明顯。

引 言

強降雪天氣形成的積雪是雪災(zāi)關(guān)鍵致災(zāi)因子之一(宮德吉和郝慕玲，1998;時興合等，2006)。2010年新年伊始，我國華北地區(qū)普降大到暴雪，降雪量達到歷史同期第一位(李青春等，2011)。與此同時，1 月3 日20 時北京南郊觀象臺觀測的積雪深度達23.4 cm，突破歷史極值。這場暴雪給人們的生產(chǎn)生活、尤其是交通出行等造成了非常嚴(yán)重的影響。2018年1月3—5日安徽省中北部出現(xiàn)了大暴雪到特大暴雪天氣,全省有51個市縣積雪超過15 cm,合肥積雪達到21 cm(楊祖祥等，2019)。受積雪影響，4日合肥市望江路公交站臺頂板在大雪中倒塌,造成一死多傷,帶來了極大的社會關(guān)注度。

雪水比，即某一時間段內(nèi)(如24 h)降雪天氣產(chǎn)生的新增積雪深度(單位：cm)與相應(yīng)的降雪融化后形成的液態(tài)水量(單位：mm)之間的比值(國際上通常稱為snow to liquid ratio或者snow ratio，簡稱SR)，是降雪過程新增積雪深度預(yù)報所必須的重要參數(shù)(崔錦等，2017)。多年來，業(yè)務(wù)新增積雪深度預(yù)報多以經(jīng)驗預(yù)報為主，預(yù)報員根據(jù)新增積雪深度與降雪量之間的經(jīng)驗關(guān)系(雪水比經(jīng)驗值),將降雪量預(yù)報換算成積雪深度預(yù)報。1878年加拿大的氣象觀測員根據(jù)多倫多安大略湖長期觀測資料分析得到了當(dāng)?shù)仄骄┧葹?0∶1。此后，相當(dāng)長的一段時間內(nèi)，國內(nèi)外氣象業(yè)務(wù)部門利用雪水比10∶1對積雪深度進行簡單的換算預(yù)報。然而，隨著觀測記錄日益增多，分析發(fā)現(xiàn)雪水比10∶1是一個相當(dāng)粗糙的近似值(Cobb and Waldstreicher，2005)。Judson and Doesken(2000)對美國大量降雪過程進行了統(tǒng)計分析，表明雪水比為3∶1～100∶1。不少研究也表明不同降雪天氣過程對應(yīng)的雪水比數(shù)值相差甚遠(Doesken and Juson,1996;Super and Holroyd,1997)。陳雙和符嬌蘭(2021)對華北地區(qū)兩次降雪過程的對比分析也表明不同天氣過程雪水比差異非常大。由此可見，要開展新增積雪深度預(yù)報，就必須弄清楚我國降雪天氣雪水比的氣候分布特征。然而目前我國對雪水比氣候特征的研究工作相對較少，起步較晚。楊琨和薛建軍(2013)僅用了兩個冬季的觀測數(shù)據(jù)進行了研究，缺少不同季節(jié)不同區(qū)域雪水比的統(tǒng)計分析。近年來，崔錦等(2015)、王一頡等(2019)、楊成芳和朱曉晴(2020)、魏凌翔等(2021)對北方部分地區(qū)的雪水比氣候特征及影響因子進行了分析，但缺乏全國范圍內(nèi)的分析。

曹冀魯(2010)指出，我國氣象業(yè)務(wù)部門僅發(fā)布降雪量預(yù)報，并不提供積雪深度預(yù)報產(chǎn)品，這容易造成公眾對預(yù)報信息理解的偏差。為了滿足日益增加的氣象預(yù)報服務(wù)需求，有必要開展精細化積雪深度預(yù)報業(yè)務(wù)，然而，目前我國新增積雪深度預(yù)報技術(shù)支撐相對較薄弱。國際上積雪深度預(yù)報技術(shù)主要有四大類：數(shù)值模式預(yù)報、氣候?qū)W、統(tǒng)計預(yù)報模型以及物理預(yù)報模型(Alcott and Steenburgh，2010)。以歐洲中期天氣預(yù)報中心全球預(yù)報模式(以下簡稱為ECMWF IFS)為例，其積雪深度預(yù)報是陸面物理過程的一部分，通過地表雪水含量和積雪密度計算得到積雪深度,由于目前模式積雪密度采用一層模式，勢必會導(dǎo)致密度偏大、積雪偏小(ECMWF，2020)。Roebber et al(2003；2007)通過分析溫度、濕度、地面風(fēng)場等物理因子與雪水比的統(tǒng)計關(guān)系，利用神經(jīng)網(wǎng)絡(luò)的方法建立了雪水比的統(tǒng)計預(yù)報模型。Cobb and Waldstreicher(2005)利用溫度和雪水比的統(tǒng)計關(guān)系，以垂直速度作為權(quán)重系數(shù)建立了雪水比物理預(yù)報模型，并結(jié)合降雪量預(yù)報開發(fā)了新增積雪深度預(yù)報技術(shù)(以下簡稱Cobb算法)，該技術(shù)已在美國氣象部門得到了很好應(yīng)用。考慮到利用統(tǒng)計模型需要大量的歷史觀測數(shù)據(jù)，2013年我國地面觀測業(yè)務(wù)開始改革，夜間降水相態(tài)等人工觀測業(yè)務(wù)取消，導(dǎo)致近年來降雪量等觀測信息的缺失，因此本文主要參考Cobb算法構(gòu)建新增積雪深度物理預(yù)報模型和客觀預(yù)報技術(shù)。

本文利用地面站觀測數(shù)據(jù)、ECMWF IFS模式預(yù)報數(shù)據(jù)、中央氣象臺智能網(wǎng)格降水、降水相態(tài)預(yù)報等資料對我國典型降雪過程的雪水比氣候特征進行分析，在此基礎(chǔ)上建立了新增積雪深度客觀預(yù)報技術(shù)，并對其預(yù)報效果進行應(yīng)用評估。

1 資料及方法介紹

利用我國1951—2012年國家站地面觀測數(shù)據(jù)集中的降雪量(SN)、風(fēng)速、積雪深度數(shù)據(jù)，基于對典型降雪事件的定義選取了我國典型降雪個例，通過計算當(dāng)日08時與前一日08時積雪量差，獲得當(dāng)日08時24 h累計新增積雪量(NSD)，利用式(1)得到雪水比(SR)，并對其進行氣候統(tǒng)計。此外，利用ECMWF-IFS模式預(yù)報的溫度、相對濕度和垂直速度等物理量(空間分辨率為0.2°，時間分辨率為3 h和6 h)、中央氣象臺智能網(wǎng)格降水和相態(tài)預(yù)報產(chǎn)品(空間分辨率為5 km，時間分辨率為3 h)以及5 km分辨率地形數(shù)據(jù)建立了雪水比和新增積雪深度客觀預(yù)報技術(shù)。進一步利用2019年10月至2020年3月地面積雪觀測、ECMWF-IFS模式積雪深度預(yù)報及新增積雪深度客觀預(yù)報數(shù)據(jù)，通過計算式(2)的TS評分對新增積雪深度客觀預(yù)報產(chǎn)品預(yù)報性能進行評估。

(1)

(2)

式中:Na為預(yù)報正確站點數(shù)，Nb為空報站點數(shù)，Nc為漏報站點數(shù)。

2 我國典型降雪天氣雪水比氣候特征分布

為避免弱降雪觀測不精確以及風(fēng)吹雪等對積雪的影響(Roebber et al,2003)，典型降雪個例選取標(biāo)準(zhǔn)為：24 h累計降雪量≥2.5 mm，累計新增積雪深度≥3 cm，同時地面風(fēng)速<9 m·s-1。需要說明是，盡管對降雪個例進行了篩選，但是我國積雪觀測資料只有08時24 h累計積雪量記錄，中間可能出現(xiàn)融雪等現(xiàn)象，因此可能會導(dǎo)致計算的雪水比比實際值要偏小。通過對我國1951—2012年國家站降水量、積雪深度以及10 m風(fēng)場等觀測數(shù)據(jù)進行篩選，共有800個站點44 567個站次滿足上述標(biāo)準(zhǔn)。根據(jù)式(1)雪水比的計算方法，利用24 h累計降雪量和新增積雪深度觀測數(shù)據(jù)，計算得到所有站次的雪水比，并對所有站次以及558個站點(該站點樣本數(shù)大于10)雪水比進行了時空特征分析。從所有站次雪水比統(tǒng)計特征分布可以看出(圖1)，雪水比分布非常廣，最小值為0.3∶1，最大值為100∶1，僅有1705個站次雪水比為10∶1(9.5∶1～10.4∶1)，占比4%左右，與Roebber et al(2013)的研究結(jié)果基本一致。按照基于雪水比的降雪性質(zhì)分類方法，可將降雪分為干雪個例(SR≥15∶1)、濕雪(SR≤9∶1)以及正常降雪(9∶1

圖1 我國44 567個站次雪水比頻次統(tǒng)計分布

楊琨和薛建軍(2013)分析指出雪水比具有明顯的地區(qū)差異。為此，進一步對雪水比時空分布進行分析。圖2a為冬半年雪水比平均值的空間分布，新疆、西北地區(qū)東部、華北、東北地區(qū)等地雪水比氣候均值≥10∶1，其中內(nèi)蒙古東北部、黑龍江及新疆、甘肅西部等地的部分地區(qū)雪水比平均值>12∶1，最大可達16∶1；青藏高原地區(qū)、黃淮、江淮、江漢以及貴州等地雪水比均值為<10∶1、>7∶1,江南等地大部地區(qū)雪水比<7∶1，云南西北部個別站點雪水比<5∶1。與此同時，從春、秋季以及冬季雪水比平均值分布可以看出(圖2b～2d),雪水比表現(xiàn)出一定的季節(jié)差異，總體而言冬季雪水比最大，其次為春季，秋季最?。欢?，黃淮及其以北大部地區(qū)平均雪水比>10∶1；秋季僅新疆、甘肅中西部、內(nèi)蒙古中東部以及東北地區(qū)大部雪水比>10∶1；此外，華北、黃淮、江淮、江漢以及青藏高原東部等地雪水比季節(jié)特征最顯著。Roebber et al(2007)研究指出，雪水比受大氣溫度、濕度、地面和地表溫度等因素影響。由于我國不同季節(jié)、不同區(qū)域大氣溫度、濕度以及地面和地表溫度等氣象條件存在較大的差異，我國雪水比存在明顯的季節(jié)性以及區(qū)域性特征。基于氣候?qū)W或者經(jīng)驗值的雪水比對新增積雪深度進行預(yù)報將會帶來較大的預(yù)報偏差。

圖2 我國(a)冬半年、(b)秋季、(c)冬季、(d)春季典型降雪天氣雪水比氣候特征分布

3 新增積雪客觀預(yù)報技術(shù)

積雪形成與云內(nèi)、云外以及地面過程等有關(guān)(Roebber et al,2003)，其中云內(nèi)微物理過程對積雪形成至關(guān)重要。Cobb算法建立的雪水比及積雪深度物理預(yù)報模型主要考慮了云內(nèi)微物理過程對積雪深度的影響。研究表明：云內(nèi)-18～-12℃溫度層最有利于樹枝狀冰晶的形成，對應(yīng)的降雪雪密度低，雪水比大，高于或低于該溫度則雪密度增加，雪水比明顯減小(Libbrecht,2006)。在飽和情況下，垂直速度與單位時間降水率直接有關(guān)，大的垂直速度將會把云內(nèi)過冷卻水帶到更高的冰晶層內(nèi)，從而有利于冰晶增長，并進一步影響雪水比，因此基于Cobb算法構(gòu)建的預(yù)報模型重點引入了云內(nèi)溫度、垂直速度等變量計算雪水比(Cobb and Waldstreicher,2005)。本文參考Cobb算法利用ECMWF IFS模式輸出的高空溫度、相對濕度、垂直速度以及高度場等變量建立了雪水比預(yù)報模型，算法具體如下：

首先需要在模式輸出的濕度廓線內(nèi)識別出云，Cobb算法(Barnwell,2011)利用相對濕度>90%來識別云，根據(jù)中央氣象臺多年預(yù)報經(jīng)驗表明，模式預(yù)報的相對濕度達到85%左右就有可能出現(xiàn)降雪，因此本文將云定義為相對濕度≥85%?？紤]到地形的影響，不同海拔高度云識別的上下邊界選擇的氣壓層不一樣，云識別范圍的下邊界氣壓層位勢高度必須≥地形高度(H)，上邊界最高氣壓層選取標(biāo)準(zhǔn)如下：

當(dāng)H≥4000 m，最高氣壓層為100 hPa；

當(dāng)2000 m≤H<4000 m，最高氣壓層為200 hPa；

當(dāng)1000 m≤H<2000 m，最高氣壓層為300 hPa；

當(dāng)H<1000 m，最高氣壓層為400 hPa。

其次，計算云內(nèi)各層雪水比權(quán)重系數(shù)Wi，具體如下：

(3)

式中：ωmax為整個云體內(nèi)最大垂直上升速度，ω是任意云層內(nèi)的平均垂直上升速度，φ2、φ1為云層上下邊界的位勢高度。

SRi(T)為云內(nèi)某一氣壓層上的雪水比，由雪水比和溫度的統(tǒng)計關(guān)系得到。最后，計算整個云內(nèi)加權(quán)的雪水比WSRi，并將各層雪水比累加得到最終雪水比WSR，具體如下：

(4)

WSR=∑WSRi

(5)

由于ECMWF IFS模式空間分辨率為9 km，為了得到更精細的積雪深度預(yù)報產(chǎn)品，利用雙線性插值方法將模式溫度、相對濕度、垂直速度以及高度場插值到5 km分辨率。結(jié)合5 km分辨率地形數(shù)據(jù)以及上述變量，根據(jù)式(3)～式(5)分析得到5 km分辨率的雪水比預(yù)報產(chǎn)品。

在雪水比預(yù)報模型的基礎(chǔ)上，根據(jù)雪水比計算式(1)，用雪水比乘以降雪量即可得到新增積雪深度。由于數(shù)值模式降水量(曹勇等，2016)和相態(tài)預(yù)報(董全等，2020a；2020b)存在一定的預(yù)報偏差，預(yù)報的精細化程度也無法滿足業(yè)務(wù)預(yù)報需求。近年來，中央氣象臺發(fā)展了智能網(wǎng)格降水量和相態(tài)預(yù)報技術(shù)。從預(yù)報效果來看，格點化定量降水預(yù)報TS評分和均方根誤差均較ECMWF IFS模式降水預(yù)報有一定提高(曹勇等，2016)。基于集合預(yù)報系統(tǒng)的最優(yōu)概率閾值法構(gòu)建的相態(tài)預(yù)報技術(shù)較集合預(yù)報系統(tǒng)控制成員及ECMWF IFS模式預(yù)報顯著提高了降水相態(tài)預(yù)報能力(董全等，2020a；2020b)。因此,本文降雪量預(yù)報數(shù)據(jù)來源于中央氣象臺智能網(wǎng)格降水以及相態(tài)預(yù)報產(chǎn)品，通過判斷逐3 h或6 h的降水相態(tài)類型，將出現(xiàn)降雪時段的降水量確定為相應(yīng)時段的降雪量，結(jié)合雪水比預(yù)報模型即可得到相應(yīng)時段的新增積雪深度預(yù)報，具體技術(shù)流程見圖3。根據(jù)上述技術(shù)流程，發(fā)展了雪水比、新增積雪深度預(yù)報技術(shù)，相應(yīng)的產(chǎn)品空間分辨率為5 km，1～3 d預(yù)報時效時間分辨率為3 h、6 h、24 h，4～10 d預(yù)報時效時間分辨率為6 h、24 h。

圖3 新增積雪預(yù)報技術(shù)流程框架圖

考慮到地形的影響，不同于Cobb算法預(yù)報模型(以下簡稱Cobb方案)對云識別范圍定義(云識別范圍為1000～500 hPa)(Cobb and Waldstreicher,2005)，本文所采用的算法(以下簡稱為改進后的Cobb方案)在云識別方案中對不同地形高度區(qū)域的云識別范圍進行了規(guī)定，具體見上文。為了對比不同云識別方案對新增積雪預(yù)報性能的影響，本文選取典型個例進行分析。2018年1月3日，西北地區(qū)東部、華北西部、黃淮西部、江淮、江漢等地出現(xiàn)了大范圍強降雪，大部地區(qū)出現(xiàn)了中到大雪，陜西關(guān)中平原、山西南部、河南西部和南部、安徽中北部、湖北中北部等地出現(xiàn)暴雪到大暴雪，局地特大暴雪(圖4a)。需要說明的是，考慮到有些區(qū)域出現(xiàn)了雨雪相態(tài)轉(zhuǎn)化，由于積雪只有24 h累計量觀測，為了保證分析的準(zhǔn)確性，分析時已將24 h內(nèi)出現(xiàn)了雨雪相態(tài)轉(zhuǎn)化的站點過濾掉。受其影響，上述大部地區(qū)新增積雪深度達到3～10 cm，其中，陜西關(guān)中平原、山西南部、河南西部和南部、湖北北部、安徽中北部等地的部分地區(qū)新增積雪深度高達10～30 cm，局地超過30 cm(圖5a)。通過積雪深度觀測及降雪量計算得到實況24 h雪水比(圖6a)，實況雪水比表現(xiàn)出明顯的區(qū)域差異，河南中南部、湖北中北部、蘇皖北部等地雪水比為3∶1～10∶1，局地大于10∶1，而西北地區(qū)東部、山西中南部、河南北部等地大部分地區(qū)雪水比在10∶1～25∶1，局地大于30∶1。

圖4 (a)2018年1月4日08時觀測的降雪量(彩色點，單位：mm)及(b)中央氣象臺智能網(wǎng)格36 h預(yù)報時效降雪量預(yù)報(填色)

圖5 (a)2018年1月4日08時觀測的24 h累計新增積雪深度(單位：cm)及(b)改進后的Cobb方案,(c)Cobb方案36 h時效24 h累計新增積雪預(yù)報(填色)

圖6 (a)2018年1月4日08時觀測的24 h累計雪水比實況以及(b)改進后的Cobb方案,(c)Cobb方案的雪水比預(yù)報

圖4b顯示，對于上述地區(qū)的降雪，中央氣象臺智能網(wǎng)格預(yù)報的降雪范圍和強度與實況基本一致，僅在湖北西南部等地降雪預(yù)報較實況偏大。湖北西南部降雪偏大的主要原因是智能網(wǎng)格相態(tài)預(yù)報為純雪，而實況該地區(qū)為降雨。對比改進后的Cobb方案與Cobb方案相應(yīng)的新增積雪預(yù)報可見(圖5b,5c)，除了在湖北西南部因相態(tài)預(yù)報偏差出現(xiàn)明顯空報之外，二者預(yù)報的新增積雪與實況基本一致，但改進后的Cobb方案對河南東南部、安徽北部以及湖北北部大于20cm新增積雪深度預(yù)報更接近實況，且在河南等地預(yù)報出了30 cm以上的積雪，與實況相符。Cobb方案對上述地區(qū)的新增積雪深度明顯偏小，對山西中部、河南東北部等地新增積雪深度高估。從TS評分也可以看出(圖7)，改進后的Cobb方案各量級新增積雪深度均較Cobb方案的TS評分提高明顯，尤其對大于10 cm及以上量級的預(yù)報TS評分提高率大于18%。

圖7 2018年1月4日08時36 h時效基于改進后的Cobb方案和Cobb方案新增積雪深度預(yù)報TS評分

對比分析Cobb方案與改進后的Cobb方案雪水比預(yù)報可以發(fā)現(xiàn)(圖6b,6c)，兩個方案預(yù)報的雪水比空間分布與實況較為一致，但Cobb方案對山西、河南北部雪水比預(yù)報明顯較實況偏大，對于河南南部等地雪水比預(yù)報略偏小，而改進后的Cobb方案對河南、山西等地的預(yù)報更接近實況，較為準(zhǔn)確地預(yù)報了河南東南部20 cm以上新增積雪，并減少了山西中部和河南北部積雪深度高估情況。

可見，兩個方案新增積雪深度以及雪水比預(yù)報差異最大主要位于河南東南部等地。為進一步說明二者偏差來源，以河南東南部為例，對其溫度、濕度以及垂直速度的垂直結(jié)構(gòu)進行分析。從ECMWF IFS模式2日20時預(yù)報的物理量廓線可以看出(圖8)，最有利于冰晶增長的溫度層結(jié)(-18～-12℃)位于500～400 hPa，且最大垂直速度基本與該層結(jié)重合，根據(jù)陳雙和符嬌蘭(2021)、Lib-brecht(2006)等研究表明，1 000～400 hPa層內(nèi)溫度和垂直速度的垂直分布更有利于密度較小、雪水比大、積雪深度大的降雪產(chǎn)生，Cobb方案僅考慮1000～500 hPa層內(nèi)的溫度和抬升速度對雪水比的影響，由于該層結(jié)內(nèi)溫度普遍高于-12℃，過冷卻水含量較高，更有利于密度大、雪水比小的降雪產(chǎn)生，這也就是改進后的Cobb方案雪水比較Cobb方案大、預(yù)報的新增積雪深度與實況更為接近的原因。

圖8 ECMWF IFS模式2018年1月2日20時預(yù)報的2日20時至4日20時區(qū)域(32.3°～33.3°N、114.3°～115.8°E)平均溫度(紅色線及填色,單位：℃)、相對濕度(綠色實線,單位：%)及垂直速度(黑色線,單位：Pa·s-1)時間-高度剖面對比

2019年11月17日，吉林大部、黑龍江等地普降大到暴雪，吉林東部、黑龍江東部出現(xiàn)大暴雪，局地特大暴雪，大部地區(qū)降雪量預(yù)報與實況基本一致，黑龍江中西部略偏大(圖略)。圖9顯示，吉林大部、黑龍江大部24 h新增積雪深度超過5 cm，吉林大部、黑龍江東部新增積雪深度10 cm以上，局地達到30 cm?？傮w而言，兩個方案新增積雪深度預(yù)報對吉林中東部和黑龍江東部積雪深度預(yù)報與實況基本一致，基于改進后的Cobb算法準(zhǔn)確預(yù)報了吉林東部局地30 cm積雪，但二者對30 cm積雪范圍預(yù)報較實況偏大。此外，對黑龍江中西部積雪預(yù)報也略偏大；一部分是由于降雪量預(yù)報較實況略偏大(圖略)，另一部分原因則與雪水比預(yù)報偏差有關(guān)。進一步對比兩種方案的雪水比預(yù)報發(fā)現(xiàn)(圖10)，基于Cobb算法的雪水比較實況明顯偏大，尤其是在黑龍江中西部，基于改進后的Cobb算法對雪水比偏大的誤差有明顯改進。從圖10a藍色方框內(nèi)的溫度、濕度及垂直速度剖面可以看出(圖略)，降水期間，云層伸展高度達到了400 hPa及以上，最有利于冰晶增長的溫度層結(jié)位于925～700 hPa，最大抬升速度主要位于500 hPa及以下，同時，黑龍江中西部地形復(fù)雜，大部分地區(qū)地形高度在500 m以上，因此基于改進后的Cobb算法較基于Cobb算法的雪水比小，且更接近實況。從新增積雪深度預(yù)報TS評分可以看出(圖11)，基于改進后的Cobb算法新增積雪深度預(yù)報TS高于基于Cobb算法的新增積雪深度預(yù)報。

圖9 同圖5，但為2019年11月18日08時

圖10 同圖6，但為2019年11月18日08時

圖11 同圖7，但為2019年11月18日08時

綜上所示，基于改進后的Cobb方案相對Cobb方案對雪水比及新增積雪預(yù)報改進明顯。需要指出的是，Cobb方案及改進后的Cobb方案僅考慮了與云內(nèi)微物理過程有關(guān)的溫濕廓線對新增積雪深度的影響，并未考慮融雪等近地面過程。由于目前常規(guī)積雪深度觀測僅有24 h累計量，期間發(fā)生的融雪量無法觀測。因此，本文所探討的雪水比有可能比實際雪水比大，這可能是上述個例中本方案預(yù)報的雪水比在部分地區(qū)較實況雪水比偏大的可能原因。目前，即使是預(yù)報性能最好的ECMWF IFS模式，也不能較好地把降雪伴隨的地面融化過程刻畫好，未來還需要增加地面融雪預(yù)報技術(shù)的研發(fā)。

4 新增積雪客觀預(yù)報產(chǎn)品檢驗評估

為更好評估新增積雪客觀預(yù)報產(chǎn)品預(yù)報性能，利用全國國家站08時實況積雪深度觀測資料對2019年10月至2020年3月期間新增積雪客觀預(yù)報產(chǎn)品以及ECMWF IFS模式新增積雪預(yù)報產(chǎn)品的預(yù)報情況進行了檢驗。用當(dāng)日08時積雪深度觀測量減去前一日08時積雪深度觀測量得到當(dāng)日08時24 h累計新增積雪深度實況資料。ECMWF IFS模式積雪預(yù)報(sd)數(shù)據(jù)是通過積雪當(dāng)量預(yù)報(SWE)以及雪密度D計算得到，具體見式(6)。進一步參考實況新增積雪深度計算方法得到ECMWF IFS模式新增積雪深度預(yù)報數(shù)據(jù)。

(6)

式中：D為降雪密度(單位：kg·m-3)，SWE為積雪深度當(dāng)量(單位：10-7m·kg·m-3)，sd為積雪深度(單位，cm)。

利用雙線性插值方法將ECMWF IFS模式和基于改進后的Cobb算法新增積雪預(yù)報數(shù)據(jù)插值到實況新增積雪站點上，并計算新增積雪TS評分。參照降雪量等級及積雪觀測業(yè)務(wù)規(guī)范，評分等級分別為1、3、5、10、20、30 cm，并進一步按照下式計算了相對ECMWF IFS模式新增積雪預(yù)報產(chǎn)品的TS評分提高率(RIT)。

(7)

式中：Tsa為基于改進后的Cobb算法新增積雪預(yù)報產(chǎn)品TS評分，而Tsb為ECMWF IFS模式新增積雪預(yù)報產(chǎn)品的TS評分。

從提前1～3 d的預(yù)報評分可以看出(圖12)，基于改進后的Cobb算法新增積雪深度預(yù)報產(chǎn)品對絕大部分量級新增積雪深度預(yù)報評分要優(yōu)于ECMWF IFS模式，提前2 d預(yù)報評分提高率基本達到了10%以上(圖12a)，尤其是大于20 cm新增積雪深度預(yù)報，評分相對模式提高率超過100%(圖12b)，可見，該方法相對數(shù)值模式改進明顯。對于30 cm以上新增積雪深度預(yù)報評分，提前1～2 d該方法較模式略有降低。通過分析發(fā)現(xiàn)，2019年10月1日至2020年3月20日期間，30 cm以上新增積雪觀測記錄主要位于西藏邊境線附近，考慮到該地區(qū)觀測極少，模式預(yù)報的降雪強度通常較觀測強，易造成空報，中央氣象臺智能網(wǎng)格降水預(yù)報對此處降雪量進行了一定程度的消空，可能造成對個別強降雪過程的漏報，從而導(dǎo)致新增積雪深度預(yù)報偏小。

圖12 2019年10月1日至2020年3月20日短期時效(36～84 h)(a)基于改進后的Cobb算法和ECMWF IFS模式新增積雪深度預(yù)報TS評分，(b)改進后的Cobb算法相對于ECMWF IFS模式預(yù)報TS評分提高率

5 結(jié)論與討論

本文利用歷史觀測數(shù)據(jù)對我國典型降雪過程的雪水比氣候特征進行了統(tǒng)計研究，在此基礎(chǔ)上，構(gòu)建了基于改進后的Cobb算法的雪水比及新增積雪深度客觀預(yù)報技術(shù)，并對其進行了應(yīng)用評估，具體結(jié)論如下：

(1)我國典型降雪過程的雪水比統(tǒng)計分布較廣，最小值為0.3∶1，最大值為100∶1，僅有4%左右的雪水比為10∶1，雪水比位于9∶1～15∶1的樣本占比也僅有30%左右。雪水比具有明顯的季節(jié)和區(qū)域特征，北方地區(qū)雪水比較南方地區(qū)大，冬季雪水比較秋、春季大。

(2)基于改進后的Cobb方案建立了雪水比及新增積雪預(yù)報技術(shù)。較原Cobb方案，改進后的Cobb方案考慮了地形對不同區(qū)域降雪云邊界識別及云內(nèi)溫度及垂直抬升速度廓線的影響，檢驗表明改進后的Cobb方案雪水比和新增積雪預(yù)報性能較原方案改進明顯。

(3)基于改進后的Cobb算法新增積雪預(yù)報性能較好，準(zhǔn)確率較ECMWF IFS模式預(yù)報有明顯提高，短期時效內(nèi)(1～3 d)大部分量級預(yù)報提高率均在10%以上，尤其是對于20 cm以上的新增積雪深度，提高率可達100%及以上。

本文利用改進后的Cobb算法構(gòu)建了雪水比和新增積雪預(yù)報技術(shù)，該技術(shù)已在中央氣象臺得到很好應(yīng)用，并加入了“智慧冬奧2022天氣預(yù)報示范計劃”，為冬季強降雪預(yù)報以及2022年北京冬奧會賽事氣象服務(wù)保障提供了有力的技術(shù)支撐。應(yīng)用評估表明，盡管該預(yù)報技術(shù)較模式積雪預(yù)報有明顯提升，但仍有較大的改進空間，一方面降水相態(tài)預(yù)報難度較大，當(dāng)出現(xiàn)較復(fù)雜的相態(tài)轉(zhuǎn)換過程，降水相態(tài)預(yù)報可能會出現(xiàn)一定的誤差，進而會導(dǎo)致新增積雪預(yù)報出現(xiàn)較大偏差；此外，降雪過程近地面溫度或地表溫度高于0℃，就會出現(xiàn)一定的融雪現(xiàn)象，而目前研發(fā)的新增積雪預(yù)報技術(shù)并未考慮到地面融雪過程，這也會影響到新增積雪預(yù)報準(zhǔn)確率(胡寧等，2021；楊成芳和劉暢，2019)。因此，未來需要進一步提升降水相態(tài)預(yù)報性能，同時，探索近地面融雪過程對新增積雪深度的影響及相應(yīng)的預(yù)報技術(shù)，從而更好地提高新增積雪預(yù)報準(zhǔn)確率。