王留杰,張行南,2,方園皓,夏達忠,2

(1.河海大學(xué)水文水資源學(xué)院,江蘇南京210098；2.河海大學(xué)水資源高效利用與工程安全國家工程研究中心,江蘇南京210098)

中國區(qū)域地面氣象要素數(shù)據(jù)集在長江上游流域的適用性評估

王留杰1,張行南1,2,方園皓1,夏達忠1,2

(1.河海大學(xué)水文水資源學(xué)院,江蘇南京210098；2.河海大學(xué)水資源高效利用與工程安全國家工程研究中心,江蘇南京210098)

從中國區(qū)域地面要素數(shù)據(jù)集中選取1980年～2010年中國區(qū)域的降水和氣溫兩個要素，與實測降水、氣溫資料進行對比分析,采用確定性系數(shù)、納什系數(shù)、平均誤差和均方根誤差4個指標(biāo)分別從時間變化特征和空間分布特征兩個方面對地面要素集的降水資料和氣溫資料在長江上游流域內(nèi)數(shù)據(jù)質(zhì)量進行系統(tǒng)性評估。結(jié)果表明,氣象網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)總體誤差較小,納什系數(shù)接近于1,擬合度高,具有很好的一致性；氣溫數(shù)據(jù)具有很好的一致性,而降水?dāng)?shù)據(jù)在不同子流域數(shù)據(jù)質(zhì)量不同。

數(shù)據(jù)質(zhì)量評估；中國區(qū)域地面氣象要素數(shù)據(jù)集；中國地面氣候資料日值數(shù)據(jù)集

中國區(qū)域地面氣象要素數(shù)據(jù)集(China Meteorolo-gical Forcing Dataset,CMFD)是中國科學(xué)院青藏高原研究所開發(fā)的一套近地面氣象與環(huán)境要素再分析數(shù)據(jù)集。該數(shù)據(jù)集是以國際上現(xiàn)有的Princeton再分析資料、GLDAS資料、GEWEX-SRB輻射資料以及TRMM降水資料為背景場,融合了中國氣象局常規(guī)氣象觀測數(shù)據(jù)制作而成[1-3]。目前,已有學(xué)者用該數(shù)據(jù)針對不同流域、不同研究和不同的氣象要素進行應(yīng)用與評估。如：Bao-Lin Xue等[4]人發(fā)現(xiàn)CMFD可很好地模擬出青藏高原的那曲河流域夏季能量通量和地表溫度等。闞寶云等[5]分析了2003年～2009年基于衛(wèi)星觀測的降水?dāng)?shù)據(jù)CMORPH、TMPA、3B42V6、中國科學(xué)院青藏高原研究所的融合數(shù)據(jù)ITPCAS和基于地面臺站的APHRODIT(2003年～2007年)4套降水?dāng)?shù)據(jù)集在葉爾羌河上游流域的時空分布特征,空間分布上,4套降水資料在研究流域的差異較大；時間序列上,4套降水資料與流域站點降水存在著不同程度的差異。趙捷等[6]基于站點氣象數(shù)據(jù)及CMFD利用Spline插值方法對黑河流域潛在蒸發(fā)量(PET)進行估算發(fā)現(xiàn)Spline空間插值獲得的PET空間分布與柵格數(shù)據(jù)計算獲得的結(jié)果基本一致。Donglin Guo等[7]模擬了青藏高原的永久凍土和季節(jié)性凍土的狀況,模擬的趨勢是接近地表凍土區(qū)域的,只是模擬值稍微偏大。Chen,Y等[8]利用CMFD為Noah驅(qū)動模擬中國干旱半干旱地區(qū)的地表面溫度,其平均誤差比用GLDAS數(shù)據(jù)作為驅(qū)動明顯減小,取得了滿意的模擬結(jié)果。從目前CMFD的應(yīng)用結(jié)果來看,缺少對CMFD降水和氣溫數(shù)據(jù)在長江上游流域的質(zhì)量評估,而降水和氣溫是陸面模式的重要驅(qū)動輸入。鑒此,本文選取1980年～2010年的中國區(qū)域地面要素數(shù)據(jù)集中的降水和氣溫兩個要素與實測資料,采用確定性系數(shù)、平均誤差和均方根誤差等指標(biāo),分別從時間變化特征和空間分布特征兩個方面對地面要素集降水資料和氣溫資料在長江上游流域內(nèi)數(shù)據(jù)質(zhì)量進行系統(tǒng)性評估。

1 資料和方法

1.1 資料

中國區(qū)域?qū)崪y站點數(shù)據(jù)是由國家氣象信息中心提供的1951年～2010年的中國地面氣候資料日值數(shù)據(jù)集(V3.0)(obs)。該數(shù)據(jù)集包含了824個基準(zhǔn)、基本氣象站,1951年1月以來該站氣壓、氣溫、降水量、蒸發(fā)量、相對濕度、風(fēng)向風(fēng)速、日照時數(shù)和0 cm地溫要素集。中國地面氣候資料是日值數(shù)據(jù),選取長江上游77個氣象站1980年～2010年的逐日平均氣溫及降水量數(shù)據(jù),以此為基礎(chǔ)整合合成月尺度降水?dāng)?shù)據(jù)和氣溫數(shù)據(jù)。

中國區(qū)域地面氣象要素數(shù)據(jù)集時間跨度為1980年～2010年,其時間分辨率為 3 h,水平空間分辨率 0.1°,水平空間覆蓋范圍70°E-140°E、15°N-55°N中國陸地區(qū)域,包含近地面氣溫、近地面氣壓、近地面空氣比濕、近地面全風(fēng)速、地面向下短波輻射、地面向下長波輻射、地面降水率,共 7 個要素(變量)[9-11]。中國區(qū)域地面氣象要素數(shù)據(jù)集是按3 h時間尺度生成各氣象要素,本文在此基礎(chǔ)合成月尺度降水?dāng)?shù)據(jù)和氣溫數(shù)據(jù),并篩選出與站點相對應(yīng)的77個網(wǎng)格數(shù)據(jù)。

1.2 評估方法

根據(jù)地面觀測站點與氣象數(shù)據(jù)網(wǎng)格的空間關(guān)系篩選出對應(yīng)的網(wǎng)格。以月為尺度進行數(shù)據(jù)的對比驗證,采用的指標(biāo)主要包括變化趨勢、擬合優(yōu)度即確定系數(shù)(R2)來表示相關(guān)性、平均偏差(EAVE)、納什系數(shù)(ENS)和均方根誤差(ERMS)等。即

(1)

(2)

(3)

(4)

2 區(qū)域月尺度數(shù)據(jù)變化趨勢

2.1 變化過程時段統(tǒng)計特征

在時間尺度上，計算了實測降水與網(wǎng)格降水在整個長江上游流域的月降水量、實測氣溫與網(wǎng)格氣溫在整個長江上游流域的逐月平均值,并繪于圖1和圖2。從圖1中可以看出,1980年～2010年網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)總體具有很好的一致性；只有1986年7月和1996年7月2個月中,網(wǎng)格數(shù)據(jù)略大于實測數(shù)據(jù),分別大了11.7%和25.53%。從圖2可以看出,網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)具有非常好的一致性,在月尺度上網(wǎng)格數(shù)據(jù)比實測數(shù)據(jù)偏小,偏差在10.9%。

圖1 長江上游流域網(wǎng)格與實測月降水量變化

圖2 長江上游流域網(wǎng)格與實測逐月平均氣溫變化

2.2 趨勢性檢驗

趨勢性檢驗的方法采用seasonal-kendall(SK)檢驗法[12],SK檢驗是一種僅考慮數(shù)據(jù)相對排列的非參數(shù)檢驗方法,該檢驗可用于資料系列存在漏測值、未檢出值以及變量分布與正態(tài)分布無關(guān)的時間序列,季節(jié)上可為個？月[13]。SK趨勢檢驗輸出的結(jié)果主要包括變化率τ值和顯著性水平p值,τ值為正代表上升,為負代表下降0.05

2.2.1 網(wǎng)格數(shù)據(jù)月降水變化過程的趨勢性檢驗

基于長江上游流域的11個子流域,計算了1980年～2010年的網(wǎng)格與實測月降水?dāng)?shù)據(jù)在各個子流域的平均月降水,對各子流域平均月降水序列的趨勢性進行了檢驗(見表1)。

表1 各子流域網(wǎng)格和實測降水?dāng)?shù)據(jù)的SK趨勢檢驗結(jié)果

注：τ>0代表上升,τ<0代表下降；**表示在0.05顯著性水平下存在趨勢；*表示在0.1顯著性水平下存在趨勢。

由表1可看出,實測數(shù)據(jù)只有雅礱江流域的降水表現(xiàn)為顯著性增加,其他子流域無顯著性變化；網(wǎng)格數(shù)據(jù)顯在大渡河和沱江流域降水表現(xiàn)為顯著性減少,其他子流域無顯著性變化。

2.2.2 網(wǎng)格數(shù)據(jù)月氣溫變化過程的趨勢性檢驗

采用與降水變化相同的趨勢性檢驗方法,同樣基于長江上游流域的11個子流域,根據(jù)1980年～2010年網(wǎng)格與實測月氣溫數(shù)據(jù)計算了各個子流域的平均月降水,并對其序列的趨勢性進行了檢驗,結(jié)果見表2。由表2可見,整個長江上游流域氣溫都有明顯的上升趨勢；網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)一致,也表現(xiàn)出在整個長江上游流域氣溫有明顯的上升趨勢[15]。

3 網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)的相關(guān)分析

在空間尺度上從整個長江流域上游與子流域兩個角度進行網(wǎng)格與實測數(shù)據(jù)的相關(guān)分析,逐月計算網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)在流域上的確定系數(shù)R2(見圖3)。

表2 各子流域網(wǎng)格和實測氣溫數(shù)據(jù)的SK趨勢檢驗結(jié)果

圖3 長江上游流域確定性系數(shù)

3.1 降水?dāng)?shù)據(jù)的相關(guān)分析

從圖3可見,1980年～2010年網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)月降水具有很高的擬合度,尤其是中低位的數(shù)據(jù),確定性系數(shù)近似為1。表3顯示,11個子流域的確定性系數(shù)都大于0.8,且除沱江流域外的其他流域都在0.9以上,更有大部分子流域的確定性系數(shù)在0.95以上。

表3 降水資料質(zhì)量評估

與實測數(shù)據(jù)相關(guān)性的時間變化特征。計算網(wǎng)格數(shù)據(jù)和實測降水資料在長江上游流域的確定系數(shù)的逐月序列,并對確定系數(shù)時間系列進行了趨勢分析(見圖4和表3、表4)。結(jié)果表明：對于整個長江上游流域,網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)具有很好的相關(guān)性；網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)的相關(guān)性具有季節(jié)性,1月份和12月份相對較低,10月份相對較高?？傮w上,確定系數(shù)均在0.85以上。

圖4 不同月份降水、氣溫數(shù)據(jù)與實測數(shù)據(jù)在長江上游流域的確定性系數(shù)

表4 網(wǎng)格數(shù)據(jù)與實測降水?dāng)?shù)據(jù)確定系數(shù)時間序列SK檢驗

3.2 氣溫數(shù)據(jù)的相關(guān)分析

首先計算實測氣溫與網(wǎng)格氣溫在長江上游流域的逐月確定系數(shù),并畫出其逐月變化曲線[16](見圖5)。

圖5 長江上游流域網(wǎng)格與實測氣溫數(shù)據(jù)逐月確定系數(shù)

從圖5可以看出,1980年～2010年網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)總體具有較好的一致性,相關(guān)系數(shù)均在0.85以上。除了1996年1月份的相關(guān)系數(shù)為0.86,其余月份的相關(guān)系數(shù)都在0.9以上。根據(jù)表5可以看到,11個子流域中,每個子流域的確定性系數(shù)高達0.99,更有部分流域確定性系數(shù)接近于1。同樣,對相關(guān)系數(shù)時間系列進行了趨勢分析,發(fā)現(xiàn)溫度的確定系數(shù)基本上不隨月份變化,不具有明顯的季節(jié)性,全年1月～12月,確定系數(shù)均在0.95以上,全年都具有很強的相關(guān)性(見圖3)。

表5 氣溫資料質(zhì)量評估

4 誤差分析

4.1 降水?dāng)?shù)據(jù)的誤差分析

以實測降水?dāng)?shù)據(jù)為基準(zhǔn),采用EAVE、ENS和ERMS3個指標(biāo),對網(wǎng)格數(shù)據(jù)進行誤差分析,主要對子流域部分進行誤差分析(見表3)。逐月計算網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)在長江上游流域所有站點的降水?dāng)?shù)據(jù)的誤差指標(biāo),同時也計算了誤差指標(biāo)在長江上游流域的逐月變化序列在各月份的均值(見圖6)。分析計算結(jié)果表明：

圖6 網(wǎng)格降水?dāng)?shù)據(jù)EAVE和ERMS的季節(jié)變化

(1)網(wǎng)格數(shù)據(jù)的EAVE的絕對值在大部分流域都在1左右,與實測數(shù)據(jù)偏差不大,除了嘉陵江、金沙江下游和長江上游干流流域以外其他的九個流域的EAVE均為負值。EAVE具有明顯的季節(jié)性,由于受降水的影響,呈現(xiàn)以年為周期的變化。

(2)網(wǎng)格數(shù)據(jù)的ENS均在0.75以上,沱江流域的ENS為0.772,大渡河流域的ENS為0.806以外,其他九個子流域的納什效率系數(shù)均在0.95以上。納什系數(shù)在年內(nèi)分布均在0.9左右,無明顯季節(jié)性差異。

(3)網(wǎng)格數(shù)據(jù)的ERMS,大渡河和沱江流域的明顯偏大都在30以上,其他流域的ERMS都在15左右,在嘉陵江、金沙江上游和雅礱江流域ERMS只有7左右。根據(jù)SK檢驗,ERMS具有明顯的季節(jié)性,受降水季節(jié)的影響,夏季偏差較大,冬季偏差較小,具有非常明顯的季節(jié)性。

4.2 氣溫數(shù)據(jù)的誤差分析

采用與降水誤差相同的分析方法,同樣以實測氣溫數(shù)據(jù)為基準(zhǔn),采用EAVE、ENS和ERMS3個指標(biāo),對網(wǎng)格數(shù)據(jù)進行誤差分析,主要從子流域部分進行誤差分析(見表5)。逐月計算網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)在長江上游流域所有站點的氣溫數(shù)據(jù)的誤差指標(biāo),同時也計算了誤差指標(biāo)在長江上游流域的逐月變化序列在各月份的均值。結(jié)果表明：

(1)網(wǎng)格數(shù)據(jù)的EAVE的絕對值在大部分流域都在1左右,最大的偏差是涪江流域為-6.52,與實測數(shù)據(jù)偏差不大,除了赤水河和沱江流域以外其他的九個流域的EAVE均為負值。EAVE年內(nèi)分布比較均勻,不具有明顯的季節(jié)性變化。

(2)網(wǎng)格數(shù)據(jù)的ENS在每一個子流域中都為0.999,并且在年內(nèi)分布均在0.9左右,無明顯季節(jié)性差異。

(3)網(wǎng)格數(shù)據(jù)的ERMS,在涪江與沱江流域相對較大,分別為6.56和7.87,赤水河流域ERMS值只有0.54,并且ERMS分布并無明顯的季節(jié)性變化。

5 對網(wǎng)格降水與氣溫數(shù)據(jù)質(zhì)量的總體評價

從整個長江上游流域逐月平均值來看(見圖1),網(wǎng)格月降水?dāng)?shù)據(jù)與實測數(shù)據(jù)在1980年～2010年均具有較好的一致性；在1980年～2010年,網(wǎng)格月氣溫數(shù)據(jù)與實測數(shù)據(jù)時間變化上具有高度一致性(見圖2)。從相關(guān)性看,網(wǎng)格與實測的月降水的確定性系數(shù)均大于0.8,具有較好的一致性；除沱江流域以外,其他子流域的確定性系數(shù),均大于0.9；而網(wǎng)格月氣溫數(shù)據(jù)的確定性系數(shù)均在0.99以上。從各項誤差指標(biāo)來看,網(wǎng)格數(shù)據(jù)的EAVE的絕對值在大部分流域都在1左右,與實測數(shù)據(jù)相差不大,降水?dāng)?shù)據(jù)的ENS均在0.75以上,氣溫數(shù)據(jù)的ENS在每一個子流域中都為0.999,降水和氣溫均具有很高的擬合度,對于降水?dāng)?shù)據(jù)的ERMS,在大渡河和沱江流域明顯偏大都在30以上,其他流域的ERMS都在15左右,誤差在接受范圍,氣溫數(shù)據(jù)的ERMS最大只有7.87。

因此,對于整個長江上游流域以及每一個子流域,網(wǎng)格數(shù)據(jù)與實測數(shù)據(jù)具有很高的一致性,能夠真實地反映長江上游流域降水和氣溫的時空分布與變化。

[1]MA Yaoming, SU Zhongbo, KOIKE T, et al. On measuring and remote sensing surface energy partitioning over the Tibetan Plateau—from GAME/Tibet to CAMP/Tibet[J]. Physics & Chemistry of the Earth Parts A/B/C, 2003, 28(1- 3): 63- 74．

[2]ABATZOGLOU J T. Development of gridded surface meteorological data for ecological applications and modelling[J]. International Journal of Climatology, 2013, 33(1)： 121- 131.

[3]SHEFFIELD J, GOTETI G, WOOD E F. Development of a 50-Year High-Resolution Global Dataset of Meteorological Forcings for Land Surface Modeling[J]. Journal of Climate, 2006, 19(13)：3088- 3111.

[4]XUE B L, WANG L, YANG K, et al. Modeling the land surface water and energy cycles of a mesoscale watershed in the central Tibetan Plateau during summer with a distributed hydrological model[J]. Journal of Geophysical Research Atmospheres, 2013, 118(16): 8857- 8868．

[5]闞寶云, 蘇鳳閣, 童凱, 等. 四套降水資料在喀喇昆侖山葉爾羌河上游流域的適用性分析[J]. 冰川凍土, 2013, 35(3): 710- 722．

[6]趙捷, 徐宗學(xué), 左德鵬. 黑河流域潛在蒸散發(fā)量時空變化特征分析[J]. 北京師范大學(xué)學(xué)報： 自然科學(xué)版, 2013(Z1)： 164- 169．

[7]GUO Donglin, WANG Huijun. Simulation of permafrost and seasonally frozen ground conditions on the Tibetan Plateau, 1981—2010[J]. Journal of Geophysical Research: Atmospheres, 2013, 118(11)： 5216- 5230.

[8]CHEN Y, YANG K, HE J, et al. Improving land surface temperature modeling for dry land of China[J]. Journal of Geophysical Research Atmospheres, 2011, 116(D20)： 999- 1010．

[9]YANG K, HE J, TANG W J, et al. On downward shortwave and longwave radiations over high altitude regions: Observation and modeling in the Tibetan Plateau[J]. Agricultural and Forest Meteorology, 2010, 150(1): 38- 46．

[10]HE J. Development of surface meteorological dataset of China with high temporal and spatial resolution (in Chinese)[D]. Institute of Tibetan Plateau Research, Chinese Academy of Sciences, 2010．

[11]何杰, 陽坤. 中國區(qū)域高時空分辨率地面氣象要素驅(qū)動數(shù)據(jù)集[DB/OL]. 中國科學(xué)院青藏高原研究所,寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心, 2011. http:∥westdc.westgis.ac.cn/data/7a35329c-c53f- 4267-aa07-e0037d913a21

[12]HIRSCH R M, SLACK J R, SMITH R A. Techniques of Trend Analysis for Monthly Water Quality Data[J]. Water Resources Research, 1982, 18(1)： 107- 121．

[13]SMITH R A, HIRSCH R M, SLACK J R. A study of trends in total phosphorus measurements at NASQAN stations[R]. United States Geological Survey, Water-supply Paper 2190, 1982．

[14]王文, 汪小菊, 王鵬. GLDAS月降水?dāng)?shù)據(jù)在中國區(qū)的適用性評估[J]. 水科學(xué)進展, 2014(6): 769- 778．

[15]王順久. 長江上游川江段氣溫、 降水及徑流變化趨勢分析[J]. 資源科學(xué), 2009, 31(7)： 1142- 1149．

[16]汪小菊, 王文, 王鵬. GLDAS月氣溫數(shù)據(jù)在中國區(qū)的適用性評估[J]. 水電能源科學(xué), 2014(11): 10- 13, 57．

(責(zé)任編輯 陳 萍)

Applicability Assessment of China Meteorological Forcing Dataset in Upper Yangtze River Basin

WANG Liujie1, ZHANG Xingnan1,2, FANG Yuanhao1, XIA Dazhong1,2

(1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, Jiangsu, China;2. National Engineering Research Center of Water Resources Efficient Utilization and Engineering Safety,Hohai University, Nanjing 210098, Jiangsu, China)

To assess the applicability of China Meteorological Forcing Dataset (CMFD) in upper Yangtze River Basin, the time variation and spatial distribution of precipitation and temperature between 1980 and 2010 are chosen and then compared with Chinese sites-based observations in terms of goodness of fit, Nash-Sutcliffe coefficient, mean bias error and root mean square error. The results show that the precipitation and temperature of CMFD are generally good consistency with those of observations in terms of a lower bias, a higher NS coefficient and a higher goodness of fit. In addition, the temperature data show a good consistency while the quality of CMFD precipitation data is different in different sub-basins．

data quality assessment; China Meteorological Forcing Dataset; Chinese Surface Climate Dataset

2016- 08- 10

國家自然科學(xué)基金資助項目(51420105014)；中國長江三峽集團公司資助項目(0799556)

王留杰(1991—),男,河南鄭州人,碩士研究生,主要從事水文物理規(guī)律模擬及水文預(yù)報研究．

P338.8；P468.4

A

0559- 9342(2017)03- 0018- 05