楊榮芳，楊彬云，康曉甫

（1.河北省氣象技術(shù)裝備中心，河北 石家莊 050021；2.商洛市氣候適應(yīng)型城市重點(diǎn)實(shí)驗(yàn)室，陜西 商洛 726000 3.河北省氣象科學(xué)研究所，河北 石家莊 050021）

商洛地區(qū)受氣候影響災(zāi)害性天氣頻發(fā)，由于該地區(qū)地形復(fù)雜，布設(shè)氣象站點(diǎn)困難，觀測(cè)實(shí)況資料相對(duì)缺乏，不能準(zhǔn)確監(jiān)測(cè)降水的空間分布和強(qiáng)度變化，給該地區(qū)的預(yù)報(bào)預(yù)警研究帶來(lái)困難。部署大量雨量計(jì)成本高、系統(tǒng)維護(hù)困難，數(shù)據(jù)質(zhì)量難以保證，因此引入天氣雷達(dá)以獲取較為準(zhǔn)確的降雨分布。天氣雷達(dá)可以獲得較高時(shí)空分辨率的降水分布，但是它探測(cè)面積有限，受下墊面和地理位置的影響嚴(yán)重，例如要處理波束遮擋、地物雜波、天氣垂向結(jié)構(gòu)變異性等諸多潛在觀測(cè)誤差。為了彌補(bǔ)臺(tái)站及天氣雷達(dá)對(duì)降水觀測(cè)的局限性和設(shè)備缺陷，填補(bǔ)資料缺乏的空白區(qū)域，引進(jìn)高空間分辨率的熱帶降水測(cè)量衛(wèi)星TRMM（Tropical Rainfall Measuring Mission）測(cè)雨雷達(dá)PR，可以在一定程度上彌補(bǔ)兩種觀測(cè)方式的不足。

熱帶降雨測(cè)量計(jì)劃TRMM 衛(wèi)星[1]，是1997 年由美國(guó)國(guó)家宇航局（NASA）與日本國(guó)家空間發(fā)展局（NASDA）聯(lián)合研制發(fā)射，2001 年8 月TRMM 衛(wèi)星軌道高度由350 km 抬升至402.5 km ，可覆蓋全球50°N～50°S 的范圍。

國(guó)內(nèi)外學(xué)者在不同地域?qū)RMM 衛(wèi)星降水的適用性和精度做了大量統(tǒng)計(jì)分析，研究表明TRMM衛(wèi)星降水產(chǎn)品在天山、浙江及京津冀地區(qū)均有良好表現(xiàn)，與地面氣象站點(diǎn)觀測(cè)降水相關(guān)性較高，隨著時(shí)空尺度的增大，其精度有一定提高[2-4]，在流域精度評(píng)價(jià)與結(jié)果分析中，TRMM 降水產(chǎn)品能更好地反映流域真實(shí)降水信息，作為其他降水產(chǎn)品的有效補(bǔ)充[5-7]。衛(wèi)星具有探測(cè)范圍廣，時(shí)空分辨率較高的優(yōu)勢(shì)，因此TRMM 衛(wèi)星降水資料在干旱特征分布、水文研究、暴雨研究、預(yù)報(bào)模式的數(shù)據(jù)同化、熱帶氣旋等方面得到較好應(yīng)用。如利用TRMM 衛(wèi)星高時(shí)空分辨率降水對(duì)區(qū)域尺度的干旱進(jìn)行監(jiān)測(cè)和特征分析[8-9]，唐國(guó)強(qiáng)等[10]通過(guò)TRMM 衛(wèi)星產(chǎn)品評(píng)估，證明TRMM 衛(wèi)星數(shù)據(jù)在贛江流域具有替代地面站點(diǎn)觀測(cè)的潛力。盧怡等[11]利用TRMM 的微波成像儀（TMI）觀測(cè)到的亮溫資料分析了亮溫與熱帶氣旋強(qiáng)度之間的關(guān)系。傅云飛等[12]基于TRMM 衛(wèi)星探測(cè)資料對(duì)東亞降水云結(jié)構(gòu)特征進(jìn)行了研究。Cesar 等[13]對(duì)TRMM 版本質(zhì)量在巴西南部進(jìn)行了評(píng)估，結(jié)果表明TRMM 極大地提高了對(duì)極端降水事件的定量識(shí)別，Savtchenko 等[14]利用17 a TMPA降水產(chǎn)品分析了在美國(guó)加利福尼亞州地區(qū)的降水變異。商洛地區(qū)山地多平地少，地勢(shì)起伏大，氣象站點(diǎn)分布稀少，很難直接通過(guò)插值手段獲取高時(shí)空分辨率的柵格化降水?dāng)?shù)據(jù)。因此引入基于TRMM 衛(wèi)星的遙感降水觀測(cè)數(shù)據(jù)彌補(bǔ)和訂正地面觀測(cè)站點(diǎn)的不足，以此提高商洛市降水的觀測(cè)質(zhì)量。

TRMM 降水?dāng)?shù)據(jù)精度及分布規(guī)律因研究區(qū)域所處地形、氣候條件等不同而有所差異，本文以商洛地區(qū)為研究區(qū)域，通過(guò)對(duì)比分析TRMM 衛(wèi)星降水?dāng)?shù)據(jù)與商洛地區(qū)7 個(gè)國(guó)家站點(diǎn)的實(shí)測(cè)降水?dāng)?shù)據(jù)，評(píng)估衛(wèi)星降水?dāng)?shù)據(jù)在商洛地區(qū)的精度及其與季節(jié)、地形等相關(guān)因子的關(guān)系，分析TRMM 降水?dāng)?shù)據(jù)在該區(qū)域的探測(cè)誤差，并采用比值系數(shù)法對(duì)TRMM 月降水?dāng)?shù)據(jù)產(chǎn)品進(jìn)行訂正，為TRMM 衛(wèi)星降水產(chǎn)品在商洛地區(qū)的研究與應(yīng)用提供可信依據(jù)。

1 研究區(qū)概況

商洛位于陜西東南部，地跨長(zhǎng)江、黃河兩大流域，總體上西北部海拔高，東南部海拔低。位于暖溫帶和北亞熱帶過(guò)渡地帶，受冬夏季風(fēng)和青藏高原環(huán)流的影響，商洛氣候呈現(xiàn)四季分明，雨熱同季，冬干夏濕，干濕分明的氣候特征，干旱、連陰雨、暴雨、冰雹、霜凍等災(zāi)害性天氣時(shí)有發(fā)生。商洛自動(dòng)站分布基本覆蓋鎮(zhèn)辦及大部分自然村（圖1），站點(diǎn)大部分部署在山區(qū)溝壑之中自然條件較好區(qū)域，平均間距在10 km 以上，嚴(yán)重影響對(duì)災(zāi)害性天氣的精細(xì)結(jié)構(gòu)探測(cè)能力。

圖1 商洛地形及國(guó)家地面站點(diǎn)分布

2 數(shù)據(jù)與方法

2.1 數(shù)據(jù)預(yù)處理

首先從降水測(cè)量計(jì)劃NASA 網(wǎng)站（https：//pmm.nasa.gov/）下載TRMM 衛(wèi)星降水?dāng)?shù)據(jù)產(chǎn)品資料，時(shí)間為2006—2015 年，時(shí)間分辨率為月，空間分辨率為0.25°×0.25°，版本為TRMM 3B43 V7，是由TRMM 3B42 數(shù)據(jù)產(chǎn)品、NOAA 氣候異常監(jiān)測(cè)系統(tǒng)（CAMS）的全球格點(diǎn)雨量測(cè)量器資料、全球降水氣候中心（GPCC）的全球降水資料聯(lián)合反演的降水產(chǎn)品。TRMM 數(shù)據(jù)是衛(wèi)星觀測(cè)月平均降水量產(chǎn)品，但數(shù)據(jù)記錄是每月的平均小時(shí)降水量（mm/h），本文通過(guò)Matlab 編程解析HDF 格式的降水資料，從中提取TRMM 降水速率層，利用經(jīng)緯度范圍從全球格網(wǎng)數(shù)據(jù)中裁剪出商洛地區(qū)的降水?dāng)?shù)據(jù)，通過(guò)小時(shí)累積計(jì)算得到研究時(shí)段內(nèi)各月區(qū)域格點(diǎn)平均降水量，對(duì)地面氣象站點(diǎn)和TRMM 衛(wèi)星格點(diǎn)數(shù)據(jù)按經(jīng)緯度進(jìn)行時(shí)空匹配[15]，得到商洛地區(qū)各氣象站點(diǎn)所在格網(wǎng)的TRMM 月降水?dāng)?shù)據(jù)，國(guó)家站點(diǎn)月平均降水量由日降水量累積計(jì)算得到。

地面氣象站點(diǎn)日降水量數(shù)據(jù)來(lái)源于全國(guó)綜合氣象信息共享平臺(tái)（CIMISS）數(shù)據(jù)環(huán)境，時(shí)間為2006—2015 年，根據(jù)數(shù)據(jù)結(jié)構(gòu)說(shuō)明文檔，批量將站點(diǎn)日降水?dāng)?shù)據(jù)資料統(tǒng)計(jì)得到月降水?dāng)?shù)據(jù)。由于區(qū)域自動(dòng)站沒(méi)有經(jīng)過(guò)質(zhì)控，經(jīng)檢查數(shù)據(jù)質(zhì)量較差，為準(zhǔn)確評(píng)估TRMM 衛(wèi)星降水量的精度，以商洛地區(qū)7 個(gè)國(guó)家地面站點(diǎn)實(shí)測(cè)降水?dāng)?shù)據(jù)作為真值來(lái)檢驗(yàn)TRMM 數(shù)據(jù)在商洛地區(qū)的有效性，通過(guò)建立TRMM 衛(wèi)星與國(guó)家地面站的關(guān)系進(jìn)一步訂正區(qū)域自動(dòng)站。

2.2 方法與指標(biāo)

本文采用統(tǒng)計(jì)分析指標(biāo)[15]：決定系數(shù)（R2）、相關(guān)系數(shù)（R）、相對(duì)偏差（Bias）和均方根誤差（RMSE）來(lái)評(píng)價(jià)TRMM 降水?dāng)?shù)據(jù)的精度。決定系數(shù)（Goodness of Fit）表示TRMM 衛(wèi)星降水回歸直線對(duì)觀測(cè)值的擬合程度，決定系數(shù)R2的值越接近1，TRMM 降水?dāng)?shù)據(jù)回歸直線對(duì)站點(diǎn)觀測(cè)值的擬合程度越好；相關(guān)系數(shù)R越大，TRMM 降水?dāng)?shù)據(jù)與站點(diǎn)降水?dāng)?shù)據(jù)的相關(guān)程度越高。相對(duì)偏差Bias用來(lái)衡量TRMM 衛(wèi)星觀測(cè)結(jié)果與雨量站點(diǎn)實(shí)測(cè)值之間的偏離程度，Bias越小準(zhǔn)確度越高；平均絕對(duì)偏差MAE和均方根誤差RMSE表示TRMM 降水?dāng)?shù)據(jù)與站點(diǎn)降水?dāng)?shù)據(jù)在數(shù)值上的差異和離散程度。其中相對(duì)偏差（Bias）的正負(fù)值表示TRMM 數(shù)據(jù)高估（低估）站點(diǎn)降水。

其中，Gi表示地面觀測(cè)站點(diǎn)的觀測(cè)值，Ti表示TRMM衛(wèi)星降水表示真實(shí)觀測(cè)值的平均值表示衛(wèi)星降水?dāng)M合值，BIAS為相對(duì)偏差，RMSE為均方根誤差。

分析TRMM 降水?dāng)?shù)據(jù)與站點(diǎn)降水?dāng)?shù)據(jù)評(píng)估結(jié)果，采用比值系數(shù)法來(lái)訂正TRMM 降水格點(diǎn)數(shù)據(jù)，以國(guó)家站點(diǎn)數(shù)據(jù)（經(jīng)過(guò)嚴(yán)格的質(zhì)量控制）作為真值，計(jì)算站點(diǎn)與對(duì)應(yīng)TRMM 產(chǎn)品對(duì)應(yīng)柵格的關(guān)系[16]；利用兩者之間的關(guān)系作為訂正因子，并將訂正因子空間化；利用空間化的訂正因子與原始TRMM 數(shù)據(jù)對(duì)之前的關(guān)系進(jìn)行反推，得到訂正后的TRMM 數(shù)據(jù)。比值系數(shù)訂正法公式如下：

3 結(jié)果分析

以商洛市國(guó)家地面氣象站點(diǎn)為代表性個(gè)例，分析2006—2015 年的TRMM 降水?dāng)?shù)據(jù)在商洛地區(qū)的精度及訂正效果。

3.1 時(shí)間尺度TRMM 衛(wèi)星降水精度

3.1.1 月尺度

由2006—2015 年研究區(qū)域內(nèi)地面站點(diǎn)與TRMM 衛(wèi)星月平均降水量的統(tǒng)計(jì)（圖2）可知，TRMM 月降水量數(shù)據(jù)與實(shí)測(cè)站點(diǎn)降水量數(shù)據(jù)在時(shí)間序列趨勢(shì)保持一致，呈明顯的周期性特征。研究區(qū)內(nèi)降水量主要集中在夏季，尤其是7、8 月，月平均降水量達(dá)150 mm，而冬季降水量總體偏小。整體而言，站點(diǎn)實(shí)測(cè)降水量較大時(shí)，TRMM 降水量產(chǎn)品往往會(huì)低估降水量。這是由于TRMM 衛(wèi)星測(cè)雨雷達(dá)PR 能源和重量的限制，無(wú)法提供大范圍降水測(cè)量，另外PR發(fā)射電磁波波長(zhǎng)大于地面雷達(dá)波長(zhǎng)，波長(zhǎng)越長(zhǎng)對(duì)降雨的感知衰減程度越小，傳感器的分辨率變大敏感度隨之降低，因此衛(wèi)星感知不到＜0.7 mm/h 的降雨，從而導(dǎo)致衛(wèi)星觀測(cè)降水結(jié)果相對(duì)偏小。當(dāng)?shù)孛嬲军c(diǎn)實(shí)測(cè)降水量為0 時(shí)，衛(wèi)星觀測(cè)往往顯示存在降雨，結(jié)果通常不為0，這是由于TRMM 衛(wèi)星降水是聯(lián)合反演均降水量，當(dāng)氣象站點(diǎn)實(shí)測(cè)降水為0 時(shí)，網(wǎng)格內(nèi)其他區(qū)域可能有降水。商洛地區(qū)TRMM 月降水量產(chǎn)品存在高估站點(diǎn)降水量的現(xiàn)象，可能與研究區(qū)的山地高原地形有關(guān)。通常地面雨量站點(diǎn)多部署在地勢(shì)平坦的地區(qū)，地形越復(fù)雜地區(qū)，探測(cè)受到降水垂直結(jié)構(gòu)、云型等多重因素干擾，衛(wèi)星探測(cè)降水量誤差往往越大。另外，TRMM 衛(wèi)星觀測(cè)降水產(chǎn)品空間分辨率為0.25°，分辨率比較低，在地形比較復(fù)雜的山區(qū)，雨量觀測(cè)站點(diǎn)降水與TRMM 衛(wèi)星觀測(cè)降水?dāng)?shù)據(jù)匹配不準(zhǔn)。

圖2 2006—2015 年各月平均降水量

圖3 為地面站點(diǎn)月降水量與TRMM 衛(wèi)星月降水量的散點(diǎn)斜率圖，其中2006—2015 年7 個(gè)國(guó)家站點(diǎn)的有效擬合點(diǎn)數(shù)為840 個(gè)。由圖3 可知，擬合結(jié)果的決定系數(shù)相對(duì)較高，R2為0.805，表明TRMM 月降水量數(shù)據(jù)與地面實(shí)測(cè)站點(diǎn)基本吻合，時(shí)間序列可以較好地描述地面降水的趨勢(shì)；偏差Bias為0.005，＞0，說(shuō)明TRMM 降水量結(jié)果值略高估站點(diǎn)實(shí)測(cè)值，TRMM 衛(wèi)星觀測(cè)結(jié)果與雨量站點(diǎn)實(shí)測(cè)值之間的偏離程度較小。RMSE和MAE相對(duì)較大，說(shuō)明TRMM 月降水量數(shù)據(jù)與站點(diǎn)實(shí)測(cè)降水仍存在較大的偏差。

圖3 站點(diǎn)平均月降水量和TRMM 衛(wèi)星降水量的擬合結(jié)果

3.1.2 季尺度

在月降水量分析的基礎(chǔ)上，對(duì)各個(gè)月的結(jié)果按季度統(tǒng)計(jì)，分為春、夏、秋、冬4 個(gè)季度，春季為3—5月，夏季為6—8 月，秋季為9—11 月，冬季為12月—次年2 月。

圖4 為4 個(gè)不同季度雨量站點(diǎn)降水量資料與TRMM 衛(wèi)星降水量資料的擬合結(jié)果。由圖4 可知，不同季度擬合結(jié)果不同，其中秋季擬合結(jié)果的R2最大，表明TRMM 衛(wèi)星在秋季探測(cè)降水的能力最強(qiáng)，春季和夏季的R2結(jié)果也相對(duì)較大，而冬季的R2則明顯較低，R2為0.75，說(shuō)明降水量較小時(shí)TRMM 衛(wèi)星觀測(cè)與實(shí)測(cè)站點(diǎn)的相關(guān)程度較低。R2表示的是兩組數(shù)據(jù)之間的密切程度，反映的是趨勢(shì)上的關(guān)系，并不能準(zhǔn)確的反映數(shù)據(jù)之間的偏差情況。冬季站點(diǎn)降水量與TRMM 衛(wèi)星降水量的決定系數(shù)最低，但是由于冬季降水量整體很小，因此冬季TRMM 降水量的均方根誤差和絕對(duì)誤差相較其他季節(jié)較小。

圖4 TRMM 衛(wèi)星與站點(diǎn)季降水?dāng)?shù)據(jù)散點(diǎn)（a 為春季，b 為夏季，c 為秋季，d 為冬季）

3.2 誤差指標(biāo)

表1 為國(guó)家站降水量與TRMM 月降水量數(shù)據(jù)擬合的統(tǒng)計(jì)結(jié)果。擬合優(yōu)度R2均在0.8 左右，說(shuō)明國(guó)家站與TRMM 衛(wèi)星擬合程度較好，R2值最大的站為商州，其站點(diǎn)年平均降水量為734.84 mm，而R2值最小的站為商南和柞水，其年均降水量分別為842.39 和794.18 mm。站點(diǎn)降水量與TRMM 衛(wèi)星觀測(cè)值存在較小的偏差Bias，從Bias的數(shù)值可以看出，商州、洛南和鎮(zhèn)安3 個(gè)站為負(fù)值外，其余站均為正值，表明商州、洛南和鎮(zhèn)安3 個(gè)站上TRMM 衛(wèi)星觀測(cè)的降水結(jié)果低估了實(shí)際的降水量，而其余站均高估了實(shí)際的降水量，總體呈現(xiàn)高估實(shí)際降水量的現(xiàn)象。RMSE和MAE在數(shù)值上基本保持一致，最大值出現(xiàn)在商南站，此站位于陜西省東南部，地處中緯度偏南地帶，北部和西北部屬暖溫帶氣候區(qū)，中部和南部屬于北亞熱帶氣候區(qū)，由于北有秦嶺天然屏障，阻擋寒潮不易侵入，雨量充沛，受地形影響較大。柞水和鎮(zhèn)安位于秦嶺南麓，商洛市西南，受地形影響誤差也較大，地形復(fù)雜地區(qū)站點(diǎn)的代表性較差。TRMM降水量數(shù)據(jù)一個(gè)網(wǎng)格的大小約為27.5 km×27.5 km，而一個(gè)雨量站往往只能表示其站點(diǎn)2 km 范圍內(nèi)的降水情況，站點(diǎn)匹配時(shí)也存在偏差。

表1 國(guó)家站各雨量站點(diǎn)擬合結(jié)果統(tǒng)計(jì)

3.3 TRMM 降水?dāng)?shù)據(jù)訂正

衛(wèi)星反演降水量數(shù)據(jù)具有覆蓋范圍廣、時(shí)間序列長(zhǎng)等優(yōu)勢(shì)，能夠在一定程度上彌補(bǔ)雨量站點(diǎn)不足，但其精度還需進(jìn)一步訂正。針對(duì)誤差指標(biāo)分析結(jié)果，利用2006—2015 年商洛各國(guó)家站降水?dāng)?shù)據(jù)對(duì)研究區(qū)內(nèi)TRMM 降水產(chǎn)品分別進(jìn)行了比值系數(shù)訂正，對(duì)訂正結(jié)果進(jìn)行誤差對(duì)比分析，結(jié)果見(jiàn)表2。利用比值系數(shù)訂正得到的TRMM 降水?dāng)?shù)據(jù)精度有很大提高，訂正后的TRMM 與氣象站點(diǎn)的決定系數(shù)R2均接近于1，平均相對(duì)誤差、偏差及均方根誤差均在一定程度上減小，分析比較表明，利用氣象站點(diǎn)對(duì)TRMM數(shù)據(jù)進(jìn)行訂正的方法具有一定的可靠性。

表2 國(guó)家站點(diǎn)比值系數(shù)訂正前后誤差比較

4 結(jié)論

對(duì)比商洛地區(qū)國(guó)家氣象站2006—2015 年實(shí)測(cè)月降水?dāng)?shù)據(jù)與氣象站所在格點(diǎn)TRMM 月降水?dāng)?shù)據(jù)，從月和季時(shí)間尺度上對(duì)TRMM 降水?dāng)?shù)據(jù)精度進(jìn)行了評(píng)估。采用決定系數(shù)、相對(duì)誤差等統(tǒng)計(jì)分析方法，研究了精度評(píng)價(jià)指標(biāo)的時(shí)間分布特征，同時(shí)分析了地形、降水強(qiáng)度等因素對(duì)精度的影響，對(duì)TRMM 遙感降水量產(chǎn)品的時(shí)空適用性和站點(diǎn)誤差進(jìn)行了探討。利用比值系數(shù)法對(duì)TRMM 月降水量進(jìn)行了訂正。主要結(jié)論如下：

（1）TRMM 月降水量數(shù)據(jù)與實(shí)測(cè)站點(diǎn)數(shù)據(jù)相關(guān)程度較高，國(guó)家站決定系數(shù)R2＞0.8，表明國(guó)家站降水和TRMM 月降水量有較好的擬合性；對(duì)月、季和年尺度的降水量結(jié)果進(jìn)行比較，季降水量結(jié)果的擬合結(jié)果較月降水量有所提高，TRMM 衛(wèi)星觀測(cè)降水與站點(diǎn)實(shí)測(cè)降水吻合度隨著時(shí)間累積有所越高。

（2）研究區(qū)內(nèi)TRMM 降水量產(chǎn)品存在明顯低估實(shí)際降水量的現(xiàn)象，利用比值系數(shù)法對(duì)TRMM 降水量產(chǎn)品進(jìn)行偏差訂正。驗(yàn)證結(jié)果表明，訂正后的TRMM 月降水結(jié)果偏差更小，決定系數(shù)更高。

（3）通過(guò)研究TRMM 降水產(chǎn)品在商洛的精度，選擇與氣象站相關(guān)性較好、均方根誤差較小的山區(qū)，利用實(shí)測(cè)降水資料與TRMM 降水產(chǎn)品建立關(guān)系模型，利用校準(zhǔn)系數(shù)對(duì)該區(qū)域內(nèi)所有的TRMM 格點(diǎn)資料進(jìn)行訂正，最終求得該區(qū)域平均面雨量。

綜上所述，TRMM 衛(wèi)星降水產(chǎn)品數(shù)據(jù)在商洛地區(qū)的精度較高，適用性較強(qiáng)。通過(guò)實(shí)測(cè)站對(duì)TRMM 數(shù)據(jù)進(jìn)行精度誤差訂正，為其在商洛地區(qū)的水文模擬預(yù)報(bào)、農(nóng)業(yè)及氣候研究的應(yīng)用提供參考價(jià)值。