降水融合格點產(chǎn)品在陜西2019年暴雨過程中的檢驗

劉菊菊，陳小婷，肖貽青，馬曉華，梁 綿

(1.陜西省氣象臺，西安 710014；2.陜西省氣象局秦嶺和黃土高原生態(tài)環(huán)境氣象重點實驗室，西安 710016)

隨著社會發(fā)展，各行各業(yè)對格點化的時、空連續(xù)的氣象數(shù)據(jù)產(chǎn)品要求越來越高[1]，多源數(shù)據(jù)融合氣象格點產(chǎn)品可滿足以上需求。國家氣象信息中心于2017年6月開始下發(fā)高時空分辨率的國家智能網(wǎng)格實況分析產(chǎn)品，目前對其準確性和真實性的了解相對有限[2]。降水的智能網(wǎng)格實況分析產(chǎn)品包括二源和三源融合格點降水產(chǎn)品。其中，二源融合格點降水產(chǎn)品融合了地面、衛(wèi)星兩種降水數(shù)據(jù)[3]，三源融合格點降水產(chǎn)品融合了地面、衛(wèi)星、雷達三種降水產(chǎn)品的數(shù)據(jù)[4]。由于衛(wèi)星、雷達的降水產(chǎn)品精度較低，因而有必要對融合產(chǎn)品的適用性進行檢驗評估。陜西省地形復(fù)雜，暴雨是陜西夏季主要災(zāi)害性天氣之一，融合降水產(chǎn)品在陜西暴雨過程中的適用性更值得探究。

高分辨率格點預(yù)報是目前中國氣象局主推業(yè)務(wù)和未來天氣預(yù)報的發(fā)展方向[5]。在國內(nèi)外大力發(fā)展無縫隙精細化預(yù)報的背景下，格點實況產(chǎn)品的研究受到重視[6]。目前，針對格點預(yù)報檢驗方法，業(yè)務(wù)中主要采用評估降水預(yù)報準確率的TS評分、漏報率、空報率等傳統(tǒng)檢驗指標[7-10]。陳小婷等[11]用TS評分、預(yù)報偏差對2014年6—9月歐洲中心集合預(yù)報10種降水產(chǎn)品在陜西的預(yù)報效果進行了分級檢驗以及時空演變綜合評估。盧小鳳等[12]采用傳統(tǒng)檢驗指標對2019年前汛期和后汛期西江流域的日面雨量預(yù)報進行了檢驗評估。李明娟等[13]基于傳統(tǒng)檢驗指標，對2019年4—10月陜西省暴雨天氣過程的省臺主觀訂正、ECMWF細網(wǎng)格和模式動態(tài)交叉取優(yōu)要素預(yù)報方法的預(yù)報產(chǎn)品進行檢驗評估。此外，也有學(xué)者利用基于對象診斷的空間檢驗方法MODE(method for object-based diagnostic evaluation)檢驗降水的空間準確率，俞劍蔚等[2]評估了國家級格點實況分析產(chǎn)品在江蘇地區(qū)的適用性，認為格點實況產(chǎn)品較準確地反映了降水空間分布，隨著降水量級增大，格點實況產(chǎn)品的降水漏報較多。劉湊華等[14]對MODE檢驗方法中目標識別、目標配對算法進行了改進，并應(yīng)用到T639模式降水預(yù)報檢驗中。

本文利用多種預(yù)報檢驗指標檢驗高分辨率二源、三源融合降水分析實時產(chǎn)品在陜西省2019年暴雨天氣過程中的適用性，分析降水誤差隨海拔高度的變化特征，以期得到不同地形和降水量級下融合降水產(chǎn)品的評估特征，為精細化格點預(yù)報和訂正提供可靠的基礎(chǔ)數(shù)據(jù)，發(fā)揮其在氣象防災(zāi)減災(zāi)工作中的基礎(chǔ)性作用。

1 資料和方法

1.1 融合產(chǎn)品資料

研究分析的數(shù)據(jù)資料為國家氣象信息中心下發(fā)的高時空分辨率中國區(qū)域融合降水分析系統(tǒng)(CMA multi-source precipitation analysis system，CMPAS-V2.1)的二源和三源降水融合分析實時產(chǎn)品，空間分辨率為0.05°×0.05°。其中，二源融合降水產(chǎn)品主要利用相對成熟的“概率密度匹配+最優(yōu)插值”兩步數(shù)據(jù)融合方法，對經(jīng)過質(zhì)量控制后的全國3萬多個氣象站觀測的小時降水量和國家衛(wèi)星氣象中心實時發(fā)布的1 小時FY2系列衛(wèi)星反演降水產(chǎn)品進行融合，生成了逐小時降水產(chǎn)品[3]。三源降水融合分析實時產(chǎn)品首先用概率密度函數(shù)匹配法訂正雷達估測、衛(wèi)星反演降水產(chǎn)品的系統(tǒng)性偏差；其次采用貝葉斯模型平均法將雷達和衛(wèi)星降水產(chǎn)品結(jié)合形成覆蓋中國的背景場；最后用最優(yōu)插值方法融入地面觀測數(shù)據(jù)[4]，其地面觀測資料采用了中國近4萬個自動氣象站逐時降水量數(shù)據(jù)，該數(shù)據(jù)經(jīng)過了實時質(zhì)量控制[15]。衛(wèi)星反演降水資料選取美國氣候預(yù)測中心研發(fā)的實時衛(wèi)星反演降水系統(tǒng) CMORPH(CPC MOR-PHing technique)的降水產(chǎn)品[16]。雷達估測降水資料采用中國氣象局氣象探測中心研制的雷達定量降水估測產(chǎn)品，根據(jù)統(tǒng)計不同雷達實時反射率因子與雨強關(guān)系的方法估測降水，并采用卡爾曼濾波、平均校準等方法對2 400個國家級自動氣象站降水數(shù)據(jù)進行校正[17]。

1.2 站點實況資料

將融合實況產(chǎn)品作為某種意義上的預(yù)報場，與自動氣象站實況觀測進行對比，將站點觀測值作為大氣狀態(tài)的“真值”，從而檢驗評估融合實況產(chǎn)品。所采用的站點觀測資料來源于同期的陜西省383個地面自動氣象站(99個國家站和經(jīng)過質(zhì)量控制并剔除缺測站后的284個區(qū)域站 )的08—08時(北京時)日降水量數(shù)據(jù)。

1.3 方法

融合產(chǎn)品中融合了陜西省地面自動氣象站資料(國家站和經(jīng)過質(zhì)量控制的區(qū)域站)，因此本文采用非獨立樣本檢驗方法，通過與383個氣象站資料對比，評估融合實況產(chǎn)品的準確性。選取陜西省2019年出現(xiàn)的8次區(qū)域性(≥5站(國家站))暴雨過程，對融合降水產(chǎn)品進行檢驗分析。在檢驗中，格點到站點的轉(zhuǎn)換選用領(lǐng)域法[18]，即以383個站點的實況數(shù)據(jù)作為標準，將所有站的經(jīng)、緯度按照融合實況產(chǎn)品格點數(shù)據(jù)的分辨率進行近似處理，選擇最接近站點的格點值作為檢驗值。檢驗降水的評估指標有等級平均誤差(Ek)、均方根誤差(R)、等級均方根誤差(Rk)、等級TS評分(Tk)、等級降水準確率(Pk)，具體計算方法如下。

(1)

(2)

(3)

(4)

(5)

其中，Oi為站點觀測值，Gi為融合產(chǎn)品插值到檢驗站點得到的數(shù)值，N為參與檢驗的總樣本數(shù)(站次數(shù))。k代表降水分級檢驗等級，Uk為第k個降水等級區(qū)間的上界，Lk為第k個降水等級區(qū)間的下界。式中[]代表邏輯轉(zhuǎn)數(shù)值的運算符，邏輯值為正時取1，否則取0。i代表第i個樣本。

2 降水實況分析

2019年陜西省共出現(xiàn)8次區(qū)域性暴雨過程，分別是6月27日、7月21日、7月28日、8月2日、8月3日、8月8日、9月9日和9月14日(表1)。由于大暴雨的樣本數(shù)太少，將其并入暴雨，即降水量≥50 mm的均為暴雨樣本。其中6月27日暴雨過程主要出現(xiàn)在陜南中西部，共出現(xiàn)暴雨43站，最大降雨量為105.9 mm。7月21日暴雨主要出現(xiàn)在陜北地區(qū)和寶雞、漢中西部，出現(xiàn)暴雨36站，其中最大降水量為108.2mm。7月28日陜北北部有區(qū)域性暴雨，關(guān)中和陜南則出現(xiàn)了局地分散暴雨，最大降雨量108.3 mm。8月2日暴雨主要位于關(guān)中陜南地區(qū)，最大降雨量為101.9 mm。8月3日暴雨主要集中在陜北，陜南有分散暴雨，最大降水量為商洛商南131.3 mm。8月8日暴雨主要出現(xiàn)在關(guān)中東部和陜南西部，最大降雨量為107.1 mm。9月9日暴雨主要出現(xiàn)在陜南，共出現(xiàn)暴雨54站，最大降雨量為139.0 mm。 9月14日暴雨主要出現(xiàn)在關(guān)中和陜南中部，共80站達暴雨，最大降雨量為106.2 mm。7月28日、8月2日、8月8日暴雨落區(qū)比較分散，其他過程暴雨落區(qū)相對集中。依據(jù)降水性質(zhì)，將8次暴雨過程分為對流性降水、穩(wěn)定性降水以及混合性降水(穩(wěn)定性降水過程中有對流發(fā)展)過程3類[19]。8次暴雨過程大部分為混合性降水，挑選7月28日、9月14日分別作為對流性、穩(wěn)定性暴雨過程的代表，與8次過程進行對比檢驗分析。

表1 2019年陜西區(qū)域性暴雨(08時—08時)個例信息

3 降水檢驗分析

3.1 空間對比

通過對比二源、三源融合降水產(chǎn)品降水量空間分布與站點實況的降水量分布(圖略)可知，格點產(chǎn)品基本能反映出8次暴雨過程的主要降水落區(qū)，二者的降水量空間分布基本一致，但是強降水落區(qū)有偏差。二源產(chǎn)品與站點實況相比，6月27日二源產(chǎn)品的中雨范圍在陜北北部偏小，南部偏大，陜南東南部暴雨范圍偏小。二源產(chǎn)品7月21日在關(guān)中陜南的暴雨范圍偏大。7月28日在陜北北部的暴雨范圍偏小。8月2日二源產(chǎn)品在延安南部的暴雨量級偏小，表現(xiàn)為大雨。8月3日二源產(chǎn)品在陜北西部的大雨量級偏小，關(guān)中及陜南西部局地的大雨和暴雨偏小一個量級。8月8日二源產(chǎn)品在陜北的小雨范圍偏大，關(guān)中中部的中雨、陜南西部和關(guān)中東部的暴雨范圍偏小。9月14日二源產(chǎn)品在陜北的小雨范圍偏大；關(guān)中陜南的暴雨范圍略偏小。三源產(chǎn)品與二源產(chǎn)品降水量空間分布基本一致，除了三源產(chǎn)品7月28日在陜南中西部暴雨的落區(qū)更接近實況，9月9日二源融合產(chǎn)品降水空間分布比三源產(chǎn)品更接近降水實況，三源產(chǎn)品對于陜北西部的大雨局地偏大。兩種融合產(chǎn)品和站點實況的8次平均降水量空間分布基本一致，均表現(xiàn)為陜西大部為中雨，陜南西部以及關(guān)中和陜北局地為大雨(圖1，見第5頁)。

圖1 2019年陜西8次暴雨過程降水量(單位為mm)空間分布(a、b、c分別表示8次過程平均和7月28日、9月14日站點實況降水量；d、e、f分別表示8次過程平均和7月28日、9月14日二源融合降水產(chǎn)品降水量；g、h、i分別表示8次過程平均和7月28日、9月14日三源融合降水產(chǎn)品降水量；審圖號為GS(2019)3333)(文見第4頁)

3.2 等級平均誤差和等級降水準確率

由表2可知，二源和三源融合降水產(chǎn)品的小雨、中雨的平均誤差值為正，大雨、暴雨的平均誤差值為負值，表明相對于站點實況，融合降水產(chǎn)品的小雨和中雨數(shù)值偏大，大雨和暴雨的數(shù)值偏小。并且隨著降水量級的增大，平均誤差的絕對值明顯增加。平均而言，三源產(chǎn)品的大雨和暴雨平均誤差略小于二源，小雨和中雨二源產(chǎn)品的平均誤差略小于三源產(chǎn)品。對流性暴雨過程(7月28日)的小雨、大雨、暴雨誤差遠大于其他過程，而穩(wěn)定性過程(9月14日)的大雨和暴雨誤差遠小于其他過程?？梢娙诤袭a(chǎn)品對于穩(wěn)定性過程的表現(xiàn)比對流性過程更優(yōu)。

表2 二源、三源融合降水產(chǎn)品在2019年陜西8次暴雨過程的等級平均誤差 單位：mm

二源產(chǎn)品準確率中雨最大，其次為小雨、大雨、暴雨，三源產(chǎn)品的準確率隨著降水量級的增大而減小。三源融合產(chǎn)品的小雨和大雨準確率高于二源產(chǎn)品，中雨準確率低于二源產(chǎn)品，暴雨準確率二者相差很小。對流性暴雨過程(7月28日)的小雨二源產(chǎn)品準確率更高，三源產(chǎn)品的大雨準確率更高，兩類融合產(chǎn)品的中雨和暴雨準確率一樣。穩(wěn)定性暴雨過程(9月14日)的小雨三源降水產(chǎn)品準確率更高，中雨和暴雨準確率二源產(chǎn)品明顯大于三源產(chǎn)品，大雨準確率二者相差很小。穩(wěn)定性過程的各降水量級(除了三源產(chǎn)品中雨外)準確率遠大于對流性暴雨過程(圖2，見第5頁)。可見，當降水量較大時，三源產(chǎn)品在對流性暴雨過程中的表現(xiàn)更優(yōu)，穩(wěn)定性過程可更多參考二源產(chǎn)品。

圖2 二源(a)和三源(b)降水融合產(chǎn)品不同等級降水準確率(文見第4頁)

3.3 等級TS評分

由圖3見第6頁可知，總體來看，8次暴雨過程的各量級降水三源產(chǎn)品的平均TS評分均高于二源產(chǎn)品。TS評分隨著降水量級的增大而減小，小雨TS評分最高，三源產(chǎn)品8次平均小雨TS評分達87.6%；暴雨TS評分最低，二源產(chǎn)品8次平均暴雨TS評分為56.7%。對流性降水過程(7月28日)的中雨、大雨和暴雨的TS評分明顯小于其他過程。穩(wěn)定性暴雨過程(9月14日)的中雨和大雨TS評分為8次暴雨過程中的最高。小雨的TS評分對流性暴雨過程(7月28日)大于穩(wěn)定性暴雨過程(9月14日)，三源產(chǎn)品二者差值達12.6%。其他量級降水的TS評分穩(wěn)定性過程明顯大于對流性暴雨過程，三源產(chǎn)品的中雨、大雨、暴雨在7月28日和9月14日的TS評分差值分別達7.9%、25.7%和22.3%。8次暴雨過程的平均TS評分介于穩(wěn)定性和對流性暴雨過程之間(圖4，見第6頁)。二源和三源產(chǎn)品的TS評分差值在對流性大雨和穩(wěn)定性暴雨較大，其他降水量級的差異很小。對流性大雨和暴雨三源產(chǎn)品TS評分高于二源產(chǎn)品，穩(wěn)定性的大雨和暴雨二源產(chǎn)品優(yōu)于三源產(chǎn)品?？梢姡诤袭a(chǎn)品在對流性和穩(wěn)定性暴雨過程中TS評分的差異較大，對于對流性過程的小雨和穩(wěn)定性過程的中雨及其以上量級的降水TS評分較高。在以后的應(yīng)用中，對于大雨及以上量級降水，對流性過程宜參考三源產(chǎn)品，穩(wěn)定性過程宜參考二源產(chǎn)品。

圖3 2019年陜西8次暴雨過程二源(a)、三源(b)融合降水產(chǎn)品等級TS評分(文見第4頁)

圖4 2019年陜西8次暴雨過程中降水融合產(chǎn)品的平均等級TS評分與對流性和穩(wěn)定性暴雨過程的等級TS評分(a 二源融合降水產(chǎn)品，b 三源融合降水產(chǎn)品)(文見第4頁)

3.4 均方根誤差

整體來看，降水的均方根誤差在陜北北部、關(guān)中南部和陜南大部地方偏大，部分地區(qū)均方根誤差大于10 mm，而陜北南部和關(guān)中北部均方根誤差較小，大部地區(qū)均方根誤差小于4 mm(圖5)。

圖5 2019年8次暴雨過程降水的均方根誤差(審圖號：GS(2019)3333)

由不同等級降水的等級均方根誤差可知(圖6)，中雨的均方根誤差最小，其次是小雨、大雨、暴雨。8次暴雨過程的平均均方根誤差二源和三源產(chǎn)品的差異很小。二源降水產(chǎn)品小雨的均方根誤差8月2日最小，為4.12 mm，9月9日最大，為13.71 mm。中雨的均方根誤差9月9日最小，為4.31 mm，8月8日最大，為9.69 mm。大雨和暴雨的均方根誤差9月14日最小，分別為6.22 mm和12.24 mm，可能是由于此次暴雨過程為穩(wěn)定性降水過程，造成暴雨的天氣系統(tǒng)較為穩(wěn)定。7月28日大雨和暴雨誤差相對其他過程明顯偏大，可能是由于造成暴雨的直接原因是不穩(wěn)定的中尺度對流系統(tǒng)導(dǎo)致。三源產(chǎn)品的均方根誤差特征與二源產(chǎn)品基本一致。綜上，與其他過程相比，大雨和暴雨的均方根誤差在9月14日穩(wěn)定性降水過程最小，7月28日的對流性過程相對偏大。

圖6 2019年陜西暴雨過程的等級均方根誤差(a二源融合降水產(chǎn)品，b三源融合降水產(chǎn)品)

由表3可知，小雨的均方根誤差對流性暴雨過程(7月28日)大于穩(wěn)定性暴雨過程(9月14日)。中雨的均方根誤差二者差異較小。大雨和暴雨的均方根誤差對流性暴雨過程明顯大于穩(wěn)定性暴雨過程，三源產(chǎn)品二者差值分別達4.57 mm和14.29 mm?？梢?，融合產(chǎn)品對于穩(wěn)定性暴雨過程的表現(xiàn)優(yōu)于對流性暴雨過程，尤其是大雨和暴雨量級降水的優(yōu)勢更明顯。隨著降水量級增大，均方根誤差增加，對流性過程的暴雨均方根誤差最大。

表3 8次暴雨過程的平均等級均方根誤差與對流性和穩(wěn)定性暴雨過程等級均方根誤差 單位：mm

4 海拔高度對降水誤差的影響

由以上分析可知三源降水產(chǎn)品優(yōu)于二源降水產(chǎn)品，下面重點分析三源降水產(chǎn)品的誤差與海拔高度的關(guān)系，按100 m的間距[20]，統(tǒng)計陜西383個站點在各個海拔高度區(qū)間的站點數(shù)分布。由圖7可知，觀測站點主要集中在200～1 700 m海拔高度范圍內(nèi)，故在200～1 700 m海拔高度內(nèi)，按每100 m間隔來分析降水誤差隨海拔高度的變化特征。

圖7 陜西不同海拔高度站點數(shù)分布

由表4可知，降水平均絕對誤差值與海拔高度的相關(guān)性表現(xiàn)為：6月27日、8月8日、9月14日三次暴雨過程降水絕對誤差與海拔高度呈顯著負相關(guān)，7月21日、8月3日兩次暴雨過程降水絕對誤差與海拔高度呈顯著正相關(guān)。7月28日和8月2日降水絕對誤差與海拔高度相關(guān)性較小。降水量與海拔高度相關(guān)系數(shù)同樣表現(xiàn)為6月27日、8月8日、9月14日顯著負相關(guān)，7月21日、8月3日顯著正相關(guān)。其原因可能是降水絕對誤差與海拔高度呈顯著負相關(guān)的3次暴雨過程的暴雨落區(qū)主要位于高海拔的秦巴山區(qū)，而顯著正相關(guān)的2次暴雨過程暴雨落區(qū)主要位于海拔相對較低的陜北地區(qū)。由此可見，暴雨落區(qū)隨海拔高度的變化是降水誤差隨海拔高度變化的主導(dǎo)因素。降水量與降水絕對誤差為顯著正相關(guān)，進一步說明降水量級是降水誤差的主要影響因素。

表4 2019年陜西8次暴雨過程不同要素之間的相關(guān)系數(shù)

挑選降水絕對誤差與海拔高度顯著相關(guān)的過程來分析不同等級降水均方根誤差隨海拔高度變化的特點。對于7月21日、8月3日兩次降水絕對誤差與海拔高度呈顯著正相關(guān)的暴雨過程，小雨和中雨的均方根誤差隨海拔高度變化較小，大雨和暴雨均方根誤差隨海拔高度呈增加趨勢?？梢姡饕谴笥旰捅┯昃礁`差隨海拔高度增大，導(dǎo)致2次暴雨過程降水均方根誤差與海拔高度顯著正相關(guān)。對于6月27日、8月8日、9月14日這三次降水絕對誤差與海拔高度呈顯著負相關(guān)的暴雨過程，小雨和暴雨均方根誤差隨海拔高度減小，中雨和大雨均方根誤差隨海拔高度變化趨勢不明顯。主要是小雨和暴雨均方根誤差隨海拔高度減小，導(dǎo)致3次暴雨過程降水誤差與海拔高度顯著負相關(guān)。降水均方根誤差與海拔高度顯著正、負相關(guān)的主要影響因素均有暴雨隨海拔高度的變化?？梢?，主要是暴雨誤差隨海拔高度的變化趨勢決定降水誤差與海拔高度的關(guān)系。

5 結(jié)論

(1)融合產(chǎn)品的降水量空間分布與站點實況空間分布基本一致，但強降水落區(qū)有所偏差。相對于站點實況，融合產(chǎn)品的小雨和中雨數(shù)值偏大，大雨和暴雨的數(shù)值偏小。

(2)融合產(chǎn)品的準確率隨著降水量級的增大而減小，平均誤差的絕對值和TS評分隨著降水量級的增大而增加。中雨的均方根誤差最小，其次是小雨、大雨、暴雨。二源產(chǎn)品降水準確率中雨最大，其次為小雨、大雨、暴雨。三源融合產(chǎn)品的小雨和大雨準確率高于二源產(chǎn)品。

(3)穩(wěn)定性暴雨過程(9月14日)的降水準確率遠大于對流性暴雨過程(7月28日)。平均誤差絕對值、均方根誤差和TS評分均表明，對于大雨和暴雨，融合產(chǎn)品在穩(wěn)定性暴雨過程的準確率明顯優(yōu)于對流性過程。

(4)三源融合產(chǎn)品的平均TS評分高于二源融合產(chǎn)品。對流性暴雨過程大雨及以上量級降水可參考三源融合產(chǎn)品，穩(wěn)定性過程的大雨及以上量級降水宜參考二源融合產(chǎn)品。

(5)暴雨落區(qū)隨海拔高度的變化影響降水誤差隨海拔高度變化。當暴雨落區(qū)主要位于高海拔地區(qū)時，降水誤差與海拔高度呈顯著負相關(guān)；而暴雨落區(qū)主要位于海拔相對較低地區(qū)時，降水誤差與海拔高度呈顯著正相關(guān)。