黎揚(yáng)兵，張洪波,2，楊天增，呂豐光，王雨巍，姚聰聰

基于MGWR的渭河流域TRMM降水產(chǎn)品空間降尺度分析

黎揚(yáng)兵1，張洪波1,2※，楊天增1，呂豐光1，王雨巍1，姚聰聰1

（1. 長(zhǎng)安大學(xué)水利與環(huán)境學(xué)院，西安 710054； 2. 水利部旱區(qū)生態(tài)水文與水安全重點(diǎn)實(shí)驗(yàn)室，西安 710054）

高分辨率降水?dāng)?shù)據(jù)有助于刻畫降水的時(shí)空分異特性，對(duì)流域水文、氣象和生態(tài)等過程的精準(zhǔn)模擬具有重要作用，因此對(duì)低分辨率降水產(chǎn)品開展空間降尺度，提高其分辨率十分必要。鑒于此，該研究在充分考慮熱帶降雨測(cè)量衛(wèi)星（Tropical Rainfall Measuring Mission，TRMM）降水產(chǎn)品在渭河流域適用性的基礎(chǔ)上，引入歸一化差分植被指數(shù)（Normalized Difference Vegetation Index，NDVI）、數(shù)字高程模型（Digital Elevation Model，DEM）、坡度、坡向和經(jīng)緯度等地理環(huán)境因子，構(gòu)建了多尺度地理加權(quán)回歸（Multi-scale Geographically Weighted Regression，MGWR）模型用以分析不同因子對(duì)渭河流域降水空間格局影響的尺度差異；并提出了一種針對(duì)TRMM降水產(chǎn)品的空間降尺度方法，通過精度評(píng)價(jià)驗(yàn)證降尺度結(jié)果的可靠性。結(jié)果表明：1）TRMM降水產(chǎn)品數(shù)據(jù)相較于站點(diǎn)實(shí)測(cè)數(shù)據(jù)存在一定精度誤差，年尺度上2=0.807，BIAS=2.909%，RMSE=83.477 mm，表現(xiàn)較好；季尺度上秋季2最高，為0.847，夏季RMSE最大，為62.393 mm，四季的BIAS均較低；月尺度2為0.456～0.815，BIAS的絕對(duì)值介于0～14%之間，多數(shù)月份為正值，RMSE值域范圍為3.507～39.342 mm，精度較好；總體而言，TRMM降水產(chǎn)品數(shù)據(jù)在年、季和月尺度上均表現(xiàn)出良好的整體適用性。2）不同因子在干濕年份對(duì)降水空間分異格局的影響呈現(xiàn)出不同的尺度特征，其中濕潤(rùn)年的DEM、NDVI、坡向和經(jīng)緯度對(duì)降水呈現(xiàn)局部影響，坡度影響具有全局性，而干旱年各因子均表現(xiàn)為局部影響。3）流域和站點(diǎn)尺度上，降尺度TRMM數(shù)據(jù)相較于降尺度前產(chǎn)品數(shù)據(jù)精度得到一定改善，流域尺度上，2和RMSE得到有限提高；站點(diǎn)尺度上，各站點(diǎn)統(tǒng)計(jì)指標(biāo)變化各異，但降尺度后統(tǒng)計(jì)指標(biāo)整體優(yōu)于降尺度前，并且由于時(shí)間尺度上的誤差累積，站點(diǎn)年尺度數(shù)據(jù)精度相比月尺度數(shù)據(jù)稍差。4）降尺度TRMM數(shù)據(jù)相比于降尺度前產(chǎn)品數(shù)據(jù)，空間分布更細(xì)膩，細(xì)節(jié)特征表現(xiàn)更好，且在年、月時(shí)間尺度上均具有較高的精度，可為渭河流域資料短缺地區(qū)的水文設(shè)計(jì)提供數(shù)據(jù)支撐。

降水；TRMM；降尺度；MGWR模型；尺度差異

0 引 言

降水是驅(qū)動(dòng)地表物質(zhì)循環(huán)、能量交換以及影響水循環(huán)過程的關(guān)鍵因素，同時(shí)也是衡量區(qū)域及至全球氣候、生態(tài)變化及水文循環(huán)的指示因子，已被廣泛用于區(qū)域水量平衡計(jì)算、水資源評(píng)價(jià)、水文模擬和氣候應(yīng)對(duì)等研究[1-3]。目前，降水?dāng)?shù)據(jù)的采集方法主要包括站點(diǎn)觀測(cè)、衛(wèi)星反演降水和雷達(dá)測(cè)雨等?；谡军c(diǎn)觀測(cè)的傳統(tǒng)降水測(cè)量一般可得到較為準(zhǔn)確的降水?dāng)?shù)值，但受站點(diǎn)覆蓋能力和代表性所限，很難獲得區(qū)域高精度的降水空間分布信息，所以多采用空間插值方法進(jìn)一步獲取降水的空間分布格局[4]；地面測(cè)雨雷達(dá)雖然在一定程度上減少了稀少的地面站點(diǎn)插值帶來的誤差，并具備大面積遙測(cè)的優(yōu)點(diǎn)，但其覆蓋范圍一般較小，受實(shí)際環(huán)境影響也較大，尚缺乏普遍適用性[5]；衛(wèi)星遙感技術(shù)的迅猛發(fā)展，使降雨產(chǎn)品具備更高的時(shí)空分辨率、更廣的覆蓋范圍以及相對(duì)全面的適用環(huán)境等優(yōu)點(diǎn)，現(xiàn)已成為獲取空間化降水?dāng)?shù)據(jù)的重要方法[6]。

近年來，在全球大尺度和區(qū)域小尺度上的衛(wèi)星遙感降水產(chǎn)品方興未艾，其中熱帶降雨測(cè)量衛(wèi)星（Tropical Rainfall Measuring Mission，TRMM）產(chǎn)品憑借其靈活的時(shí)效性、較大的覆蓋范圍、可融合全球地面觀測(cè)站點(diǎn)及多源衛(wèi)星數(shù)據(jù)等優(yōu)勢(shì)[7]，在諸多衛(wèi)星降水產(chǎn)品中異軍突起[8]。但遺憾的是，目前TRMM降水產(chǎn)品的空間分辨率為0.25°（約27.5 km），對(duì)于區(qū)域的精細(xì)化研究仍顯粗糙，很難準(zhǔn)確刻畫小尺度流域和復(fù)雜流域的降水分布規(guī)律。因此，開展TRMM衛(wèi)星降水產(chǎn)品的空間降尺度研究，滿足新時(shí)期智慧水利的精細(xì)化需求，推進(jìn)區(qū)域降水時(shí)空變異性特征及相關(guān)研究尤為必要。

降尺度研究可以實(shí)現(xiàn)數(shù)據(jù)信息由大尺度向小尺度、低空間分辨率向高空間分辨率的轉(zhuǎn)化，其主要方法可概括為動(dòng)力降尺度、統(tǒng)計(jì)降尺度以及動(dòng)力統(tǒng)計(jì)降尺度相結(jié)合3類[9]。統(tǒng)計(jì)降尺度法因具有便于操作、計(jì)算量相對(duì)較小、方法眾多及應(yīng)用靈活等特點(diǎn)，已成為國(guó)內(nèi)外廣泛應(yīng)用且發(fā)展相對(duì)成熟的一種方法[10]。Immerzeel等[11]通過分析TRMM 3B43降水?dāng)?shù)據(jù)和年尺度歸一化植被指數(shù)（Normalized Difference Vegetation Index，NDVI）之間的相關(guān)關(guān)系，首次建立區(qū)域降水降尺度模型，得到1 km空間分辨率的伊比利半島TRMM降水?dāng)?shù)據(jù)。Jia等[12]進(jìn)一步考慮降水的空間分布不只受到單一變量因素影響，在Immerzeel等的研究基礎(chǔ)上，將數(shù)字高程模型（Digital Elevation Model，DEM）納入變量考慮范圍，建立了TRMM 3B43與DEM、NDVI之間的多元線性回歸（Multiple Linear Regression，MLR）模型，將柴達(dá)木盆地0.25°分辨率降水?dāng)?shù)據(jù)提升至1 km。同時(shí)，鑒于降水空間分布受多種因素協(xié)同作用，不少學(xué)者考慮采用各種典型機(jī)器學(xué)習(xí)方法來構(gòu)建降水降尺度模型，以期在概念混雜和物理機(jī)理不明確的情況下更合理地展示降水與影響因子之間的相關(guān)關(guān)系。例如，Shi等[13]以2001—2010年中國(guó)大陸TRMM 3B43年降水?dāng)?shù)據(jù)為研究對(duì)象，采用隨機(jī)森林（Random Forest，RF）方法進(jìn)行降尺度研究，提升其空間分辨率至1 km，并發(fā)現(xiàn)非線性模型對(duì)降水的空間分布表征能力更強(qiáng)。杜方洲等[14]以NDVI和地理因子為變量，使用深度學(xué)習(xí)算法構(gòu)建了降尺度模型，實(shí)現(xiàn)了TRMM 3B43產(chǎn)品在東北地區(qū)的精度提升，并將數(shù)據(jù)外延至50°N以上地區(qū)（產(chǎn)品未覆蓋區(qū)）。

綜上分析，不難發(fā)現(xiàn)上述方法大多基于全局尺度建立關(guān)系，并未體現(xiàn)降水與變量因子之間的空間非平穩(wěn)性和尺度依賴，亦未考慮數(shù)據(jù)的局部特性[15]。鑒于降水分布在空間上存在顯著的異質(zhì)性，已有部分學(xué)者就降水與地表特征參數(shù)間的空間非平穩(wěn)性關(guān)系展開了研究。Chen等[16]應(yīng)用地理加權(quán)回歸（Geographically Weighted Regression，GWR）方法捕捉不同位置的空間變化，并進(jìn)一步檢測(cè)了降水空間非平穩(wěn)性的特質(zhì)，對(duì)華北地區(qū)TRMM 3B43數(shù)據(jù)進(jìn)行降尺度，得到了1 km分辨率的降水?dāng)?shù)據(jù)。Zeng等[17]在分析降水與NDVI關(guān)系的基礎(chǔ)上，建立了時(shí)空地理加權(quán)回歸（Geographically and Temporally Weighted Regression，GTWR）降尺度模型，得到降水豐富地區(qū)的高空間分辨率數(shù)據(jù)，并在廣東省進(jìn)行了評(píng)估與驗(yàn)證。上述各類GWR模型雖然刻畫了變量因子間的空間非平穩(wěn)關(guān)系，但不同因子對(duì)降水影響的作用尺度差異并未被充分考慮。在實(shí)際中，不同變量因子對(duì)降水空間分布的作用尺度可能存在差異，表現(xiàn)出不同的異質(zhì)性和尺度范圍[18]。對(duì)于上述關(guān)鍵問題，F(xiàn)otheringham等[19]提出了多尺度地理加權(quán)回歸（Multi-scale Geographically Weighted Regression，MGWR），其可通過對(duì)不同參數(shù)設(shè)置不同帶寬結(jié)構(gòu)，更好地呈現(xiàn)不同因子對(duì)降水作用機(jī)制的尺度差異。

目前，國(guó)內(nèi)關(guān)于TRMM數(shù)據(jù)的降尺度研究已非常豐富，不一而足，但在中緯度的半濕潤(rùn)、半干旱地區(qū)相對(duì)較少。而渭河流域位于濕潤(rùn)氣候向半干旱氣候的過渡區(qū)域，其復(fù)雜的地形地貌，給降水的空間變異帶來了較強(qiáng)的影響。同時(shí)，已有區(qū)域研究多基于站點(diǎn)降水?dāng)?shù)據(jù)，以點(diǎn)擴(kuò)面，缺少對(duì)流域衛(wèi)星降水產(chǎn)品的尺度影響驗(yàn)證和降尺度研究[20]，使得對(duì)于準(zhǔn)確刻畫渭河流域水資源形成機(jī)制及空間異質(zhì)性、有效識(shí)別渭河流域水旱災(zāi)害形成機(jī)理及風(fēng)險(xiǎn)區(qū)劃、科學(xué)認(rèn)識(shí)氣候變化條件下渭河流域生態(tài)系統(tǒng)演化模式的支撐尚有不足。故本文擬以渭河流域?yàn)檠芯繀^(qū)，在對(duì)TRMM降水產(chǎn)品進(jìn)行適用性評(píng)估的基礎(chǔ)上，以NDVI、DEM、坡度、坡向和經(jīng)緯度為協(xié)變量，利用MGWR模型構(gòu)建0.25°分辨率上TRMM降水產(chǎn)品與協(xié)變量之間的空間非平穩(wěn)性多尺度關(guān)系，從而剖析協(xié)變量對(duì)降水影響的尺度效應(yīng)。同時(shí)，依據(jù)MGWR模型原理與多尺度關(guān)系，開展TRMM產(chǎn)品向1 km空間分辨率的降尺度研究，并采用多個(gè)精度評(píng)價(jià)指標(biāo)，以站點(diǎn)實(shí)測(cè)數(shù)據(jù)為基準(zhǔn)，對(duì)尺度下延的產(chǎn)品數(shù)據(jù)進(jìn)行精度評(píng)估，從而形成一種適用于渭河流域的降尺度方法，為流域?qū)崪y(cè)數(shù)據(jù)匱乏地區(qū)的水文分析與設(shè)計(jì)提供更好的數(shù)據(jù)支撐。

1 研究區(qū)概況與數(shù)據(jù)源

1.1 研究區(qū)概況

渭河流域位于中國(guó)西北核心地帶，西邊始于鳥鼠山北側(cè)，東邊到潼關(guān)入黃河，北邊以白于山為界，南邊終止于秦嶺山區(qū)，全域橫跨陜西、甘肅和寧夏。作為黃河第一大支流，流域處于104°00′E～110°20′E，33°50′N～37°18′N之間，總面積達(dá)13.49×104km2，約占黃河流域總面積的18%（圖1）[21]。流域地勢(shì)總體表現(xiàn)為西高東低，最大高差達(dá)3 000 m以上，起伏較大；流域水系呈不對(duì)稱羽狀分布，兩岸支流眾多，主要的兩大支流為涇河和北洛河。

注：DEM為數(shù)字高程模型。

渭河流域?qū)儆诘湫偷拇箨懶约撅L(fēng)氣候，四季變化較為明顯，春暖干旱少雨水，夏熱多雨有伏旱，秋涼濕潤(rùn)日照足，冬冷干燥水稀少。流域降雨時(shí)空分布不均，且年內(nèi)、年際變化較大，多年平均降雨量為500～800 mm。降水的不均勻性使得流域水資源分布呈現(xiàn)出時(shí)空上的差異性變化[22]。

1.2 數(shù)據(jù)源及預(yù)處理

1.2.1 數(shù)據(jù)來源

本研究旨在基于MGWR模型探究一種全新的降尺度方法，目標(biāo)實(shí)現(xiàn)對(duì)數(shù)據(jù)覆蓋的時(shí)段并不敏感，同時(shí)充分考慮站點(diǎn)數(shù)據(jù)的全面性和有效性，進(jìn)而采用2003—2016年的衛(wèi)星遙感數(shù)據(jù)和氣象站點(diǎn)數(shù)據(jù)。其中，衛(wèi)星遙感數(shù)據(jù)包括來自TRMM 3B43（降水?dāng)?shù)據(jù)）、MOD13A3（NDVI）和SRTM3（DEM）的相關(guān)數(shù)據(jù)集，氣象站點(diǎn)數(shù)據(jù)包括位于渭河流域內(nèi)及其周邊的21個(gè)站點(diǎn)的經(jīng)緯度、海拔和逐日降水量。數(shù)據(jù)類型及來源見表1，氣象站點(diǎn)名稱及空間分布如圖1所示。

表1 數(shù)據(jù)類型及來源

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

下載得到的原始TRMM產(chǎn)品數(shù)據(jù)為多幅倒置的.HDF格式影像，利用MATLAB軟件編碼，對(duì)其進(jìn)行批量格式轉(zhuǎn)換（轉(zhuǎn)為.tif格式）、影像旋轉(zhuǎn)、投影變換和單位變換（變?yōu)楹撩祝┑忍幚?，得?003－2016年逐月TRMM產(chǎn)品柵格數(shù)據(jù)。

NDVI數(shù)據(jù)按照MODIS中國(guó)區(qū)域行列號(hào)進(jìn)行下載，為同一時(shí)期多幅.HDF格式影像。首先利用NASA提供的MRT（MODIS Reprojection Tool）軟件對(duì)其進(jìn)行波段提取、拼接和投影變換等處理，得到NDVI柵格數(shù)據(jù)，之后采用BAT轉(zhuǎn)換進(jìn)行批量操作，得到研究時(shí)段內(nèi)逐月柵格數(shù)據(jù)集。在此基礎(chǔ)上，使用Python進(jìn)行影像裁剪，無效值和異常像元剔除，然后在ArcGIS平臺(tái)上運(yùn)用IDW方法填充影像空值[23]，從而獲得更優(yōu)質(zhì)的1 km分辨率NDVI月值數(shù)據(jù)集；最后對(duì)逐月柵格數(shù)據(jù)集進(jìn)行最大值合成，得到2003－2016年1 km年尺度NDVI數(shù)據(jù)集。

將DEM數(shù)據(jù)在ArcGIS平臺(tái)上進(jìn)行投影變換和裁剪，再將其重采樣至1 km空間分辨率，并計(jì)算生成1 km分辨率的坡度和坡向數(shù)據(jù)。

氣象站點(diǎn)實(shí)測(cè)降水?dāng)?shù)據(jù)為逐日降水量。為保證數(shù)據(jù)質(zhì)量，對(duì)其進(jìn)行數(shù)據(jù)檢查，包括進(jìn)行無效值剔除和缺失值插補(bǔ)等處理。之后對(duì)逐日降水量進(jìn)行時(shí)間尺度上的累加，得到月和年尺度氣象站點(diǎn)的實(shí)測(cè)降水量。

2 研究方法

2.1 MGWR模型

MGWR是一種局部回歸模型，是在傳統(tǒng)GWR模型基礎(chǔ)上進(jìn)行改進(jìn)得到，優(yōu)點(diǎn)是可以消除GWR模型建模過程中變量因子均為同一帶寬的局限性，揭示變量因子作用尺度的差異性。Yu等[24]對(duì)該模型局部參數(shù)的統(tǒng)計(jì)推斷方法進(jìn)行了進(jìn)一步完善與改進(jìn)，現(xiàn)已被廣泛地應(yīng)用于不同學(xué)科領(lǐng)域的實(shí)證研究。

MGWR模型的表達(dá)式[19]如下：

MGWR模型在回歸過程中對(duì)所有協(xié)變量選擇差異性帶寬，使預(yù)設(shè)模型建模過程能在不同的空間尺度上運(yùn)行。模型采用后退擬合算法（Back-Fitting Algorithm）進(jìn)行模型校準(zhǔn)，以完成參數(shù)的估計(jì)迭代。主要做法是利用經(jīng)典GWR模型的估計(jì)值初始化擬合過程，以迭代的方式確定最優(yōu)帶寬和新的局部參數(shù)估計(jì)，然后不斷更新參數(shù)估計(jì)，當(dāng)差值達(dá)到收斂的指定閾值時(shí)，迭代結(jié)束，實(shí)現(xiàn)校準(zhǔn)。

本研究為實(shí)現(xiàn)TRMM數(shù)據(jù)空間尺度轉(zhuǎn)換，分析各變量因子的作用尺度差異，構(gòu)建了TRMM降水產(chǎn)品數(shù)據(jù)與NDVI、DEM、坡度和坡向之間的MGWR模型，其模擬計(jì)算過程均在MGWR 2.2軟件中進(jìn)行。主要是將0.25°分辨率的TRMM數(shù)據(jù)、NDVI、DEM、坡度和坡向以及經(jīng)緯度數(shù)據(jù)處理成Excel文件作為模型輸入，然后設(shè)置模型運(yùn)行參數(shù)，其中核函數(shù)和帶寬選擇準(zhǔn)則沿用經(jīng)典GWR設(shè)置（選擇常用的二次核函數(shù)（Bisquare）作為空間核函數(shù)類型，更正的赤池信息準(zhǔn)則（Akaike information criterion，AICc）為帶寬選擇準(zhǔn)則），同時(shí)設(shè)定最優(yōu)帶寬搜索類型為黃金分割算法（Golden Section），進(jìn)而實(shí)現(xiàn)模型模擬。

2.2 降水降尺度轉(zhuǎn)換

由于低分辨率遙感數(shù)據(jù)攜帶的信息量偏少，所以在降尺度轉(zhuǎn)換過程中需要引入額外信息，以此來實(shí)現(xiàn)尺度轉(zhuǎn)換[25]。而降水作為大氣運(yùn)動(dòng)過程的基本輸出通量，其強(qiáng)弱和分布受到諸多因素的影響[26]，所以結(jié)合前人的研究成果[27-29]，引入NDVI、DEM、坡度和坡向數(shù)據(jù)，并結(jié)合經(jīng)緯度作為TRMM數(shù)據(jù)降尺度模型的變量因子。其中，就地理因子而言，研究區(qū)域的經(jīng)度決定了海陸距離的遠(yuǎn)近，經(jīng)度越高，距離越小，降水發(fā)生概率越大，而緯度決定著區(qū)域的溫度和濕度，從而影響降水的發(fā)生[27]；從地形因子來看，地形起伏導(dǎo)致水汽對(duì)流活動(dòng)的發(fā)生，促使低層水汽上升并形成局地氣候，進(jìn)而產(chǎn)生豐富的降水[28]；植被與降水形成互饋，對(duì)降水?dāng)?shù)據(jù)的空間變異信息進(jìn)行補(bǔ)充，水分的多寡制約著植被的長(zhǎng)勢(shì)[29]。

考慮到降水量和植被、地理及地形等因素之間可能存在空間非平穩(wěn)性，且NDVI、高程、坡度和坡向?qū)邓诳臻g上的影響存在顯著的作用尺度差異，所以本研究選用MGWR模型來實(shí)現(xiàn)TRMM降水產(chǎn)品的空間尺度轉(zhuǎn)換。具體的降尺度轉(zhuǎn)換步驟如下：

1）前置處理。為使TRMM數(shù)據(jù)完全覆蓋流域邊界及選用的周邊臨近站點(diǎn)，并考慮地形因素對(duì)降水的影響，故沿著渭河流域邊界向外拓展0.5°的緩沖區(qū)作為降尺度研究區(qū)域。

2）數(shù)據(jù)準(zhǔn)備。將預(yù)處理得到的1 km分辨率的逐月NDVI、DEM、坡度和坡向數(shù)據(jù)進(jìn)行重采樣，使其達(dá)到0.25°空間分辨率。并將重采樣前后數(shù)據(jù)分別進(jìn)行分組，包括高分辨率降尺度數(shù)據(jù)集（1 km）和低分辨率建模數(shù)據(jù)集（0.25°）。

3）構(gòu)建MGWR模型。由于降水從降落到地表，然后被植被吸收并反映到植被蓋度的變化需要一定的時(shí)間，故植被NDVI對(duì)降水的響應(yīng)多表現(xiàn)為一定的滯后性[30]。相關(guān)研究成果已表明，其滯后期約為1個(gè)月[30-32]。因此，本研究選取滯后1個(gè)月的NDVI數(shù)據(jù)作為當(dāng)月的降尺度模型變量因子。以0.25°分辨率的TRMM數(shù)據(jù)為響應(yīng)變量，低分辨率建模數(shù)據(jù)集（0.25°）中滯后1個(gè)月的NDVI、DEM、坡度、坡向?yàn)閰f(xié)變量，建立MGWR模型。

4）模型結(jié)果處理。由式（1）知，經(jīng)過MGWR模型運(yùn)算，得到的輸出結(jié)果為常數(shù)項(xiàng)（經(jīng)緯度信息）、變量因子（NDVI、DEM、坡度和坡向）系數(shù)和殘差。對(duì)輸出結(jié)果進(jìn)行處理，主要步驟包括：先將上述輸出結(jié)果柵格化，對(duì)柵格化后的殘差結(jié)果進(jìn)行反距離權(quán)重插值得到1 km分辨率的殘差值，將柵格化后的各變量因子系數(shù)和常數(shù)項(xiàng)重采樣至1 km分辨率。

5）尺度轉(zhuǎn)換。根據(jù)MGWR模型原理，將高分辨降尺度數(shù)據(jù)集（1 km）中的數(shù)據(jù)和對(duì)應(yīng)的1 km分辨率的各變量因子系數(shù)、常數(shù)項(xiàng)和殘差回代到式（1）中，即可得到1 km空間分辨率降尺度TRMM數(shù)據(jù)。

重復(fù)上述步驟，得到2003－2016年逐月1 km空間分辨率降尺度數(shù)據(jù)集。

2.3 精度評(píng)價(jià)

以氣象站點(diǎn)實(shí)測(cè)數(shù)據(jù)為“真值”，引入擬合優(yōu)度（2，也稱決定系數(shù)）、相對(duì)偏差（BIAS）和均方根誤差（Root Mean Square Error，RMSE）3項(xiàng)指標(biāo)對(duì)TRMM降水產(chǎn)品數(shù)據(jù)和基于TRMM的降尺度模擬降水?dāng)?shù)據(jù)進(jìn)行適用性和精度分析，主要計(jì)算式[33]為

其中：

式中T為TRMM產(chǎn)品降水量或降尺度降水量，mm；P為氣象站點(diǎn)實(shí)測(cè)降水量，mm；T為TRMM產(chǎn)品降水量或降尺度降水量的平均值，mm；`為站點(diǎn)降水量的平均值，mm；為樣本數(shù)。

擬合優(yōu)度2用以反映TRMM產(chǎn)品降水量或降尺度降水量與氣象站點(diǎn)實(shí)測(cè)降水量之間的一致性，取值范圍為[0，1]，數(shù)值越逼近于1，表明2種數(shù)據(jù)的相關(guān)性越高；BIAS可以評(píng)估以上數(shù)據(jù)之間的偏離程度，其值越接近0，說明偏離越小，數(shù)據(jù)越精確；RMSE用來衡量數(shù)據(jù)誤差的整體水平，其值越小，表明兩者之間的誤差越小。

3 結(jié)果與分析

3.1 TRMM數(shù)據(jù)適用性分析

3.1.1 年尺度適用性

在年尺度上，選擇渭河流域21個(gè)氣象站點(diǎn)2003－2016年的年實(shí)測(cè)降水量作為自變量，站點(diǎn)對(duì)應(yīng)格點(diǎn)位置的年TRMM產(chǎn)品降水量作為因變量，進(jìn)行一元線性回歸分析，并計(jì)算2組數(shù)據(jù)之間的精度評(píng)價(jià)指標(biāo)，所得結(jié)果如圖2a所示。同時(shí)，將TRMM降水量和站點(diǎn)實(shí)測(cè)降水量處理為年均降水量，并繪制時(shí)間序列圖，如圖2b。由圖2a可知，年降水量的擬合優(yōu)度2=0.807，相對(duì)偏差BIAS=2.909%，均方根誤差RMSE=83.477 mm，整體來看，除RMSE值偏大外，其余評(píng)價(jià)指標(biāo)均表現(xiàn)良好，這表明TRMM產(chǎn)品數(shù)據(jù)與氣象站點(diǎn)實(shí)測(cè)數(shù)據(jù)的一致性較高。從圖2b可以看出，TRMM產(chǎn)品數(shù)據(jù)與站點(diǎn)實(shí)測(cè)數(shù)據(jù)年際變化趨勢(shì)總體一致，時(shí)序同步，但TRMM產(chǎn)品的年尺度數(shù)據(jù)對(duì)站點(diǎn)實(shí)測(cè)數(shù)據(jù)存在一定程度的高估現(xiàn)象，其年均降水比站點(diǎn)實(shí)測(cè)降水平均高16.68 mm。

注：R2為擬合優(yōu)度；BIAS為相對(duì)偏差；RMSE為均方根誤差。下同。

3.1.2 季尺度適用性

將研究區(qū)內(nèi)2003－2016年TRMM產(chǎn)品月尺度降水量和21個(gè)氣象站點(diǎn)月尺度實(shí)測(cè)降水量分別按春季（3－5月）、夏季（6－8月）、秋季（9－11月）、冬季（12－翌年2月）進(jìn)行統(tǒng)計(jì)，并對(duì)2組數(shù)據(jù)進(jìn)行一元線性回歸擬合及相關(guān)性分析，得到散點(diǎn)圖和相關(guān)指標(biāo)統(tǒng)計(jì)結(jié)果，如圖3所示。結(jié)果表明，秋季（2=0.847）擬合優(yōu)度最高，春季（2=0.781）和夏季（2=0.726）次之，冬季（2=0.584）最差。由于渭河流域降水主要集中于汛期，且常出現(xiàn)強(qiáng)降水，而TRMM衛(wèi)星雷達(dá)難于捕捉強(qiáng)降水信號(hào)，導(dǎo)致數(shù)據(jù)異常值較多，故夏季均方根誤差最大（RMSE= 62.393 mm）。從相對(duì)偏差結(jié)果來看，所有季節(jié)TRMM產(chǎn)品降水量都存在一定高估現(xiàn)象，但除春季略大外（5.19%），其余季節(jié)數(shù)據(jù)精度相對(duì)較高。綜合3個(gè)精度指標(biāo)的評(píng)價(jià)結(jié)果，認(rèn)為TRMM產(chǎn)品數(shù)據(jù)在渭河流域季尺度上表現(xiàn)出了相對(duì)較好的精度。

圖3 季尺度TRMM數(shù)據(jù)適用性評(píng)估

3.1.3 月尺度適用性

對(duì)渭河流域21個(gè)氣象站點(diǎn)的TRMM產(chǎn)品降水量和實(shí)測(cè)降水量在月尺度上進(jìn)行適用性分析，計(jì)算兩者的2、BIAS和RMSE，得到各月統(tǒng)計(jì)指標(biāo)值（如表2）。從表2可以看出，TRMM產(chǎn)品數(shù)據(jù)和站點(diǎn)實(shí)測(cè)月降水量總體相關(guān)性較好，2在0.456～0.815之間浮動(dòng)，其中最小的月份是6月，最大的月份為9月；從BIAS指標(biāo)值來看，BIAS的絕對(duì)值總體介于0～14%之間，2月最高（BIAS=13.061%），9月最低（BIAS=0.252%），另外除4月和8月為負(fù)值外，其余月份均為正值，表明TRMM產(chǎn)品降水量在絕大數(shù)情況下對(duì)站點(diǎn)降水存在一定程度的高估；RMSE指標(biāo)值域范圍為3.507～39.342 mm，較年、季尺度上數(shù)值偏低，說明月尺度上的TRMM產(chǎn)品數(shù)據(jù)和站點(diǎn)實(shí)測(cè)降水量之間的誤差相對(duì)偏小。綜合精度評(píng)估各類統(tǒng)計(jì)指標(biāo)值可發(fā)現(xiàn)，月尺度TRMM產(chǎn)品降水量與站點(diǎn)實(shí)測(cè)降水量之間有較好的一致性，說明其在渭河流域上具有較強(qiáng)的適用性。

表2 TRMM月數(shù)據(jù)適用性評(píng)估

綜上，認(rèn)為在時(shí)空尺度上，研究區(qū)域內(nèi)TRMM降水產(chǎn)品數(shù)據(jù)適用性和精度較好，采用TRMM降水產(chǎn)品數(shù)據(jù)在渭河流域進(jìn)行降尺度模擬研究具有一定的科學(xué)基礎(chǔ)和合理性。

3.2 基于MGWR模型的因子作用尺度分析

為了探究地形、地理和植被因子對(duì)降水空間分布的作用尺度差異，本文選取研究期內(nèi)的2003年作為典型濕潤(rùn)年，2008年作為典型干旱年，以納入空間因素的MGWR與GWR同時(shí)建?；貧w，從而得到各變量因子對(duì)降水空間格局產(chǎn)生影響時(shí)的作用尺度，相關(guān)參數(shù)如表3所示。

表3 多尺度地理加權(quán)回歸與地理加權(quán)回歸間的因子帶寬差異

注：NDVI為歸一化差分植被指數(shù)。

Note: NDVI is Normalized Difference Vegetation Index.

由MGWR模型原理可知，帶寬大小可以反映出作用尺度的差異，故可依據(jù)帶寬大小將變量因子分為局地影響因子和全局影響因子，而GWR模型只能表征每個(gè)協(xié)變量因子作用尺度的平均值。由MGWR建模因子參數(shù)，可發(fā)現(xiàn)濕潤(rùn)年（2003年）的變量因子中存在作用尺度較大的值（即坡度因子帶寬為357，接近于全局尺度），而干旱年（2008年）各變量尺度差異不明顯（表3）。就濕潤(rùn)年而言，地理位置、NDVI、DEM和坡向?qū)邓酗@著影響，其作用尺度較?。ň鶠?3），而坡度的影響較為平緩（帶寬為357），接近于全局尺度；對(duì)于干旱年來說，各變量因子均顯現(xiàn)出差異不大的作用尺度，整體上表現(xiàn)為局地影響。整體來看，渭河流域降水隨空間變化在不同緯度區(qū)域、地勢(shì)高度、植被覆蓋區(qū)及坡向間差異顯著。而坡度對(duì)降水的空間分異特征影響體現(xiàn)在任何尺度范圍內(nèi)，這主要是因?yàn)槿饔蚱露确旨?jí)中，平坡占據(jù)了流域總面積的三分之二，致使其對(duì)濕潤(rùn)年降水的影響呈現(xiàn)出了全局性。

為進(jìn)一步分析回歸系數(shù)空間格局，統(tǒng)計(jì)了MGWR模型在濕潤(rùn)年和干旱年各變量因子系數(shù)的描述性結(jié)果，將回歸結(jié)果進(jìn)行可視化，得到了系數(shù)空間分布，如圖4所示。

注：紅線圈畫區(qū)域?yàn)榛貧w系數(shù)通過0.05顯著性檢驗(yàn)的區(qū)域。

從圖4可以看出，在濕潤(rùn)年的回歸結(jié)果中，常數(shù)項(xiàng)的回歸系數(shù)整體上顯著，其系數(shù)取值介于605.686～790.843之間，呈現(xiàn)出較大的變動(dòng)范圍，遂可見常數(shù)項(xiàng)所表征的地理位置信息對(duì)降水呈現(xiàn)顯著的正向作用。系數(shù)高值區(qū)主要集中于渭河下游的南岸山區(qū)秦嶺北麓和涇河谷地地區(qū)，低值則主要分布在涇河和北洛河上游（黃土丘陵溝壑區(qū)）及渭河干流上游?；跐駶?rùn)年實(shí)際的植被覆蓋狀態(tài)，可知NDVI對(duì)降水的影響在涇河上游黃土丘陵溝壑區(qū)呈負(fù)向作用，且該區(qū)域回歸系數(shù)不顯著，其余地區(qū)表現(xiàn)為正向作用，回歸系數(shù)取值在-54.737～209.323之間，且系數(shù)高值集中趨勢(shì)與常數(shù)項(xiàng)大致相同。DEM和坡向?qū)邓挠绊懢尸F(xiàn)出由負(fù)向至正向過渡的趨勢(shì)，其中，DEM的回歸系數(shù)在關(guān)中平原、涇河谷地和北洛河谷地不顯著，系數(shù)取值在?0.111～0.042之間；坡向回歸系數(shù)顯著區(qū)主要位于涇河谷地、渭河干流中上部（六盤山分界區(qū)域）以及涇河上游黃土丘陵溝壑區(qū)，系數(shù)取值范圍為?0.348～0.149，且兩者回歸系數(shù)的正負(fù)高值均主要分布在系數(shù)顯著區(qū)。坡度在渭河流域表現(xiàn)為全局變量，其回歸系數(shù)在流域整體上并不顯著，系數(shù)取值范圍為?1.417～?0.259，呈西北向東南梯度遞增，均值為?0.791，坡度每抬升1%，降水值平均下降0.791%。

對(duì)于干旱年，從變量因子的回歸系數(shù)來看，除坡度外，其余因子的顯著區(qū)域都較濕潤(rùn)年少，表現(xiàn)出較小的絕對(duì)值。常數(shù)項(xiàng)和NDVI回歸系數(shù)在全域上均滿足顯著性要求，而坡度和坡向的回歸系數(shù)在渭河流域全域上顯著區(qū)極少。對(duì)于DEM，其對(duì)降水的影響整體上顯著，不顯著區(qū)域基本同濕潤(rùn)年一致，主要集中于關(guān)中平原、涇河谷地和北洛河谷地區(qū)域。其中，表征地理位置信息的常數(shù)項(xiàng)系數(shù)取值范圍介于347.814～501.623之間；代表植被覆蓋狀態(tài)的NDVI在全流域上對(duì)降水的影響表現(xiàn)為正向作用，NDVI值每上升1%，降水值抬升28.089%～243.776%，并且常數(shù)項(xiàng)和NDVI系數(shù)的空間分布趨勢(shì)具有較好的一致性。而DEM、坡度和坡向在空間上都表現(xiàn)出正負(fù)兩向作用，DEM回歸系數(shù)呈由西北至東南逐步由負(fù)過渡為正的趨勢(shì)，而坡度回歸系數(shù)正值主要位于流域東、西部?jī)蓚?cè)，負(fù)值主要位于中部河谷和DEM低值區(qū)，坡向回歸系數(shù)同坡度呈現(xiàn)相似的規(guī)律，僅有別于正值分布在流域西北和東南區(qū)。

從上述回歸結(jié)果來看，不論干濕年，常數(shù)項(xiàng)和NDVI的回歸系數(shù)都整體顯著，是反映降水空間分異的主要作用因素。對(duì)比所有變量因子回歸系數(shù)的均值絕對(duì)值，可以看出，常數(shù)項(xiàng)的影響程度在所有變量中最大，其表征的地理位置對(duì)降水的空間分異起著決定性的作用；NDVI的影響程度較大，反映出植被覆蓋狀態(tài)對(duì)降水分布格局的影響；而DEM回歸系數(shù)不顯著的區(qū)域集中于它的低值區(qū)，這是因?yàn)楹０紊邥?huì)對(duì)氣流產(chǎn)生膨脹和冷卻效應(yīng)，使其相對(duì)濕度增加，導(dǎo)致DEM高值區(qū)降水發(fā)生[34]。坡度和坡向回歸系數(shù)整體顯著性差，且坡度在干濕年顯示出了不同的尺度效應(yīng)，究其原因，2003年在大尺度環(huán)流背景下發(fā)生特大洪災(zāi)[35]，大氣的水汽輸送和垂直運(yùn)動(dòng)受到坡度的影響，減弱了空間降水的產(chǎn)生和分配[36]；而濕潤(rùn)年和干旱年水汽來源不一致，也會(huì)驅(qū)動(dòng)坡向上的降水量產(chǎn)生明顯差異[37]。

3.3 降尺度結(jié)果與精度評(píng)價(jià)

3.3.1 總體精度評(píng)價(jià)

基于21個(gè)氣象站點(diǎn)2003－2016年逐月降水量數(shù)據(jù)對(duì)TRMM降水產(chǎn)品降尺度前后數(shù)據(jù)進(jìn)行總體精度檢驗(yàn)，得到TRMM降尺度前后數(shù)據(jù)與所有站點(diǎn)逐月降水量的密度散點(diǎn)圖，一共匹配了3 528對(duì)數(shù)據(jù)，具體如圖5所示。

注：N為樣本數(shù)。 Note: N is sample size.

由圖5可以看出，TRMM數(shù)據(jù)降尺度前后精度相差不大，各站點(diǎn)月降水量分布在0～600 mm之間，且可以看到數(shù)據(jù)在0～50 mm區(qū)間聚集，密度較大。降尺度后的數(shù)據(jù)比原始產(chǎn)品數(shù)據(jù)的離散程度要略小。同時(shí)，從精度評(píng)價(jià)指標(biāo)的結(jié)果對(duì)比來看，降尺度數(shù)據(jù)的擬合優(yōu)度2、BIAS以及RMSE的值變化有限，但均略優(yōu)于TRMM原始產(chǎn)品數(shù)據(jù)。

3.3.2 站點(diǎn)精度評(píng)價(jià)

由以上分析可知，降尺度后的TRMM數(shù)據(jù)在流域尺度上具有較好的精度，但總體精度檢驗(yàn)難免會(huì)掩蓋單一站點(diǎn)降尺度數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)之間的差異。因此，在總體檢驗(yàn)的基礎(chǔ)上，對(duì)單個(gè)氣象站點(diǎn)降尺度結(jié)果進(jìn)行評(píng)估尤為必要。

計(jì)算21個(gè)氣象站點(diǎn)月、年尺度降尺度前后TRMM數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)之間的統(tǒng)計(jì)指標(biāo)值，得到統(tǒng)計(jì)指標(biāo)的空間分布，如圖6所示。由圖6可知：1）在月和年尺度上，站點(diǎn)降尺度后的統(tǒng)計(jì)指標(biāo)值相較于原始產(chǎn)品指標(biāo)值有一定改善，但原始產(chǎn)品數(shù)據(jù)本身精度較差的站點(diǎn)，其降尺度結(jié)果也較為不理想；2）不同時(shí)間尺度下，統(tǒng)計(jì)指標(biāo)空間分布存在明顯差異，且隨著時(shí)間尺度增加，TRMM降水?dāng)?shù)據(jù)誤差不斷累積，致使年尺度統(tǒng)計(jì)指標(biāo)值要差于月尺度統(tǒng)計(jì)指標(biāo)值。但總體而言，降尺度后的TRMM數(shù)據(jù)在各站點(diǎn)上確有較好的精度，在原始產(chǎn)品精度較好的站點(diǎn)（如崆峒，鳳翔和臨洮等）上，降尺度前后指標(biāo)值改善較大（2增加約6%；BIAS減少約4%，更趨近于0；RMSE減少約3mm），也體現(xiàn)了模型在流域上的適用性。

圖6 降水量降尺度前后統(tǒng)計(jì)指標(biāo)空間分布

3.3.3 年尺度TRMM降尺度結(jié)果分析

在年尺度上構(gòu)建MGWR模型，利用降尺度方法得到各年1 km分辨率TRMM降水量，并對(duì)降尺度前后數(shù)據(jù)求平均，得到多年平均值；同時(shí)，為避免降水在時(shí)間上存在不連續(xù)性和高度非線性特征[38]，同前文一致，選擇2003年作為濕潤(rùn)年、2008年作為干旱年進(jìn)行典型年份降水降尺度結(jié)果分析，從而得到年尺度TRMM降水空間分布情況，如圖7所示。

圖7 年尺度TRMM降水空間分布

由圖7可知：1）TRMM產(chǎn)品數(shù)據(jù)空間分辨率過低，一個(gè)柵格單元覆蓋面積較大，柵格內(nèi)部降水的空間異質(zhì)性無法得到體現(xiàn)。而降尺度后的TRMM數(shù)據(jù)空間分辨率得到很大提升，其空間分布細(xì)節(jié)展現(xiàn)更好。以流域部分氣象站點(diǎn)（西峰、麥積和鳳翔等）為例，可發(fā)現(xiàn)這些氣象站點(diǎn)降尺度前靠近柵格單元邊界，且主要以此像元值為降水值，無法體現(xiàn)臨近像元值的作用，而降尺度后數(shù)據(jù)融合臨近像元值，柵格單元面積更小，格點(diǎn)之間平滑過渡，能更細(xì)致地體現(xiàn)流域內(nèi)降水的空間分布特征；2）濕潤(rùn)年降水量空間分布趨勢(shì)與多年年均降水相近，總體上表現(xiàn)為西北部降水偏小、東南部較大，且自西北向東南逐漸增加的趨勢(shì)；而干旱年降水西北部偏小，降水空間分布在地理上未呈現(xiàn)出較明顯的變化趨勢(shì)。同時(shí)可看到，濕潤(rùn)年降水主要出現(xiàn)在渭河干流中下游、兩大支流的中游以及秦嶺北麓，而干旱年降水雨量較少，大的降水主要分布在渭河下游的南岸山區(qū)秦嶺北麓。

利用研究區(qū)域21個(gè)氣象站點(diǎn)提取濕潤(rùn)年和干旱年降尺度前后TRMM降水量，結(jié)合站點(diǎn)實(shí)測(cè)降水量值，計(jì)算得到不同統(tǒng)計(jì)指標(biāo)值，如表4所示。

表4 典型年TRMM降水量降尺度前后精度對(duì)比

由表4可知，干旱年和濕潤(rùn)年的降尺度結(jié)果與TRMM產(chǎn)品數(shù)據(jù)相比，擬合優(yōu)度2均有所提升，且相對(duì)偏差BIAS和均方根誤差RMSE亦存在一定幅度下降，這表明降尺度后TRMM數(shù)據(jù)較原始產(chǎn)品數(shù)據(jù)精度得到了一定的提高。同時(shí)研究發(fā)現(xiàn)，濕潤(rùn)年降尺度后的統(tǒng)計(jì)指標(biāo)2和BIAS值都優(yōu)于干旱年的對(duì)應(yīng)指標(biāo)值，只有均方根誤差略高，遂可見降尺度方法在濕潤(rùn)年的降尺度結(jié)果略優(yōu)于干旱年份。

3.3.4 月尺度TRMM降尺度精度評(píng)估

為了進(jìn)一步實(shí)現(xiàn)對(duì)降尺度前后降水量的精度檢驗(yàn)和對(duì)比分析，分別統(tǒng)計(jì)歷年降尺度前后TRMM各月降水量，計(jì)算其與站點(diǎn)實(shí)測(cè)月降水量的2、BIAS和RMSE，其精度評(píng)價(jià)指標(biāo)結(jié)果及對(duì)比分析如表5所示。由表可以看出，除1月外，經(jīng)過降尺度處理后的TRMM降水量與站點(diǎn)實(shí)測(cè)降水量的擬合優(yōu)度在各月份均有不同程度的提高；均方根誤差都存在一定幅度的降低。其中，降尺度前后RMSE在5－9月一直保持在較大值，這主要是因?yàn)榱饔蛉甑膹?qiáng)降雨都集中在該時(shí)期，而雷達(dá)捕捉強(qiáng)降水信號(hào)的能力又較差，致使數(shù)據(jù)探測(cè)出現(xiàn)偏差[39]；降尺度后的相對(duì)偏差在各月份表現(xiàn)不一，大多數(shù)月份的結(jié)果較降尺度前有所下降，但1月和7－9月出現(xiàn)增大，原因在于冬季降水形式多為降雪，而夏季常發(fā)生強(qiáng)降雨，TRMM衛(wèi)星對(duì)于這兩種形式的降水捕捉較差[26]，數(shù)據(jù)存在偏差。而降尺度處理可能會(huì)進(jìn)一步低估實(shí)測(cè)降水[40]，從而出現(xiàn)高估實(shí)際降水情況向低估轉(zhuǎn)移?？偟膩碚f，在月尺度上TRMM數(shù)據(jù)降尺度會(huì)改善其數(shù)據(jù)精度，更好地反映研究區(qū)域的真實(shí)降水信息。

表5 TRMM產(chǎn)品降尺度前后定量評(píng)價(jià)指標(biāo)對(duì)比

結(jié)合不同尺度的分析結(jié)果，認(rèn)為降尺度后的數(shù)據(jù)相較于原始TRMM產(chǎn)品數(shù)據(jù)，在提高降水面域化數(shù)據(jù)空間分辨率的同時(shí)，還對(duì)精度有一定的改善，使其能夠更加全面地反映渭河流域降水分布的空間格局。

4 結(jié) 論

本文以渭河流域?yàn)檠芯繀^(qū)，對(duì)熱帶降雨測(cè)量衛(wèi)星（Tropical Rainfall Measuring Mission，TRMM）降水產(chǎn)品在流域內(nèi)的適用性進(jìn)行了分析，引入數(shù)字高程模型（Digital Elevation Model，DEM）、歸一化差分植被指數(shù)（Normalized Difference Vegetation Index，NDVI）、經(jīng)緯度、坡度和坡向等“額外信息”變量因子，從因子作用尺度空間差異性的角度，構(gòu)建了基于多尺度地理加權(quán)回歸（Multi-scale Geographically Weighted Regression，MGWR）模型的低空間分辨率TRMM降水產(chǎn)品和變量因子間的函數(shù)關(guān)系，并將其應(yīng)用于高分辨率變量因子，實(shí)現(xiàn)了TRMM數(shù)據(jù)空間尺度下延，同時(shí)以站點(diǎn)實(shí)測(cè)數(shù)據(jù)為基準(zhǔn)，對(duì)降尺度后的TRMM降水?dāng)?shù)據(jù)進(jìn)行了精度驗(yàn)證。所得結(jié)論如下：

1）以研究區(qū)域站點(diǎn)實(shí)測(cè)降水?dāng)?shù)據(jù)為參考，TRMM衛(wèi)星產(chǎn)品降水?dāng)?shù)據(jù)在年、季和月時(shí)間尺度上均表現(xiàn)出較好的適用性。年尺度TRMM降水與站點(diǎn)實(shí)測(cè)降水之間2=0.807，BIAS=2.909%，RMSE=83.477 mm，表現(xiàn)較好；季尺度上秋季（2=0.847）擬合優(yōu)度最高，夏季均方根誤差最大（RMSE=62.393 mm），四季的相對(duì)偏差均較低；月尺度2介于0.456～0.815之間，6月最小，9月最大，BIAS的絕對(duì)值介于0～14%之間，多數(shù)月份為正值，表明對(duì)降水存在一定的高估現(xiàn)象，RMSE指標(biāo)值域范圍為3.507～39.342 mm，較年、季尺度上數(shù)值偏低。總體而言，TRMM降水?dāng)?shù)據(jù)雖存在一定誤差，但精度相對(duì)較高，可以用于流域降尺度模擬研究。

2）MGWR模型通過對(duì)不同變量因子采用不同的帶寬，顯示出了不同變量作用尺度的差異性。研究表明，渭河流域干濕年份TRMM降水空間分異格局由各類變量因子的多種空間尺度過程共同決定。其中，濕潤(rùn)年的地理位置、NDVI、DEM和坡向?qū)邓捻憫?yīng)十分敏感，其作用尺度呈現(xiàn)局地性（帶寬均為43），而坡度對(duì)降水的影響較為平緩（帶寬為357），接近于全局尺度；干旱年的各變量因子均顯現(xiàn)出差異不大的作用尺度（坡度因子帶寬為52，其余因子帶寬為43），整體上表現(xiàn)為局地影響。從MGWR回歸結(jié)果來看，不論干或濕年份，常數(shù)項(xiàng)和NDVI的回歸系數(shù)在流域尺度上均顯著，是反映降水空間分異的主要作用因素。

3）在流域尺度和站點(diǎn)尺度上對(duì)降尺度數(shù)據(jù)進(jìn)行精度檢驗(yàn)。從流域尺度來看，降尺度后的TRMM數(shù)據(jù)相較實(shí)測(cè)降水?dāng)?shù)據(jù)計(jì)算得到的擬合優(yōu)度和均方根誤差得到有限提高。從站點(diǎn)尺度來看，產(chǎn)品數(shù)據(jù)精度較差的站點(diǎn)降尺度結(jié)果也較為不理想，同時(shí)在時(shí)間尺度上，由于誤差的累積，年尺度數(shù)據(jù)的統(tǒng)計(jì)指標(biāo)值差于月尺度統(tǒng)計(jì)指標(biāo)值。

4）在時(shí)間尺度上對(duì)TRMM降尺度數(shù)據(jù)進(jìn)行分析，發(fā)現(xiàn)濕潤(rùn)年的降尺度結(jié)果略優(yōu)于干旱年份，2和BIAS值均優(yōu)于干旱年；同時(shí)由降尺度前后數(shù)據(jù)空間分布可以看出，降尺度后細(xì)節(jié)特征表現(xiàn)更好；月尺度上降尺度后數(shù)據(jù)統(tǒng)計(jì)指標(biāo)值較降尺度前有所提升?？傮w而言，降尺度后數(shù)據(jù)精度得到一定改善，也能更好地反映研究區(qū)域真實(shí)的降水信息與分布格局。

總體而言，針對(duì)TRMM衛(wèi)星降水產(chǎn)品，本文基于MGWR模型提出了一種行之有效的降水?dāng)?shù)據(jù)空間降尺度方法，能夠提高降水?dāng)?shù)據(jù)空間分辨率，展現(xiàn)更多的細(xì)節(jié)特征，并且保障了數(shù)據(jù)精度和空間分布的一致性，可為渭河流域的水文模擬、干旱監(jiān)測(cè)以及水資源管理等提供相對(duì)可靠的高分辨率降水?dāng)?shù)據(jù)集。

[1] Dobson A P, Bradshaw A, Baker A á. Hopes for the future: Restoration ecology and conservation biology[J]. Science, 1997, 277(5325): 515-522.

[2] Goovaerts P. Geostatistical approaches for incorporating elevation into the spatial interpolation of rainfall[J]. Journal of Hydrology, 2000, 228(1/2): 113-129.

[3] 寧珊，張正勇，劉琳，等. TRMM偏最小二乘降尺度降水模型在新疆不同地貌的適應(yīng)性[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(12)：99-109.

Ning Shan, Zhang Zhengyong, Liu Lin, et al. Adaptability of precipitation estimation method based on TRMM data combined with partial least squares downscaling in different landforms of Xinjiang, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(12): 99-109. (in Chinese with English abstract)

[4] 朱會(huì)義，劉述林，賈紹鳳. 自然地理要素空間插值的幾個(gè)問題[J]. 地理研究，2004，23(4)：425-432.

Zhu Huiyi, Liu Shulin, Jia Shaofeng. Problems of the spatial interpolation of physical geographical elements[J]. Geographical Research, 2004, 23(4): 425-432. (in Chinese with English abstract)

[5] 谷黃河，余鐘波，楊傳國(guó)，等. 衛(wèi)星雷達(dá)測(cè)雨在長(zhǎng)江流域的精度分析[J]. 水電能源科學(xué)，2010，28(8)：3-6.

Gu Huanghe, Yu Zhongbo, Yang Chuanguo, et al. Application of satellite radar observed precipitation to accuracy analysis in Yangtze River Basin[J]. Water Resources and Power, 2010, 28(8): 3-6. (in Chinese with English abstract)

[6] 李哲. 多源降雨觀測(cè)與融合及其在長(zhǎng)江流域的水文應(yīng)用[D]. 北京：清華大學(xué)，2015.

Li Zhe. Multi-source Precipitation Observations and Fusion for Hydrological Applications in the Yangtze River Basin[D]. Beijing: Tsinghua University, 2015. (in Chinese with English abstract)

[7] 譚麗麗，黃峰，喬學(xué)瑾，等. TRMM在海河流域南系的降水估算精度評(píng)價(jià)及其對(duì)SWAT模型的適用性[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(6)：132-141.

Tan Lili, Huang Feng, Qiao Xuejin, et al. Evaluation of TRMM satellite-based rainfall data in southern Haihe River Basin and suitability for SWAT model[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(6): 132-141. (in Chinese with English abstract)

[8] Gao Y, Liu M. Evaluation of high-resolution satellite precipitation products using rain gauge observations over the Tibetan Plateau[J]. Hydrology and Earth System Sciences, 2013, 17(2): 837-849.

[9] 劉昌明，劉文彬，傅國(guó)斌，等. 氣候影響評(píng)價(jià)中統(tǒng)計(jì)降尺度若干問題的探討[J]. 水科學(xué)進(jìn)展，2012，23(3)：427-437.

Liu Changming, Liu Wenbin, Fu Guobin, et al. A discussion of some aspects of statistical downscaling in climate impacts assessment[J]. Advances in Water Science, 2012, 23(3): 427-437. (in Chinese with English abstract)

[10] 王曉杰. 基于TRMM的天山山區(qū)降水降尺度方法及其空間變異特征研究[D]. 石河子：石河子大學(xué)，2013.

Wang Xiaojie. Research on Precipitation Downscaling Method and Its Spatial Variability in Tianshan Mountain Area Based on TRMM[D]. Shihezi: Shihezi University, 2013. (in Chinese with English abstract)

[11] Immerzeel W, Rutten M, Droogers P. Spatial downscaling of TRMM precipitation using vegetative response on the Iberian Peninsula[J]. Remote Sensing of Environment, 2009, 113(2): 362-370.

[12] Jia S, Zhu W, L? A, et al. A statistical spatial downscaling algorithm of TRMM precipitation based on NDVI and DEM in the Qaidam Basin of China[J]. Remote sensing of Environment, 2011, 115(12): 3069-3079.

[13] Shi Y, Song L, Xia Z, et al. Mapping annual precipitation across mainland China in the period 2001–2010 from TRMM3B43 product using spatial downscaling approach[J]. Remote Sensing, 2015, 7(5): 5849-5878.

[14] 杜方洲，石玉立，盛夏. 基于深度學(xué)習(xí)的TRMM降水產(chǎn)品降尺度研究[J]. 國(guó)土資源遙感，2020，32(4)：145-153.

Du Fangzhou, Shi Yuli, Sheng Xia. Research on downscaling of TRMM precipitation products based on deep learning: Exemplified by northeast China[J]. Remote Sensing for Land and Resources, 2020, 32(4): 145-153. (in Chinese with English abstract)

[15] 張寒博，韋夢(mèng)思，覃金蘭，等. 基于不同植被指數(shù)的TRMM數(shù)據(jù)降尺度及誤差校正研究[J]. 水土保持通報(bào)，2021，41(4)：214-223.

Zhang Hanbo, Wei Mengsi, Qin Jinlan, et al. Downscaling and error correction of TRMM data based on different vegetation indices[J]. Bulletin of Soil and Water Conservation, 2021, 41(4): 214-223. (in Chinese with English abstract)

[16] Chen C, Zhao S, Duan Z, et al. An improved spatial downscaling procedure for TRMM 3B43 precipitation product using geographically weighted regression[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(9): 4592-4604.

[17] Zeng Z, Qian S, Plaza J, et al. Spatial downscaling for global precipitation measurement using a geographically and temporally weighted regression model[C]. Proceedings of The IGARSS 2020-2020 IEEE International Geoscience And Remote Sensing Symposium, Waikoloa, HI, USA, IEEE，2020: 5348-5351.

[18] 沈體雁，于瀚辰，周麟，等. 北京市二手住宅價(jià)格影響機(jī)制：基于多尺度地理加權(quán)回歸模型（MGWR）的研究[J]. 經(jīng)濟(jì)地理，2020，40(3)：75-83.

Shen Tiyan, Yu Hanchen, Zhou Lin, et al. On hedonic price of second-hand houses in Beijing based on muti-scale geographically weighted regression: Scale law of spatial heterogeneity[J]. Economic Geography, 2020, 40(3): 75-83. (in Chinese with English abstract)

[19] Fotheringham A S, Yang W, Kang W. Multiscale geographically weighted regression (MGWR)[J]. Annals of the American Association of Geographers, 2017, 107(6): 1247-1265.

[20] 趙海根，楊勝天，周旭，等. TRMM 3B42降雨數(shù)據(jù)在渭河流域的應(yīng)用分析[J]. 水文，2015，35(1)：61-67.

Zhao Haigen, Yang Shengtian, Zhou Xu, et al. Application analysis of TRMM 3B42 precipitation data in Weihe River Basin[J]. Journal of China Hydrology, 2015, 35(1): 61-67. (in Chinese with English abstract)

[21] 郭雯雯，黃生志，趙靜，等. 渭河流域潛在蒸散發(fā)時(shí)空演變與驅(qū)動(dòng)力量化分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(3)：81-89.

Guo Wenwen, Huang Shengzhi, Zhao Jing, et al. Spatio-temporal dynamics and driving forces of potential evapotranspiration in the Wei River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(3): 81-89. (in Chinese with English abstract)

[22] 張洪波，顧磊，陳克宇，等. 渭河生態(tài)水文聯(lián)系變異分區(qū)研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào): 自然科學(xué)版，2016，44(6)：210-220.

Zhang Hongbo, Gu Lei, Chen Keyu, et al. Division of eco-hydrologic connection variability in Wei River basin[J]. Journal of Northwest A&F University (Natural Science Edition), 2016, 44(6): 210-220. (in Chinese with English abstract)

[23] 張?jiān)聢A. 紅河流域TRMM降水?dāng)?shù)據(jù)空間降尺度研究[D]. 昆明：云南大學(xué)，2018.

Zhang Yueyuan. Spatial Downscaling of TRMM Precipitation Data in the Red River Basin, China[D]. Kunming: Yunnan University, 2018. (in Chinese with English abstract)

[24] Yu H, Fotheringham A S, Li Z, et al. Inference in multiscale geographically weighted regression[J]. Geographical Analysis, 2020, 52(1): 87-106.

[25] 王祎婷，謝東輝，李亞惠. 光譜指數(shù)趨勢(shì)面的城市地表溫度降尺度轉(zhuǎn)換[J]. 遙感學(xué)報(bào)，2014，18(6)：1169-1181.

Wang Yiting, Xie Donghui, Li Yahui. Downscaling remotely sensed land surface temperature over urban areas using trend surface of spectral index[J]. Journal of Remote Sensing, 2014, 18(6): 1169-1181. (in Chinese with English abstract)

[26] 范田億，張翔，黃兵，等. TRMM衛(wèi)星降水產(chǎn)品降尺度及其在湘江流域水文模擬中的應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(15)：179-188.

Fan Tianyi, Zhang Xiang, Huang Bing, et al. Downscaling of TRMM satellite precipitation products and its application in hydrological simulation of Xiangjiang River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(15): 179-188. (in Chinese with English abstract)

[27] 杜軍凱，賈仰文，李曉星，等. 基于TRMM衛(wèi)星降水的太行山區(qū)降水時(shí)空分布格局[J]. 水科學(xué)進(jìn)展，2019，30(1)：1-13.

Du Junkai, Jia Yangwen, Li Xiaoxing, et al. Study on the spatial-temporal distribution pattern of precipitation in the Taihang Mountain region using TRMM data[J]. Advances in Water Science, 2019, 30(1): 1-13. (in Chinese with English abstract)

[28] 文遷，譚國(guó)良，羅嗣林. 降水分布受地形影響的分析[J]. 水文，1997(增刊1)：64-66.

[29] 鄭杰，閭利，馮文蘭，等. 基于TRMM 3B43數(shù)據(jù)的川西高原月降水量空間降尺度模擬[J].中國(guó)農(nóng)業(yè)氣象，2016，37(2)：245-254.

Zheng Jie, Yan Li, Feng Wenlan, et al. Spatial downscaling simulation of monthly precipitation in western Sichuan Plateau based on TRMM 3B43 data[J]. Chinese Journal of Agricultural Meteorology, 2016, 37(2): 245-254. (in Chinese with English abstract)

[30] 姜琳，馮文蘭，郭兵. 雅魯藏布江流域近13年植被覆蓋動(dòng)態(tài)監(jiān)測(cè)及與降水因子的相關(guān)性分析[J]. 長(zhǎng)江流域資源與環(huán)境，2014，23(11)：1610-1619.

Jiang Lin, Feng Wenlan, Guo Bing. Analysis of dynamic monitoring of vegetation change and the correlation with precipitation factor in Yalu Tsangpo River basin during the past 13 years[J]. Resources and Environment in the Yangtze Basin, 2014, 23(11): 1610-1619. (in Chinese with English abstract)

[31] 白建軍，白江濤，王磊. 2000～2010年陜北地區(qū)植被NDVI時(shí)空變化及其與區(qū)域氣候的關(guān)系[J]. 地理科學(xué)，2014，34(7)：882-888.

Bai Jianjun, Bai Jiangtao, Wang Lei. Spatio-temporal change of vegetation NDVI and its relations with regional climate in northern Shaanxi province in 2000-2010[J]. Scientia Geographica Sinica, 2014, 34(7): 882-888. (in Chinese with English abstract)

[32] 白淑英，王莉，史建橋. 長(zhǎng)江流域NDVI對(duì)氣候變化響應(yīng)的時(shí)滯效應(yīng)[J]. 中國(guó)農(nóng)業(yè)氣象，2012，33(4)：579-586.

Bai Shuying, Wang Li, Shi Jianqiao. Time lag effect of NDVI response to climatic change in Yangtze River Basin[J]. Chinese Journal of Agrometeorology, 2012, 33(4): 579-586. (in Chinese with English abstract)

[33] 任亮. 陜西秦巴山區(qū)多源遙感降水?dāng)?shù)據(jù)降尺度研究[D]. 西安：陜西師范大學(xué)，2018.

[34] Sokol Z, Bli?ňák V. Areal distribution and precipitation-altitude relationship of heavy short-term precipitation in the Czech Republic in the warm part of the year[J]. Atmospheric Research, 2009, 94(4): 652-62.

[35] 高維英，李明，徐新田，等. 2003年渭河流域連陰雨過程中的暴雨成因分析[M]. 北京：氣象出版社，2004：188-190.

[36] 張杰，李棟梁，何金梅，等. 地形對(duì)青藏高原豐枯水年雨季降水量空間分布的影響[J]. 水科學(xué)進(jìn)展，2007，18(3)：319-326.

Zhang Jie, Li Dongliang, He Jinmei, et al. Influence of terrain on precipitation distribution in Qingzang tableland in wet and dry years[J]. Advances in Water Science, 2007, 18(3): 319-326. (in Chinese with English abstract)

[37] 顧文亞. BEMD地形分解支持的不同尺度地形下中國(guó)月降水分布網(wǎng)格化模擬研究[D]. 南京：南京信息工程大學(xué)，2021.

Gu Wenya. Study on Grid-based Simulation of Monthly Precipitation Distribution in China over the Different Scale Terrain Based on BEMD Terrain Decomposition[D]. Nanjing: Nanjing University of Information Science & Technology, 2021. (in Chinese with English abstract)

[38] 馬金輝，屈創(chuàng)，張海筱，等. 2001-2010年石羊河流域上游TRMM降水資料的降尺度研究[J]. 地理科學(xué)進(jìn)展，2013，32(9)：1423-1432.

Ma Jinhui, Qu Chuang, Zhang Haixiao, et al. Spatial downscaling of TRMM precipitation data based on DEM in the upstream of Shiyang River basin during 2001-2010[J]. Progress in Geography, 2013, 32(9): 1423-1432. (in Chinese with English abstract)

[39] 劉小嬋. TRMM降水?dāng)?shù)據(jù)的空間降尺度研究[D]. 長(zhǎng)春：東北師范大學(xué)，2015.

Liu Xiaochan. Spatial Downscaling of TRMM Precipitation Data[D]. Changchun: Northeast Normal University, 2015. (in Chinese with English abstract)

[40] 熊俊楠，李偉，劉志奇，等. 基于GWR模型的青藏高原地區(qū)TRMM數(shù)據(jù)降尺度研究[J]. 國(guó)土資源遙感，2019，31(4)：88-95.

Xiong Junnan, Li Wei, Liu Zhiqi, et al. Research on downscaling of TRMM data in the Tibetan Plateau based on GWR model[J]. Remote Sensing for Land and Resources, 2019, 31(4): 88-95. (in Chinese with English abstract)

A MGWR-based spatial downscaling for TRMM precipitation in the Weihe River Basin

Li Yangbing1, Zhang Hongbo1,2※, Yang Tianzeng1, Lyu Fengguang1, Wang Yuwei1, Yao Congcong1

(1710054,; 2.710054)

High-resolution precipitation data can be directly used to characterize the spatial-temporal differentiation features of precipitation after the accurate simulation of hydrological, meteorological, and biological systems. Therefore, it is crucial to implement the spatial downscaling for the precipitation products with the low spatial resolution. However, it is still lacking on the precision and detail features in the downscaling precipitation data, due to fail to consider the scale variations in the spatial distribution of precipitation. In this study, a spatial downscaling approach was proposed to improve the TRMM precipitation data in the Wei River basin (WRB). A Multi-scale Geographically Weighted Regression (MGWR) model was also integrated to enable the conditional relationships between the response and predictor variables that changed at the spatial scales. Therein, the goodness of fit (2), relative deviation (BIAS), and Root Mean Square Error (RMSE) were employed to verify the TRMM satellite precipitation product, compared with the actual precipitation data from meteorological stations. Normalized Difference Vegetation Index (NDVI), digital elevation model (DEM), slope, aspect, latitude and longitude were induced as the Geographic Environmental Factors (GEFs). The MGWR models with the monthly TRMM precipitation data were constructed to further investigate the scale effects of factors on precipitation distribution. The spatial downscaling of TRMM production data was then implemented using a scale conversion process. Finally, the reliable spatial downscaling was achieved in the TRMM precipitation products. The results illustrated that: 1) The TRMM precipitation data was better suited for use at different scales in the WRB. An acceptable fitness was found in the statistics of2(0.807), BIAS (2.909%), and RMSE (83.477 mm) at the annual scale. Specifically, the maximum2was 0.847 at the seasonal scale, the largest RMSE was 62.393 mm in the summer, and the BIAS values were lower in all four seasons. More importantly, the2varied between 0.456 and 0.815 on the monthly scale, with the smallest value in June and the largest value in September. The BIAS was positive in the most months, indicating that the TRMM product data generally overestimated the precipitation. The RMSE index was fallen in the range of 3.507-39.342 mm, which was lower than those on the annual and seasonal scales. 2) Different scale characteristics were found in the influence of various GEFs on the spatial pattern of precipitation divergence in wet and dry years. Slope was set as the global scale, whereas the DEM, NDVI, aspect, latitude, and longitude were the local effects on the precipitation in wet years, and all GEFs were used the local impacts in dry years. 3) The more precise data was obtained in the downscaled TRMM on the watershed and station scales, compared with the product data, indicating an increase in the2of the entire watershed of 3%, while a decrease in the RMSE of 1 mm. However, the accuracy of station downscaling precipitation data at the annual scale was worse than that at the monthly scale, due to the accumulation of errors on the time scale, as shown by the2range decreasing, while the RMSE range increasing. 4) The downscaled TRMM data presented the better detailed characteristics, the greater precision at annual and monthly scales, and a much more delicate geographical distribution than the product data. The finding can provide a strong data support for the hydrological design in the areas with less precipitation data.

precipitation;TRMM; downscaling; MGWR model; scale difference

10.11975/j.issn.1002-6819.2022.23.015

P426.61; P237

A

1002-6819(2022)-23-0141-11

黎揚(yáng)兵，張洪波，楊天增，等. 基于MGWR的渭河流域TRMM降水產(chǎn)品空間降尺度分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(23)：141-151.doi：10.11975/j.issn.1002-6819.2022.23.015 http://www.tcsae.org

Li Yangbing, Zhang Hongbo, Yang Tianzeng, et al. A MGWR-based spatial downscaling for TRMM precipitation in the Weihe River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(23): 141-151. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2022.23.015 http://www.tcsae.org

2022-06-14

2022-11-13

國(guó)家自然科學(xué)基金項(xiàng)目（51979005）；陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃項(xiàng)目（2020JM-250）

黎揚(yáng)兵，研究方向?yàn)樗乃Y源。Email：2020129016@chd.edu.cn

張洪波，博士，教授，研究方向?yàn)樗乃Y源。Email：hbzhang@chd.edu.cn