趙珊珊 周向陽 童錦施

摘要：為評(píng)估中國(guó)大陸區(qū)域典型長(zhǎng)序列、高時(shí)間分辨率降水產(chǎn)品（MSWEP，TRMM，CMFD）在不同尺度的精度水平、演變特征和空間差異，基于2 117個(gè)氣象站近5 a的逐小時(shí)數(shù)據(jù)，選取相關(guān)系數(shù)（Corr）、均方根誤差（RMSE）、相對(duì)誤差中位數(shù)（REM）和擊中率（POD）4個(gè)指標(biāo)，系統(tǒng)地評(píng)估了降水產(chǎn)品在小時(shí)、日、旬、月尺度的精度及空間分布特征，建立了相關(guān)指標(biāo)隨時(shí)間尺度增加而演變的函數(shù)模型。結(jié)果表明：3組降水產(chǎn)品精度的空間分布總體表現(xiàn)為由西北逐漸向東南區(qū)域增加的特征，且POD和REM在各尺度均為MSWEP明顯優(yōu)于 CMFD和TRMM，但Corr、RMSE在亞日和旬月尺度精度表現(xiàn)不一致；多尺度變化方面，RMSE和REM呈冪函數(shù)變化，而Corr和POD服從修正后的對(duì)數(shù)Gamma分布函數(shù)。研究成果可為降水產(chǎn)品的優(yōu)選識(shí)別、精度提升和深度融合提供依據(jù)。

關(guān)鍵詞：降水產(chǎn)品； 精度評(píng)價(jià)； 尺度演變； 空間分布； 函數(shù)模型

中圖法分類號(hào)： P412.27

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.04.015

0引 言

降水?dāng)?shù)據(jù)是氣象、水文和地質(zhì)災(zāi)害等領(lǐng)域監(jiān)測(cè)預(yù)警的重要基礎(chǔ)資料［1-3］，其精度和時(shí)空分辨率直接影響模型的預(yù)測(cè)結(jié)果［4］。降水?dāng)?shù)據(jù)可通過地面采集和遙感監(jiān)測(cè)獲得。前者精度高但受地形、氣象站點(diǎn)數(shù)目（布局）、降水空間異質(zhì)性的限制，在一些區(qū)域因站點(diǎn)分布稀疏且觀測(cè)建設(shè)較晚，從而缺乏長(zhǎng)序列、高分辨率和高精度的地面觀測(cè)降水?dāng)?shù)據(jù)［5］，并且高時(shí)間分辨率以及高空間覆蓋率數(shù)據(jù)也難以獲?。?］。隨著遙感技術(shù)的飛速發(fā)展，網(wǎng)格化降水產(chǎn)品的出現(xiàn)為降水量估算［7］、徑流模擬［8］、干旱檢測(cè)［9］等方面提供了重要支撐。因此，評(píng)價(jià)長(zhǎng)序列、高時(shí)間分辨率降水產(chǎn)品的精度并優(yōu)選數(shù)據(jù)，將促進(jìn)上述工作更好地開展。

針對(duì)降水產(chǎn)品的精度問題，很多研究基于地面實(shí)測(cè)降水?dāng)?shù)據(jù)對(duì)不同類型和用途的降水?dāng)?shù)據(jù)集進(jìn)行評(píng)估。在全球尺度，通過分析8種常用降水產(chǎn)品的日尺度數(shù)據(jù)在中國(guó)、歐洲和北美的1 382個(gè)流域的效果可知，MSWEP V2.0的面雨量比地面站網(wǎng)的多點(diǎn)數(shù)據(jù)擁有更好的水文性能［10］。在中國(guó)不同降水產(chǎn)品精度在特定區(qū)域的表現(xiàn)如下：在長(zhǎng)江流域總體表現(xiàn)為CMFD＞ERA5＞MSWEP［11］，TRMM 3B43產(chǎn)品能更好地描述鄱陽湖流域的年度降水空間分布［12］；在渭河流域，對(duì)6種降水產(chǎn)品描述極端降水的能力評(píng)價(jià)結(jié)果表明，TRMM 3B42 V7在監(jiān)測(cè)極端降水時(shí)具有更高的精度和更小的空間不確定性［13］；在柴達(dá)木盆地，對(duì)比MSWEP V6、GPM IMERG V3和TRMM 3B42降水的準(zhǔn)確性可知，MSWEP在月和年尺度上精度最好［14］。在中國(guó)中東部率水河流域?qū)PM和TRMM進(jìn)行了多時(shí)間尺度精度評(píng)估，結(jié)果表明二者無法有效觀測(cè)小時(shí)降水量［15］。在中國(guó)大陸的不同氣候帶，IMEGE、TRMM 3B42數(shù)據(jù)在日、月、年的精度可通過線性關(guān)系進(jìn)行修正［16］。

上述研究多針對(duì)某區(qū)域進(jìn)行評(píng)價(jià)，時(shí)間尺度多在日尺度及以上，少數(shù)研究評(píng)估了小時(shí)尺度上降水?dāng)?shù)據(jù)集的精度；一些評(píng)估TRMM數(shù)據(jù)精度的文獻(xiàn)并未將其與時(shí)區(qū)進(jìn)行統(tǒng)一，因?yàn)門RMM 3B42 V7的3 h數(shù)據(jù)時(shí)間為世界時(shí)加減90 min［17］。中國(guó)幅員遼闊、地形地貌復(fù)雜、氣候多樣，一些長(zhǎng)序列的降水產(chǎn)品精度如何、典型誤差指標(biāo)隨著尺度的變化如何轉(zhuǎn)換、在不同氣候分帶條件下哪一種產(chǎn)品最優(yōu)等問題仍需要進(jìn)一步評(píng)價(jià)。因此，本文選取時(shí)間分辨率高、跨度長(zhǎng)的3種典型降水產(chǎn)品（分別為MSWEP V2.8，TRMM 3B42 V7和China Meteorological Forcing Dataset（CMFD）［18］），基于中國(guó)的氣候分帶特征分別評(píng)估其在不同區(qū)域的精度水平，揭示典型誤差指標(biāo)隨小時(shí)至月尺度的演變函數(shù)關(guān)系。研究成果可為降水產(chǎn)品的優(yōu)選、進(jìn)一步深度融合以及氣象、水文、地質(zhì)等災(zāi)害的防治提供參考依據(jù)。

1研究區(qū)及數(shù)據(jù)

1.1數(shù) 據(jù)

（1） 降水產(chǎn)品及處理。

研究選取3種典型長(zhǎng)序列降水產(chǎn)品，分別為MSWEP V2.8、CMFD以及TRMM 3B42 V7，詳細(xì)信息如表1所列。

（2） 地面氣象站降水資料。

研究選取國(guó)家氣象科學(xué)數(shù)據(jù)中心提供的中國(guó)國(guó)家級(jí)地面站小時(shí)值數(shù)據(jù)，包括中國(guó)31個(gè)省級(jí)行政區(qū)（不包含港澳臺(tái)地區(qū)）［19］，這是目前國(guó)內(nèi)最佳的小時(shí)降水?dāng)?shù)據(jù)［20］。降水?dāng)?shù)據(jù)由各氣象站點(diǎn)的虹吸管或翻斗式雨量計(jì)自動(dòng)收集每小時(shí)降水?dāng)?shù)據(jù)所得，并且數(shù)據(jù)都經(jīng)過了CMA 3個(gè)級(jí)別的嚴(yán)格質(zhì)量控制，包括極值檢查、內(nèi)部一致性檢查、空間一致性檢查［21］。數(shù)據(jù)的監(jiān)測(cè)時(shí)段為2015年6月16日至2020年12月31日。由于數(shù)據(jù)下載不及時(shí)而存在少量缺失，所選取測(cè)站的逐小時(shí)降水?dāng)?shù)據(jù)實(shí)有率均大于85%，共計(jì)采用2 117個(gè)測(cè)站，其空間分布如圖1所示。

（3） 時(shí)間尺度匹配。

由于3種降水產(chǎn)品3 h數(shù)據(jù)時(shí)間為世界時(shí)，故實(shí)測(cè)數(shù)據(jù)也相應(yīng)調(diào)整為3 h世界時(shí)數(shù)據(jù)。

另外，TRMM 3 h數(shù)據(jù)不是整時(shí)數(shù)據(jù)，例如“3B42.19980101.00.7.HDF”數(shù)據(jù)的記錄時(shí)段為1998年1月1日的22：30到次日01：30，而所采用的數(shù)據(jù)為整小時(shí)記錄。為減小二者的偏差，在不同尺度降水?dāng)?shù)據(jù)匹配時(shí)，將逐小時(shí)地面數(shù)據(jù)視為均勻分布，平均分配到兩個(gè)相關(guān)聯(lián)的尺度。TRMM日數(shù)據(jù)實(shí)為當(dāng)日22：30到次日22：30，月數(shù)據(jù)也偏差1.5 h，如TRMM 7月1日數(shù)據(jù)為6月30日22：30至8月1日22：30。

1.2氣候分帶

根據(jù)地理位置、水分情況和溫度帶，將中國(guó)分為7個(gè)氣候分區(qū)［22］，分別為：Ⅰ東北濕潤(rùn)半濕潤(rùn)溫帶地區(qū)、Ⅱ華北濕潤(rùn)半濕潤(rùn)暖溫帶地區(qū)、Ⅲ華中華南濕潤(rùn)亞熱帶地區(qū)、Ⅳ華南熱帶濕潤(rùn)地區(qū)、Ⅴ內(nèi)蒙溫帶草原地區(qū)、Ⅵ西北溫帶及暖溫帶荒漠地區(qū)和Ⅶ青藏高原地區(qū)，如圖1所示。其中Ⅰ、Ⅱ、Ⅲ、Ⅳ屬于以濕潤(rùn)的季風(fēng)氣候?yàn)橹鞯臇|部季風(fēng)區(qū)，Ⅴ、Ⅵ屬于以半干旱和干旱氣候?yàn)橹鞯奈鞅备珊祬^(qū)，Ⅶ屬于具有氣候垂直變化特點(diǎn)的青藏高寒區(qū)。

1.3評(píng)價(jià)方法

1.3.1數(shù)據(jù)采樣與尺度選擇

（1） 數(shù)據(jù)采樣。

地面站點(diǎn)所對(duì)應(yīng)的降水產(chǎn)品數(shù)據(jù)，采用雨量站所在單元格的數(shù)據(jù)。

（2） 尺度選擇。

基于降水產(chǎn)品的最高時(shí)間分辨率，研究采用從3 h至1月的典型尺度，分別為：3，6，12，24 h和10 d，1月。

1.3.2精度評(píng)價(jià)因子

研究采用相關(guān)系數(shù)Corr、均方根誤差RMSE、相對(duì)誤差中位數(shù)REM和擊中率POD（probability of detection）進(jìn)行精度評(píng)價(jià)，如公式（1）～（4）所示。

點(diǎn)所對(duì)應(yīng)時(shí)間的降水值，S、G分別是衛(wèi)星、地面站點(diǎn)所對(duì)應(yīng)時(shí)間序列的平均降水值；N為降水序列長(zhǎng)度；RE為相對(duì)誤差，REM為相對(duì)誤差中位數(shù)；N11表示氣象站點(diǎn)和降水產(chǎn)品均能監(jiān)測(cè)到的降水事件個(gè)數(shù)，N10表示氣象站點(diǎn)能夠而降水產(chǎn)品不能監(jiān)測(cè)到的降水事件個(gè)數(shù)。

1.3.3區(qū)域平均

研究選取加權(quán)平均的方法計(jì)算全國(guó)及不同氣候分區(qū)的評(píng)價(jià)誤差水平，權(quán)重采用泰森多邊形法則進(jìn)行確定［23］。

1.3.4誤差的描述

基于區(qū)域各站點(diǎn)的4個(gè)誤差指標(biāo)，分別分析其特征統(tǒng)計(jì)量-中位數(shù)指標(biāo)。

1.3.5演變模型

2結(jié)果分析

2.1相關(guān)系數(shù)

2.1.1空間分布特征

3種降水產(chǎn)品在不同驗(yàn)證站點(diǎn)的相關(guān)系數(shù)空間分布特征如圖2所示。結(jié)果表明：相關(guān)系數(shù)總體呈現(xiàn)出由西北至東南逐漸增加的特征，但在不同時(shí)間尺度也有所差異。根據(jù)圖2中的氣候分帶，統(tǒng)計(jì)各站點(diǎn)相關(guān)系數(shù)的中位數(shù)，如表2所列，其特征差異主要表現(xiàn)為3個(gè)方面。

（1） 3組降水產(chǎn)品在不同氣候帶的差異顯著，相關(guān)系數(shù)總體呈現(xiàn)出分區(qū)Ⅰ-Ⅳ高于分區(qū)Ⅴ-Ⅵ的特征，即濕潤(rùn)半濕潤(rùn)地區(qū)明顯高于干旱半干旱地區(qū)。

（2） 隨時(shí)間尺度增加，精度較差的Ⅴ內(nèi)蒙溫帶草原地區(qū)和Ⅶ青藏高原地區(qū)的相關(guān)系數(shù)差異與精度較高的4個(gè)分區(qū)差異逐漸縮小，而Ⅵ西北荒漠地區(qū)的差異仍較大。

（3） 3組降水產(chǎn)品在全國(guó)區(qū)域不同時(shí)間尺度的相關(guān)系數(shù)表現(xiàn)有所不同，24 h以下尺度MSWEP數(shù)據(jù)的中位數(shù)不超過0.62，CMFD不超過0.58，TRMM則小于等于0.55，這表明基于相關(guān)系數(shù)的評(píng)估在24 h及更小尺度可優(yōu)先考慮MSWEP。在旬和月尺度，3組產(chǎn)品的相關(guān)系數(shù)明顯提升，且表現(xiàn)為CMFD高于其余兩組，原因是各產(chǎn)品融合的數(shù)據(jù)源不同。

2.1.2尺度變化關(guān)系

3組降水產(chǎn)品相關(guān)系數(shù)在全國(guó)及不同氣候分區(qū)的尺度演變特征如圖3所示。統(tǒng)計(jì)模型的4個(gè)參數(shù)和R2，如表3所列。結(jié)果表明：通過矯正后的對(duì)數(shù)Gamma分布能夠較好地描述相關(guān)系數(shù)的演變規(guī)律，模型的R2除Ⅵ西北地區(qū)的MSWEP呈較大偏差外，其余地區(qū)均大于0.94，且2/3的測(cè)試結(jié)果超過0.99。進(jìn)一步分析模型參數(shù)的特征表明，MSWEP數(shù)據(jù)總體表現(xiàn)為更大的變程、更小的初始值、更小的形狀因子和更大的尺度因子；CMFD和TRMM則表現(xiàn)為更小的變程、尺度因子和更大的初始值、形狀因子。更大的變程R意味著更小的初始值，反映出在小尺度的相關(guān)系數(shù)較小，但后期增長(zhǎng)率較大；相關(guān)形狀因子都大于1，更大的形狀因子表明MSWEP產(chǎn)品相關(guān)系數(shù)演變特征一階導(dǎo)數(shù)的正態(tài)性更好；更大的尺度因子也反映出相同形狀因子時(shí)其拖尾性更好，在更大尺度趨近于1的速率更緩。3組產(chǎn)品的各參數(shù)在Ⅲ華中華南濕潤(rùn)亞熱帶地區(qū)和Ⅳ華南熱帶濕潤(rùn)地區(qū)表現(xiàn)出總體比較接近的特征。綜合分析表明，MSWEP產(chǎn)品相關(guān)系數(shù)隨尺度的演變更加復(fù)雜。

2.2均方根誤差

2.2.1空間分布特征

3組降水產(chǎn)品在全國(guó)及不同區(qū)域的均方根誤差的空間分布總體呈現(xiàn)出由西北至東南逐漸增加的特征，如圖4所示。但在不同氣候分區(qū)，均方根誤差和相關(guān)系數(shù)的表現(xiàn)特征不同。相關(guān)系數(shù)高的濕潤(rùn)半濕潤(rùn)地區(qū)，其均方根誤差明顯高于干旱半干旱地區(qū)。均方根誤差最大的區(qū)域?yàn)棰羧A南熱帶濕潤(rùn)地區(qū)，最小的為Ⅵ西北溫帶及暖溫帶荒漠地區(qū)。全國(guó)及各氣候分區(qū)站點(diǎn)的均方根誤差典型分位數(shù)統(tǒng)計(jì)結(jié)果如表2所列。結(jié)果表明：在全國(guó)的總體表現(xiàn)特征和相關(guān)系數(shù)一致，在24 h及更小時(shí)間尺度呈現(xiàn)出MSWEP＜CMFD＜TRMM，其中位數(shù)分別由3 h的5.63，6.14 mm和6.17 mm增加至24 h的11.21，11.61 mm和12.73 mm，而旬和月尺度基本呈現(xiàn)出CMFD＜TRMM＜MSWEP，3組產(chǎn)品的中位數(shù)分別為旬尺度的19.07，23.86，21.19 mm，以及月尺度的28.95，33.12，36.04 mm。

2.2.2尺度變化關(guān)系

選取合適的曲線對(duì)3組降水產(chǎn)品的均方根誤差隨尺度的演變關(guān)系進(jìn)行擬合，結(jié)果如圖5所示。結(jié)果表明：均方根誤差隨時(shí)間尺度增加呈現(xiàn)出良好的冪函數(shù)演變特征。模型的R2可達(dá)到0.94及以上，低值區(qū)主要位于Ⅵ西北干旱區(qū)，而Ⅲ華中華南濕潤(rùn)亞熱帶、Ⅳ華南熱帶濕潤(rùn)地區(qū)和Ⅶ青藏高原地區(qū)R2可達(dá)到0.99。另一方面，3組降水產(chǎn)品的冪指數(shù)和系數(shù)均有所不同，總體表現(xiàn)為MSWEP的冪指數(shù)最高，TRMM次之，CMFD最小。需要注意的是，這種特征并不能說明哪一組數(shù)據(jù)最優(yōu)，因均方根誤差隨尺度增加的演變關(guān)系既受降水特征影響，也需要服從統(tǒng)計(jì)學(xué)規(guī)律，因此需要進(jìn)一步深入討論。

2.3擊中率

2.3.1空間分布特征

3組產(chǎn)品的擊中率總體呈現(xiàn)為濕潤(rùn)半濕潤(rùn)區(qū)域的精度高于半干旱、干旱區(qū)域的精度，且精度最差的均在Ⅵ西北溫帶及暖溫帶荒漠地區(qū)，但表現(xiàn)最優(yōu)的區(qū)域有所不同，結(jié)果如圖6所示。MSWEP在Ⅳ華南熱帶濕潤(rùn)地區(qū)最優(yōu)；CMFD最優(yōu)位于Ⅲ華中華南濕潤(rùn)亞熱帶地區(qū)；TRMM在12 h及更小時(shí)間尺度上的Ⅱ華北濕潤(rùn)半濕潤(rùn)暖溫帶地區(qū)表現(xiàn)最優(yōu)，在24 h及更大時(shí)間尺度上的青藏高原地區(qū)表現(xiàn)較好。3組降水產(chǎn)品擊中率的統(tǒng)計(jì)結(jié)果如表4所列。結(jié)果表明：① 3組產(chǎn)品的擊中率總體呈現(xiàn)出MSWEP＞CMFD＞TRMM，且隨著時(shí)間尺度的減小差異逐漸增加。② MSWEP的擊中率精度較高。在3 h尺度的全國(guó)平均中位數(shù)水平為0.77，6 h尺度增加至0.83，24 h尺度則可達(dá)到0.94。就擊中率而言，MSWEP的精度明顯高于其余兩組，相關(guān)結(jié)果將為數(shù)據(jù)的融合提供參考依據(jù)。

2.3.2變化關(guān)系

對(duì)數(shù)Gamma分布函數(shù)能夠很好地描述3組降水產(chǎn)品擊中率隨時(shí)間尺度增加的演變特征，其擬合結(jié)果和相關(guān)模型參數(shù)分別如圖7和表3所示。模型R2的最小值為0.97，一般的站點(diǎn)達(dá)到或超過0.99。所反演的模型參數(shù)表明：MSWEP產(chǎn)品擊中率的變程最小，初始值最高，形狀因子和尺度因子略低于TRMM但比較接近；CMFD則呈現(xiàn)出最大的變程、最小的初始值、最小的形狀因子和最大的尺度因子；TRMM表現(xiàn)為較大的變程、較小的初始值、最大的形狀因子和較小的尺度因子。這也反映出不同產(chǎn)品的演變特征，MSWEP的總體精度最高，TRMM和MSWEP演變特征一階導(dǎo)數(shù)的對(duì)數(shù)正態(tài)性較好、拖尾性較差，而CMFD則呈現(xiàn)出相反的特征，演變模式更加復(fù)雜。

2.4相對(duì)誤差中位數(shù)

2.4.1空間分布特征

3組降水產(chǎn)品的相對(duì)誤差中位數(shù)空間分布如圖8所示，在不同區(qū)域的相對(duì)誤差中位數(shù)統(tǒng)計(jì)結(jié)果如表2所列。結(jié)果表明：

（1） 在24 h及更小尺度，絕大多數(shù)站點(diǎn)的相對(duì)誤差中位數(shù)都為1或接近1。其原因是此時(shí)地面發(fā)生弱降水（一般為0.01 mm），而相關(guān)降水產(chǎn)品并未感知。結(jié)合降水產(chǎn)品的擊中率分析表明，除MSWEP外，CMFD和TRMM兩組產(chǎn)品的擊中率在該尺度范圍較低（尤其是12 h及更小尺度），產(chǎn)品預(yù)測(cè)值為0而實(shí)測(cè)值為0.1 mm，相對(duì)誤差表現(xiàn)為1且占比較高；并且一些降水產(chǎn)品在地面發(fā)生較弱降水時(shí)的值明顯高于地面值，呈現(xiàn)出相對(duì)誤差大于1，導(dǎo)致大多數(shù)站點(diǎn)的相對(duì)誤差中位數(shù)為1，最終導(dǎo)致不同區(qū)域站點(diǎn)分析出的相對(duì)誤差中位數(shù)多呈現(xiàn)為1。

（2） 在旬和月尺度，降水產(chǎn)品的相對(duì)誤差中位數(shù)顯著降低且呈現(xiàn)南部精度總體好于北部。在旬尺度，表現(xiàn)為MSWEP略優(yōu)于CMFD，TRMM最差，三者的中位數(shù)分別為0.44，0.46和0.76，而月尺度則CMFD略優(yōu)于MSWEP，TRMM最差，三者的中位數(shù)分別為 0.27，0.30和0.37。空間分布方面3組產(chǎn)品均表現(xiàn)為在Ⅵ西北溫帶及暖溫帶荒漠地區(qū)明顯低于另外6個(gè)分區(qū)。

2.4.2尺度變化

3組降水產(chǎn)品相對(duì)誤差中位數(shù)隨時(shí)間尺度增加的演變關(guān)系如圖9所示?？傮w而言，其變化趨勢(shì)可用冪函數(shù)進(jìn)行擬合，但在不同產(chǎn)品的模型表現(xiàn)差異較大。總體表現(xiàn)最優(yōu)的為MSWEP，R2從3 h至月尺度均大于 0.9；CMFD的表現(xiàn)次之，R2在不同尺度均大于0.8；TRMM表現(xiàn)較差，除在華中華南濕潤(rùn)地區(qū)的R2達(dá)到0.9以外，其余分區(qū)介于0.54～0.72之間。TRMM表現(xiàn)較差的原因是該產(chǎn)品在24 h及更小尺度的擊中率精度較差，從而導(dǎo)致各站點(diǎn)的相對(duì)誤差中位數(shù)多為1，對(duì)各站點(diǎn)的經(jīng)驗(yàn)頻率進(jìn)行分析時(shí)發(fā)現(xiàn)其中位數(shù)也多為1。進(jìn)一步對(duì)比MSWEP和CMFD的模型參數(shù)可知，后者的冪指數(shù)均小于前者，說明其相對(duì)誤差的衰減率隨尺度增加大于前者，其原因是CMFD融入了更多的地面氣象數(shù)據(jù)。

3討 論

3.1總體精度水平對(duì)照

（1） 亞日尺度-旬月尺度的不一致性。

研究基于高時(shí)空分辨率的地面監(jiān)測(cè)數(shù)據(jù)，通過相關(guān)系數(shù)和均方根誤差兩個(gè)指標(biāo)揭示出在亞日尺度MSWEP的精度最優(yōu)，在旬及月尺度則表現(xiàn)為CMFD最優(yōu)。盡管較多研究評(píng)估了不同產(chǎn)品的精度，但3種產(chǎn)品在日、亞日尺度與旬月尺度精度表現(xiàn)的不一致性在前期研究并未被揭示［24］。一是MSWEP的3 h基礎(chǔ)數(shù)據(jù)精度略高于后者，這取決于其主要數(shù)據(jù)源［25］。二是在更大時(shí)間尺度CMFD數(shù)據(jù)融合了更多的中國(guó)大陸區(qū)域氣象觀測(cè)站，從而導(dǎo)致其精度升高更快。

（2） 不同產(chǎn)品間擊中率與相對(duì)誤差中位數(shù)的尺度一致性。

研究揭示出擊中率在所有尺度均表現(xiàn)為MSWEP顯著優(yōu)于其余兩組，CMFD略優(yōu)于TRMM，這表明對(duì)于降水事件是否發(fā)生的識(shí)別，MSWEP的精度明顯優(yōu)于其余兩組數(shù)據(jù)，而TRMM由于融合地面數(shù)據(jù)較少精度表現(xiàn)最差。根據(jù)相對(duì)誤差定義可知，當(dāng)降水產(chǎn)品預(yù)報(bào)降水事件不發(fā)生時(shí)，其相對(duì)誤差為1（100%）；當(dāng)擊中率低于50%時(shí)，降水產(chǎn)品一半以上數(shù)據(jù)的相對(duì)誤差為1，并且在低雨量時(shí)精度較差且易導(dǎo)致相對(duì)誤差大于1，故CMFD和TRMM在亞日尺度的相對(duì)誤差中位數(shù)多為1。

（3） TRMM的時(shí)段選取與精度。

研究表明TRMM由3 h至月尺度的相關(guān)系數(shù)中位數(shù)由0.23逐漸增加至0.90以上（表2），均值在日和月尺度分別增加至0.60和0.70以上（圖3）。相關(guān)研究表明TRMM在日尺度的相關(guān)系數(shù)約為0.50［16-17］，略低于評(píng)價(jià)結(jié)果。其原因一方面是因?yàn)橄嚓P(guān)研究采用的TRMM 3B42 V7 3 h數(shù)據(jù)時(shí)間為世界時(shí)加減90 min，如00：00的3 h數(shù)據(jù)是前一天22：30至后一天01：30的累計(jì)降雨量，24 h數(shù)據(jù)為相鄰22：30至22：30的數(shù)據(jù)［17］。這和中國(guó)氣象局的日值數(shù)據(jù)集采用08：00至08：00的尺度不一致，如中國(guó)地面氣候資料日值數(shù)據(jù)集，故評(píng)比的結(jié)果存在偏差。另一方面是因?yàn)檠芯坎捎玫臍庀笳军c(diǎn)數(shù)據(jù)時(shí)間序列長(zhǎng)度以及來源不同，數(shù)據(jù)精度可能存在一定的隨機(jī)性。研究基于逐小時(shí)數(shù)據(jù)拆分為半小時(shí)數(shù)據(jù)并經(jīng)過數(shù)據(jù)檢驗(yàn)，再對(duì)應(yīng)TRMM數(shù)據(jù)所代表的時(shí)間段，從而更準(zhǔn)確評(píng)價(jià)TRMM數(shù)據(jù)，揭示出TRMM在中國(guó)大陸分區(qū)的精度略高于前期研究結(jié)果。

（4） 數(shù)據(jù)精度的空間分布特征。

研究揭示出3組降水產(chǎn)品的精度由西北至東南逐漸增加，總體表現(xiàn)為濕潤(rùn)地區(qū)最優(yōu)，西北干旱地區(qū)最差。其原因是降水產(chǎn)品對(duì)微量降水事件精度較差，如0.1 mm/3 h，0.1 mm/d等，從而產(chǎn)生上述空間分布特征。需要注意的是，CMFD數(shù)據(jù)在日、旬和月尺度精度最高的區(qū)域?yàn)榍嗖馗咴貐^(qū)，這反映出CMFD在該區(qū)域融合了更多的地面觀測(cè)數(shù)據(jù)，因?yàn)樵摂?shù)據(jù)的開發(fā)團(tuán)隊(duì)一直致力于青藏高原的相關(guān)研究。此外，對(duì)于西北干旱地區(qū)，CMFD產(chǎn)品的精度在日、旬、月的精度也明顯高于其余兩組產(chǎn)品，在月尺度的相關(guān)系數(shù)中位數(shù)達(dá)到0.8，可為研究戈壁沙漠地區(qū)的降水特征提供參考依據(jù)。

3.2尺度變化模型及意義

研究分別基于修正后的對(duì)數(shù)Gamma分布函數(shù)構(gòu)建了相關(guān)系數(shù)和擊中率由小時(shí)至月尺度的演變模型，基于冪函數(shù)構(gòu)建了均方根誤差和相對(duì)誤差中位數(shù)的演變模型，模型參數(shù)可分別指示不同的意義。

3.2.1基于冪函數(shù)的演變模型

將降水產(chǎn)品的值視為期望值，實(shí)測(cè)值視為抽樣值，均方根誤差可視為樣本的標(biāo)準(zhǔn)差。降水序列可采用兩參數(shù)Gamma分布進(jìn)行描述［26-31］，則3 h降水?dāng)?shù)據(jù)的概率分布可表示為服從形狀因子α、尺度因子β的Gamma分布，記為G（α，β）；6 h降水?dāng)?shù)據(jù)為兩個(gè)3 h降水?dāng)?shù)據(jù)之和，假設(shè)二者相互獨(dú)立，則根據(jù)Gamma分布的可加性［32］，其概率分布為G（2α，β）；依次遞推，當(dāng)尺度逐漸增加至3 h的k倍時(shí)，其概率分布為G（kα，β）。而兩參數(shù)Gamma分布的方差為αβ2，標(biāo)準(zhǔn)差為α1/2β，且其樣本的表達(dá)式和均方根誤差形式一致。故當(dāng)各尺度間相互獨(dú)立時(shí)，時(shí)間尺度為1的均方根誤差為標(biāo)準(zhǔn)差α1/2β，時(shí)間尺度為3 h的k倍時(shí)，其誤差為（kα）1/2β。可見，均方根誤差的基本增加趨勢(shì)滿足冪函數(shù)形式，且冪指數(shù)為1/2。因此，當(dāng)增加的冪指數(shù)在1/2附近時(shí)，表明樣本在不同尺度間的獨(dú)立性較好；當(dāng)冪指數(shù)小于1/2時(shí)則有更多的已知信息來減小均方根誤差。研究結(jié)果表明，MSWEP產(chǎn)品均方根誤差隨尺度增加的冪指數(shù)評(píng)價(jià)約為0.5，說明尺度間的相互獨(dú)立性較好，在更大尺度無相關(guān)監(jiān)測(cè)數(shù)據(jù)融合；而TRMM和CMFD則明顯在1 d以上的尺度開始減小，說明融合后的序列將明顯減小更大時(shí)間尺度產(chǎn)品的均方根誤差。

對(duì)于相對(duì)誤差中位數(shù)而言，其變化趨勢(shì)也能用冪函數(shù)進(jìn)行描述，但需要注意在亞日尺度由于CMFD和TRMM的擊中率精度較差，其值多為1，故模型的表現(xiàn)相對(duì)較差。其中TRMM的R2在大部分分區(qū)為0.5～0.7，CMFD除在西北地區(qū)和青藏高原R2達(dá)到0.9外，其余區(qū)域也低于0.9；而MSWEP相對(duì)誤差中位數(shù)隨時(shí)間尺度增加呈冪函數(shù)減小的模型R2普遍高于0.9。因此，當(dāng)降水產(chǎn)品擊中率高時(shí)，相對(duì)誤差中位數(shù)演變趨勢(shì)可視為冪函數(shù)；當(dāng)擊中率較低時(shí)則相對(duì)誤差中位數(shù)多為1，影響模型的精度。

3.2.2基于修正對(duì)數(shù)Gamma分布的演變模型

研究表明，基于修正后的對(duì)數(shù)Gamma分布能夠很好地描述相關(guān)系數(shù)及擊中率隨尺度增加的演變特征，其R2在全國(guó)及各分區(qū)大多達(dá)到0.99。研究模型的4個(gè)參數(shù)能夠很好反映出不同產(chǎn)品的演變特征：其中R0為初始值，表示最小時(shí)間尺度的精度（如相關(guān)系數(shù)和擊中率），1-R0表示變程，即隨著尺度變化的波動(dòng)范圍；形狀因子α決定一階導(dǎo)數(shù)的分形特征，值越大偏態(tài)越??；β表示尺度因子，能夠揭示其演變的拖尾性。當(dāng)前針對(duì)該問題的研究極少，本研究所建立的模型將為這兩個(gè)參數(shù)的分析提供新技術(shù)方法。

4結(jié) 論

（1） 3組降水產(chǎn)品精度對(duì)比方面，擊中率和相對(duì)誤差中位數(shù)在各尺度均表現(xiàn)為MSWEP明顯優(yōu)于CMFD和TRMM，CMFD略優(yōu)于TRMM；但相關(guān)系數(shù)、均方根誤差在亞日尺度和旬月尺度不一致：前者為MSWEP略優(yōu)于CMFD，后者為CMFD略優(yōu)于MSWEP，且二者明顯優(yōu)于TRMM。

（2） 空間分布方面，3組降水產(chǎn)品總體表現(xiàn)為由西北至東南區(qū)域精度逐漸增加的特征，且均對(duì)小降水量事件不夠敏感。但CMFD在西北區(qū)域和青藏高原的精度最高，可能是相關(guān)團(tuán)隊(duì)融入了更多的地面觀測(cè)數(shù)據(jù)。

（3） 降水產(chǎn)品精度隨時(shí)間演變的函數(shù)模型為：相關(guān)系數(shù)和擊中率可用修正后的對(duì)數(shù)Gamma分布進(jìn)行描述，模型參數(shù)能揭示出其最小尺度的初始值、變化范圍、分形特征和拖尾特征；均方根誤差以冪函數(shù)形式增大，冪指數(shù)能反映出降水序列的隨機(jī)性和特定尺度新數(shù)據(jù)源的融入；相對(duì)誤差中位數(shù)以冪函數(shù)形式減小，但TRMM和CFMD在亞日尺度受擊中率低影響而表現(xiàn)較差。

需要注意，上述研究結(jié)果基于2015～2020年的地面監(jiān)測(cè)資料（少量缺失），可能具有一定的隨機(jī)性和偏差；通過對(duì)TRMM數(shù)據(jù)在日及亞日尺度的劃分進(jìn)行修正，其評(píng)估精度略優(yōu)于相關(guān)文獻(xiàn)。這些將為降水產(chǎn)品的優(yōu)選識(shí)別、精度提升和深度融合提供參考依據(jù)。

參考文獻(xiàn)：

［1］TIAN Y，PETERS L C D，EYLANDER J B，et al.Component analysis of errors in satellite‐based precipitation estimates［J］.Journal of Geophysical Research：Atmospheres，2009，114（D24）：D24101.

［2］LI G，YU Z，WANG W，et al.Analysis of the spatial Distribution of precipitation and topography with GPM data in the Tibetan Plateau［J］.Atmospheric Research，2021，247：105259.

［3］朱慧琴，陳生，李曉俞，等.GPM 衛(wèi)星降水產(chǎn)品在 “7· 21” 河南極端暴雨過程中的誤差評(píng)估［J］.水利水電技術(shù)，2022，53（1）：1-13.

［4］TONG K，SU F，YANG D，et al.Tibetan Plateau precipitation as depicted by gauge observations，reanalyses and satellite retrievals［J］.International Journal of Climatology，2014，34（2）：265-285.

［5］MA Y，ZHANG Y，YANG D，et al.Precipitation bias variability versus various gauges under different climatic conditions over the Third Pole Environment （TPE） region［J］.International Journal of Climatology，2015，35（7）：1201-1211.

［6］ALIJANIAN M，RAKHSHANDEHROO G R，MISHRA A K，et al.Evaluationof satellite rainfall climatology using CMORPH，PERSIANN‐CDR，PERSIANN，TRMM，MSWEP over Iran［J］.International Journal of Climatology，2017，37（14）：4896-4914.

［7］LUO X，WU W，HE D，et al.Hydrological simulation using TRMM and CHIRPS precipitation estimates in the lower Lancang-Mekong river basin［J］.Chinese Geographical Science，2019，29：13-25.

［8］ZHU H，LI Y，HUANG Y，et al.Evaluation and hydrological application of satellite-based precipitation datasets in driving hydrological models over the Huifa river basin in Northeast China［J］.Atmospheric Research，2018，207：28-41.

［9］WU W，LI Y，LUO X，et al.Performance evaluation of the CHIRPS precipitation dataset and its utility in drought monitoring over Yunnan Province，China［J］.Geomatics，Natural Hazards and Risk，2019，10（1）：2145-2162.

［10］XIANG Y，CHEN J，LI L，et al.Evaluation of eight global precipitation datasets in hydrological modeling［J］.Remote Sensing，2021，13（14）：2831.

［11］丁光旭，郭家力，湯正陽，等.多種降水再分析數(shù)據(jù)在長(zhǎng)江流域的適用性對(duì)比［J］.人民長(zhǎng)江，2022，53（9）：72-79.

［12］LI X，LI Z，LIN Y.Suitability of TRMM products with different temporal resolution （3-hourly，daily，and monthly） for rainfall erosivity estimation［J］.Remote Sensing，2020，12（23）：3924.

［13］LIU J，XIA J，SHE D，et al.Evaluation of six satellite-based precipitation products and their ability for capturing characteristics of extreme precipitation events over a climate transition area in China［J］.Remote Sensing，2019，11（12）：1477.

［14］QI S，LV A.Applicability analysis of multiple precipitation products in the Qaidam Basin，Northwestern China［J］.Environmental Science and Pollution Research，2022，29：12557-12573.

［15］YANG X，LU Y，TAN M L，et al.Nine-Year systematic evaluation of the GPM and TRMM precipitation products in the Shuaishui river basin in East-Central China［J］.Remote Sensing，2020，12（6）：1042.

［16］MA Q，LI Y，F(xiàn)ENG H，et al.Performance evaluation and correction of precipitation data using the 20-year IMERG and TMPA precipitation products in diverse subregions of China［J］.Atmospheric Research，2021，249：105304.

［17］SHEN Y，XIONG A，WANG Y，et al.Performance of high‐resolution satellite precipitation products over China［J］.Journal of Geophysical Research：Atmospheres，2010，115（D2）：D012097.

［18］HE J，YANG K，TANG W，et al.The first high-resolution meteorological forcing dataset for land process studies over China［J］.Scientific Data，2020，7（1）：25.

［19］張強(qiáng)，趙煜飛，范邵華.中國(guó)國(guó)家級(jí)氣象臺(tái)站小時(shí)降水?dāng)?shù)據(jù)集研制［J］.暴雨災(zāi)害，2016，35（2）：182-186.

［20］MA Y，TANG G，LONG D，et al.Similarity and error intercomparison of the GPM and its predecessor-TRMM multisatellite precipitation analysis using the best available hourly gauge network over the Tibetan Plateau［J］.Remote Sensing，2016，8（7）：569.

［21］REN Z H，ZHAO P，ZHANG Q，et al.Quality control procedures for hourly precipitation data from automatic weather stations in China［J］.Meteorological Monthly，2010，36（7）：123-132.

［22］趙松喬.中國(guó)綜合自然地理區(qū)劃的一個(gè)新方案［J］.地理學(xué)報(bào)，1983，38（1）：1-10.

［23］龔浩哲，侯婷.1961～2021年舟山群島年降水序列變化趨勢(shì)及突變分析［J］.水利水電快報(bào)，2022，43（12）：21-26.

［24］XU Z，WU Z，HE H，et al.Evaluating the accuracy of MSWEP V2.1 and its performance for drought monitoring over mainland China［J］.Atmospheric Research，2019，226：17-31.

［25］BECK H E，WOOD E F，PAN M，et al.MSWEP V2 global 3-hourly 0.1 precipitation：methodology and quantitative assessment［J］.Bulletin of the American Meteorological Society，2019，100（3）：473-500.

［26］岳書旭，胡實(shí)，莫興國(guó)，等.一種基于頻率的GCM日降水偏差校正方法改進(jìn)及其在長(zhǎng)江流域的應(yīng)用［J］.地理研究，2021，40（5）：1432-1444.

［27］BARGER G L，THOM H C.Evaluation of drought hazard［J］.Agronomy Journal，1949，41（11）：519-526.

［28］HANSON L S，VOGEL R.The probability distribution of daily rainfall in the United States［C］∥ World Environmental and Water Resources Congress，Ahupua′A，2008.

［29］NEYMAN J，SCOTT E L.Some outstanding problems relating to rain modification［C］∥ Proceedings of the Fifth Berkeley Symposium on Mathematical Statistics and Probability，Volume 5：Weather Modification.Berkeley：University of California Press，1967.

［30］周向陽，張榮，雷文娟.極端降水事件概率分布識(shí)別方法對(duì)比研究［J］.自然災(zāi)害學(xué)報(bào)，2018，27（5）：158-168.

［31］ZHOU X，LIN J，LIANG X，et al.Rainfall patterns from multiscale sample entropy analysis［J］.Frontiers in Water，2022，4：885456.

［32］黃振平.水文統(tǒng)計(jì)學(xué)［M］.南京：河海大學(xué)出版社，2003.

（編輯：謝玲嫻）