嚴(yán)家寶 賈紹鳳 呂愛(ài)鋒 朱文彬

摘要：利用GRACE重力衛(wèi)星數(shù)據(jù)分析了我國(guó)及十大流域近十年水儲(chǔ)量變化趨勢(shì)、年變化特征、年內(nèi)分布特征以及時(shí)空分布規(guī)律，結(jié)合TRMM降水?dāng)?shù)據(jù)分析了水儲(chǔ)量與降水的關(guān)系。研究結(jié)果表明：近十年，中國(guó)水儲(chǔ)量變化趨勢(shì)具有空間差異性，西南大部、華北平原及黃河中下游、西北準(zhǔn)噶爾盆地一帶水儲(chǔ)量呈減少趨勢(shì)，東南部、長(zhǎng)江大部分區(qū)域、長(zhǎng)江黃河源頭以及塔里木盆地區(qū)域水儲(chǔ)量呈增加趨勢(shì)；中國(guó)水儲(chǔ)量年變化幅度較小，淮河、海河、珠江、松花江流域振幅較大，西北諸河流域振幅最小，全國(guó)除黃河、海河和西北諸河流域外，流域水儲(chǔ)量年變化與降水年變化均呈顯著正相關(guān)，東南諸河、珠江、長(zhǎng)江流域相關(guān)系數(shù)均達(dá)0.7以上；年內(nèi)分布上，我國(guó)冬春季水儲(chǔ)量虧缺，夏秋季水儲(chǔ)量盈余，3月-4月西南諸河及長(zhǎng)江流域水儲(chǔ)量虧缺嚴(yán)重，7月-9月則盈余較大，華北平原5月-7月水儲(chǔ)量有虧缺，其他月份則水儲(chǔ)量略為盈余。此外，黃河、長(zhǎng)江、東南諸河、西南諸河以及珠江流域水儲(chǔ)量與降水量年內(nèi)分布一致性較好，西北諸河流域2月-4月份水儲(chǔ)量與降水一致性較差，其他月份一致性較好，而其它流域則一致性較差。

關(guān)鍵詞：近十年；中國(guó)；十大流域；水儲(chǔ)量；時(shí)空分布

中圖分類號(hào)：TV213 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1672-1683（2016）04-0021-08

Abstract：Water storage change trend，annual change characteristics，intra-year distribution characteristics and spatial-temporal variation characteristics of China and the ten main basins in the last ten years were analyzed using the GRACE （Gravity Recovery and Climate Experiment） data.The relationship between total water storage and precipitation was analyzed by combining the TRMM （Tropical Rainfall Measuring Mission） precipitation data.The results were as follows：In the last ten years，China water storage trend showed some different characteristics.Most of China southwest part，North China Plain and the middle and lower reaches of Yellow River，the Junggar Basin of northwest China all showed a decrease in water storage.Meanwhile，southeast of China，most part of Yangtze River Basin，the origin area of Yellow and Yangtze River and Tarim Basin showed an increase.The annual change range of China water storage was very narrow.The annual change ranges of Huaihe，Haihe，Pearl，Songhuajiang river basins were wide and that of Northwest river basin was the narrowest.The water storage annual change in most of river basins had a significant positive correlation with the precipitation change except Yellow River Basin，Haihe Basin and northwest river basins.The correlation coefficients of southeast river basins，Pearl Basin and Yangtze River Basin were more than 0.7.In intra-annual distribution，China water storage was lack in winter and spring，and adequate in summer and autumn.Southwest river basins and Yangtze River Basin showed a severe deficiency in March and April，and showed more surplus from July to September.North China Plain showed a deficiency from May to July and showed a surplus in other months.Yellow River Basin，Yangtze River Basin，southeast river basins，southwest river basins and Pearl Basin showed a good consistence in intra-annual water storage and precipitation.Northwest river basins showed a good consistence except in the period of February to April.The other basins showed a bad consistence.

Key words：last ten years；China；ten main basins；water storage；spatial-temporal distribution

從垂直分層角度看，區(qū)域水儲(chǔ)量由區(qū)域地表水儲(chǔ)量、冰川積雪水儲(chǔ)量、土壤水儲(chǔ)量和地下水儲(chǔ)量組成[1]。從水量平衡角度看，陸地水儲(chǔ)量的變化是由降水、蒸散發(fā)、徑流等活動(dòng)過(guò)程的綜合反映[2]。水儲(chǔ)量的盈余和虧缺反映了某地水儲(chǔ)量相對(duì)于其長(zhǎng)期水儲(chǔ)量的正負(fù)偏離，分析水儲(chǔ)量的盈虧時(shí)空變化分布規(guī)律，對(duì)了解區(qū)域水平衡和干濕狀況具有重要意義。

在GRACE（Gravity Recovery and Climate Experiment）衛(wèi)星發(fā)射前，研究水儲(chǔ)量變化主要依賴于實(shí)測(cè)數(shù)據(jù)以及數(shù)值計(jì)算模型推算。由于實(shí)測(cè)數(shù)據(jù)的匱乏，以及空間范圍的限制，研究大尺度的水儲(chǔ)量變化受到制約[3]；數(shù)值計(jì)算模型需要提供很多參數(shù)以及輸入數(shù)據(jù)，輸出往往也限制在特定水文要素上或特定深度的水量上[4]，如GLDAS模型輸出不同深度土壤含水量的變化，并不能反映整個(gè)區(qū)域由地表到地下的水儲(chǔ)量變化情況。2002年GRACE衛(wèi)星的發(fā)射，為大尺度水儲(chǔ)量變化研究提供了新途徑[2-3]。研究者利用GRACE衛(wèi)星反演的水儲(chǔ)量進(jìn)行規(guī)律及趨勢(shì)分析[2-3，5–13]、蒸散發(fā)估計(jì)[14-16]、極端天氣分析[17–19]、地下水儲(chǔ)量變化[20-26]等方面的研究。

當(dāng)前，利用GRACE衛(wèi)星對(duì)我國(guó)水儲(chǔ)量相關(guān)方面研究大多集中在流域尺度上。翟寧[27]、許民[28]、Huang[10]等對(duì)長(zhǎng)江流域水儲(chǔ)量變化規(guī)律進(jìn)行了分析總結(jié)；Zhang[29]對(duì)長(zhǎng)江流域干旱事件進(jìn)行分析并與厄爾尼諾現(xiàn)象進(jìn)行相關(guān)分析；許民[14]利用GRACE水儲(chǔ)量數(shù)據(jù)和實(shí)測(cè)數(shù)據(jù)對(duì)黃河源區(qū)的實(shí)際蒸散發(fā)量進(jìn)行估算；尼勝楠[30]結(jié)合水文模型以及降水?dāng)?shù)據(jù)對(duì)長(zhǎng)江和黃河流域水儲(chǔ)量變化規(guī)律進(jìn)行研究；Feng[31]、冉全[21]、任永強(qiáng)[5]、Huang[32]結(jié)合陸面模型以及實(shí)測(cè)數(shù)據(jù)對(duì)海河流域進(jìn)行地下水儲(chǔ)量變化趨勢(shì)及原因進(jìn)行分析，曹艷萍[20]對(duì)黑河流域也做了相似的分析；Wang[33]利用GRACE水儲(chǔ)量變化數(shù)據(jù)對(duì)海河流域近10年發(fā)生的干旱災(zāi)害進(jìn)行研究；Moiwo[34]結(jié)合氣候模型研究喜馬拉雅山及青藏高原地區(qū)的水儲(chǔ)量變化，以及對(duì)水文、生態(tài)的影響；許朋琨[6]對(duì)青藏高原地區(qū)以及雅魯藏布江流域進(jìn)行年、季水儲(chǔ)量變化分析；Long[19]利用GRACE數(shù)據(jù)對(duì)西南洪澇和干旱事件進(jìn)行監(jiān)測(cè)，并利用神經(jīng)網(wǎng)絡(luò)算法對(duì)歷史和未來(lái)進(jìn)行反演和預(yù)測(cè)；李瓊[35]結(jié)合全球水文模型和降水信息對(duì)2010年西南干旱事件進(jìn)行分析。利用GRACE衛(wèi)星數(shù)據(jù)對(duì)全國(guó)水儲(chǔ)量變化進(jìn)行過(guò)研究的學(xué)者有段建賓[36]、邢樂(lè)林[37]、盧飛[38]、Zhao[39]等，這些學(xué)者的研究多偏重于GRACE衛(wèi)星數(shù)據(jù)處理過(guò)程以及典型區(qū)域水儲(chǔ)量變化分析，對(duì)全國(guó)及各流域尺度的詳細(xì)闡釋上相對(duì)缺乏。高歌[40]利用降水和MODIS潛在蒸散發(fā)產(chǎn)品對(duì)全國(guó)十大流域近50年的水分盈虧量進(jìn)行時(shí)空分析，由于潛在蒸散發(fā)與實(shí)際蒸散發(fā)的意義不同，加之MODIS蒸散發(fā)產(chǎn)品在中國(guó)西北大片區(qū)域缺乏覆蓋，所以與中國(guó)實(shí)際水儲(chǔ)量變化情況有差別。本文在借鑒這些研究的基礎(chǔ)上，利用GRACE數(shù)據(jù)產(chǎn)品，對(duì)中國(guó)及十大流域水儲(chǔ)量變化及其時(shí)空規(guī)律進(jìn)行分析，試圖全面了解中國(guó)水儲(chǔ)量變化的趨勢(shì)及原因。

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

1.1 數(shù)據(jù)

1.1.1 GRACE數(shù)據(jù)

GRACE衛(wèi)星任務(wù)是由美德合作的雙星項(xiàng)目，于2002 年發(fā)射升空。GRACE任務(wù)通過(guò)探測(cè)2顆衛(wèi)星之間的距離變化來(lái)反演地球重力場(chǎng)的變化[41]。由時(shí)變重力場(chǎng)可得到地球表面密度變化，一般情況下，直接將地表密度變化轉(zhuǎn)為等效水高變化。

目前，有多個(gè)機(jī)構(gòu)發(fā)布GRACE數(shù)據(jù)產(chǎn)品。本文選用美國(guó)航天局（NASA）噴氣推進(jìn)實(shí)驗(yàn)室（Jet Propulsion Laboratory，JPL）發(fā)布的陸地等效水高格網(wǎng)數(shù)據(jù)月產(chǎn)品數(shù)據(jù)[42]，該數(shù)據(jù)是在官方三個(gè)數(shù)據(jù)處理中心（CSR、GFZ、JPL）發(fā)布的Level-2 RL05球諧系數(shù)產(chǎn)品的基礎(chǔ)上經(jīng)過(guò)一系列處理[43]得到的，數(shù)據(jù)精度為1°×1°。本文選用2002年4月到2015年1月三個(gè)數(shù)據(jù)中心產(chǎn)品各141個(gè)月（部分月份數(shù)據(jù)缺失）的數(shù)據(jù)進(jìn)行研究。

1.1.2 TRMM降水?dāng)?shù)據(jù)

TRMM（Tropical Rainfall Measuring Mission）衛(wèi)星實(shí)驗(yàn)計(jì)劃是美國(guó)和日本合作開(kāi)展的熱帶降雨測(cè)量計(jì)劃，1997年發(fā)射于日本。TRMM衛(wèi)星進(jìn)行全球尺度的熱帶降雨觀測(cè)，是第一次用于定量測(cè)量熱帶降雨的空間衛(wèi)星計(jì)劃，以期更深入的了解熱帶降雨對(duì)全球氣候變化的影響。

TRMM衛(wèi)星提供三個(gè)等級(jí)的資料，在研究中使用的較多的是三級(jí)融合產(chǎn)品3B42和3B43，分別為3小時(shí)降雨數(shù)據(jù)和月降雨數(shù)據(jù)，數(shù)據(jù)空間分辨率為0.25°×0.25°。本文選用2002年4月至2014年12月的TRMM 3B43數(shù)據(jù)。

1.2 研究方法

本文對(duì)GRACE三個(gè)數(shù)據(jù)中心各月數(shù)據(jù)產(chǎn)品先乘以增益系數(shù)，再減去各自2002年4月-2015年1月的平均值，然后對(duì)同月份的三個(gè)中心數(shù)據(jù)進(jìn)行算術(shù)平均，得到141個(gè)月用等效水高表示的水儲(chǔ)量產(chǎn)品，最后對(duì)這141個(gè)月產(chǎn)品進(jìn)行包含缺失值的線性回歸分析，得到2002年4月到2015年1月的線性變化趨勢(shì)。對(duì)各流域進(jìn)行分區(qū)統(tǒng)計(jì)，對(duì)缺失數(shù)據(jù)進(jìn)行三次樣條曲線插值，并對(duì)結(jié)果進(jìn)行顯著性檢驗(yàn)。

在分析年變化特征時(shí)，利用中心點(diǎn)年滑動(dòng)平均算法對(duì)水儲(chǔ)量和降水?dāng)?shù)據(jù)進(jìn)行滑動(dòng)，中心點(diǎn)年滑動(dòng)平均算法公式如下：

TRMM降水?dāng)?shù)據(jù)緯度范圍為50°S～50°N，松花江流域北部部分區(qū)域在數(shù)據(jù)覆蓋范圍之外，筆者對(duì)2003年1月至2014年12月共144個(gè)月的TRMM在松花江流域剩余部分區(qū)域的降水量進(jìn)行統(tǒng)計(jì)，并與對(duì)應(yīng)的GPCP2.2（Global Precipitation Climatology Project）月降水產(chǎn)品在松花江流域的降水進(jìn)行線性回歸分析，得到回歸系數(shù)為1.002，相關(guān)系數(shù)為0.996，說(shuō)明以松花江剩余部分區(qū)域的TRMM產(chǎn)品計(jì)算整個(gè)松花江流域降水是可行的。

2 結(jié)果與分析

2.1 近十年中國(guó)水儲(chǔ)量變化趨勢(shì)

2.1.1 變化趨勢(shì)空間分布特點(diǎn)

2002年4月到2015年1月，全國(guó)水儲(chǔ)量變化趨勢(shì)空間分布表現(xiàn)出明顯的西北-東南走向的條帶分布規(guī)律，空間分布同F(xiàn)eng[31]和Zhao[39]的結(jié)果一致。東北部松花江流域北部大部分區(qū)域，中部從西北塔里木盆地、到可可西里、長(zhǎng)江和黃河源頭、再到我國(guó)東南部大部分區(qū)域，以及最南部云南南部的一小片地區(qū)水儲(chǔ)量增加，其中新疆東南部、阿爾金山自然保護(hù)區(qū)附近、昆侖山沿線區(qū)域以及長(zhǎng)江中下游、岳陽(yáng)-武漢-九江-蕪湖段區(qū)域增幅明顯，水儲(chǔ)量增量在2～3 cm/a；而把三個(gè)水儲(chǔ)量增加區(qū)域隔斷的有兩個(gè)水儲(chǔ)量減少條帶，一是北方從西北諸河流域北部到黃河流域、東部海河流域、遼河流域、淮河北部，二是西南諸河流域，從西藏南部到云貴高原，甚至往東南發(fā)展到海南島，其中喜馬拉雅山沿線與印度、尼泊爾交界區(qū)域水儲(chǔ)量減少量為2～3 cm/a，天山山脈北部、烏魯木齊東部部分區(qū)域水儲(chǔ)量減少量為3～6 cm/a。

水儲(chǔ)量減少原因存在區(qū)域差異。西南諸河流域水儲(chǔ)量減少的主要原因是降水減少[44-45]和冰川質(zhì)量減少[46]；黃淮海區(qū)域是我國(guó)的產(chǎn)糧基地，農(nóng)業(yè)灌溉消耗大量的地表和地下水，致使地下水嚴(yán)重超采[31，47]；新疆西北部水儲(chǔ)量減少的一方面原因?yàn)楸p少[46]，由全球耕地分布圖[50]和地下水灌溉耕地分布圖[49]可以看到，在新疆西北部（灌溉）耕地與水儲(chǔ)量減少區(qū)域基本吻合，筆者認(rèn)為水儲(chǔ)量減少的另一方面原因?yàn)檗r(nóng)業(yè)灌溉以及蒸發(fā)的影響。

2.1.2 全國(guó)及十大流域變化趨勢(shì)

對(duì)2002年4月到2015年1月全國(guó)和十大流域線性變化趨勢(shì)進(jìn)行統(tǒng)計(jì)（表1）。中國(guó)近10年水儲(chǔ)量呈微弱負(fù)增長(zhǎng)，平均年減少速率約為0.65 mm等效水高。各流域上，東南諸河、松花江、長(zhǎng)江以及珠江呈增加趨勢(shì)，平均年增加速率分別為3.22 mm/a、2.74 mm/a、1.96 mm/a和1.72 mm/a，其中東南諸河和松花江增長(zhǎng)趨勢(shì)顯著；西南諸河、海河、淮河、黃河、遼河以及西北諸河均呈顯著減少趨勢(shì)，平均年減少速率分別為6.48 mm/a、6.41 mm/a、4.16 mm/a、3.03 mm/a、1.32 mm/a和0.86 mm/a。

2.2 水儲(chǔ)量年波動(dòng)特征

圖1為中國(guó)及十大流域水儲(chǔ)量年滑動(dòng)和降水距平累積年滑動(dòng)變化曲線。全國(guó)來(lái)看，近十年水儲(chǔ)量變化較為平穩(wěn)，振幅為2.44 cm；2003年初到2004年初，松、遼、西北諸河、黃、淮、海河流域水儲(chǔ)量急劇增加，其中淮、海、松花江流域增幅較大，分別達(dá)到14.56 cm、9.92 cm和7.09 cm；黃、淮、海河流域從2004年到2014年水儲(chǔ)量持續(xù)波動(dòng)減少，其中海河流域和淮河流域在2004年初到2007年初、海河流域在2013年初到2014年、淮河流域在2012年初到2014年減少幅度較大，并在2014年達(dá)到極低值；2002年底到2004年，珠江流域及東南諸河流域水儲(chǔ)量則急劇減少，減少幅度分別為12.49 cm和7.03 cm；松花江流域、遼河流域、西北諸河流域、長(zhǎng)江流域、東南諸河流域在2003年到2014年水儲(chǔ)量在均值（圖中0值）附近上下波動(dòng)，其中松花江流域波動(dòng)幅度最大，達(dá)到11.95 cm，其次為東南諸河，為8.83 cm，波動(dòng)最小的為西北諸河，波動(dòng)幅度為2.33 cm；西南諸河流域近十年水儲(chǔ)量處于下降趨勢(shì)，在2005年初和2008年底有兩個(gè)快速下降階段，2009年末2010年初水儲(chǔ)量達(dá)到極低值，2010年到2014年，水儲(chǔ)量在極低值附近波動(dòng)。

從GRACE水儲(chǔ)量年滑動(dòng)值和降水距平年滑動(dòng)值的線性相關(guān)分析結(jié)果看，中國(guó)東南部相關(guān)系數(shù)較高，東南諸河流域、珠江流域、長(zhǎng)江流域的相關(guān)系數(shù)分別達(dá)到0.75、0.72和0.72，全國(guó)及十大流域除黃河流域、海河流域、西北諸河流域外均通過(guò)0.05的相關(guān)性顯著性檢驗(yàn)。黃河流域、海河流域以及西北諸河流域相關(guān)性較差的原因在本文2.3節(jié)進(jìn)行解釋。

2.3 水儲(chǔ)量年內(nèi)分布特征

2.3.1 水儲(chǔ)量年內(nèi)分布空間變化特征

圖2為2002年到2015年平均各月水儲(chǔ)量空間分布?？傮w上來(lái)看，冬季（12月-2月）和春季（3月-5月）全國(guó)水儲(chǔ)量較為虧缺，夏季（6月-8月）和秋季（9月-11月）水儲(chǔ)量較為盈余。從各月變化看出，12月-次年1月水儲(chǔ)量虧缺范圍覆蓋全國(guó)除黃淮海交匯區(qū)域以外的大部分區(qū)域，其中長(zhǎng)江流域、東南諸河流域、西南諸河流域以及珠江流域虧損較多；2月-3月，西北塔里木盆地及其西北部區(qū)域、東北、華北以及東部沿海水儲(chǔ)量由微弱虧缺轉(zhuǎn)為微弱盈余，西南諸河流域云南范圍虧損加重；4月-6月，黃淮海區(qū)域水儲(chǔ)量由微弱盈余變?yōu)樘潛p，并逐步加重，東南諸河流域水儲(chǔ)量逐漸由虧損變?yōu)橛啵⒅鸩较蜷L(zhǎng)江流域及內(nèi)陸延伸，西南虧損范圍逐步縮小并減輕；7月-9月，全國(guó)大部分區(qū)域水儲(chǔ)量為盈余狀態(tài)，其中長(zhǎng)江流域、東南、西南水儲(chǔ)量盈余量較大，黃淮海區(qū)域水儲(chǔ)量逐步變?yōu)橛?，西北、東北和東南水儲(chǔ)量開(kāi)始轉(zhuǎn)為微弱虧缺； 10月-12月，水儲(chǔ)量虧缺范圍由西北、東北、東南向中國(guó)中部延伸，水儲(chǔ)量虧缺范圍逐步覆蓋全國(guó)除黃淮海交匯區(qū)域以外的大部分區(qū)域，長(zhǎng)江下游虧缺量較重。

2.3.2 全國(guó)及十大流域水儲(chǔ)量年內(nèi)分布特征及其與降水的對(duì)應(yīng)關(guān)系

圖3為2002年到2015年全國(guó)及十大流域水儲(chǔ)量年內(nèi)分布。全國(guó)來(lái)看， 6月-10月水儲(chǔ)量為盈余狀態(tài)，其中8月份水儲(chǔ)量最高，達(dá)到3.03 cm，11月-次年5月，水儲(chǔ)量處于虧缺狀態(tài)，1月出現(xiàn)最小值，為-1.78 cm。流域上，松花江流域2月-5月和8月，水儲(chǔ)量均為盈余狀態(tài)，其余月份均為虧損狀態(tài)，3月份水儲(chǔ)量最多，10月份水儲(chǔ)量最少；遼河流域、海河流域、淮河流域4月-7月份水儲(chǔ)量虧缺，6月虧缺量最大，其余月份基本處于盈余狀態(tài)；西北諸河流域2月-8月份水儲(chǔ)量有盈余，9月-次年1月水儲(chǔ)量為虧缺，7月和11月分別為最大盈余和最大虧缺月份；黃河流域、長(zhǎng)江流域、西南諸河流域、珠江流域、東南諸河流域年內(nèi)分布基本一致，夏秋季水儲(chǔ)量盈余，冬春季水儲(chǔ)量虧缺，中國(guó)東南部盈虧轉(zhuǎn)換時(shí)間較早、中部及南部盈虧變換時(shí)間次之、西北部及西南部盈虧轉(zhuǎn)換時(shí)間較晚。

年內(nèi)振幅上，全國(guó)振幅為4.81 cm，全國(guó)十大流域依次從高到低分別為珠江流域（18.36 cm）、西南諸河流域（17.27 cm）、長(zhǎng)江流域（12.31cm）、淮河流域（9.84 cm）、東南諸河流域（6.41cm）、海河流域（5.97 cm）、黃河流域（4.13 cm）、遼河流域（3.65 cm）、松花江流域（3.49 cm）和西北諸河流域（1.81 cm）。從虧損轉(zhuǎn)為盈余月份上看，松花江流域和遼河流域盈虧月份較為反復(fù)，松花江在2月由虧損轉(zhuǎn)為盈余，8月出現(xiàn)虧損轉(zhuǎn)盈余異常，遼河流域無(wú)明顯盈虧轉(zhuǎn)換月份；其他流域由虧損轉(zhuǎn)為盈余起始月份從早到晚依次為2月（西北諸河流域）、4月（東南諸河流域）、6月（珠江流域、長(zhǎng)江流域）、7月（黃河流域、西南諸河流域）、8月（海河流域、淮河流域）；由盈余轉(zhuǎn)虧損起始月份分別為9月（西北諸河流域）、10月（東南諸河流域）、12月（黃河流域、長(zhǎng)江流域、西南諸河流域、珠江流域）、次年4月（海河流域、淮河流域）。

由圖3中降水與水儲(chǔ)量的對(duì)比可以看出，黃河流域、長(zhǎng)江流域、東南諸河流域、西南諸河流域以及珠江流域水儲(chǔ)量年內(nèi)分布與降水量年內(nèi)分布一致性較好，這些流域除東南諸河流域外，水儲(chǔ)量盈虧狀態(tài)轉(zhuǎn)換均較降水有1月-2月的滯后；松花江流域、遼河流域、海河流域和淮河流域水儲(chǔ)量年內(nèi)分布與降水年內(nèi)分布呈負(fù)相關(guān)，其中海河流域呈顯著負(fù)相關(guān)；西北諸河流域2月-4月份水儲(chǔ)量與降水一致性較差外，其他月份一致性較好。水儲(chǔ)量與降水的正相關(guān)是正常的，因?yàn)樗畠?chǔ)量來(lái)源于降水，降水偏豐應(yīng)該水儲(chǔ)量也偏大。不太正常的是北方松、遼、海、淮流域水儲(chǔ)量與降水呈負(fù)相關(guān)，其原因可能主要有兩個(gè)：一是北方春夏灌溉消耗了大量地表水和地下水，雖然春季降水開(kāi)始增多但也抵償不了灌溉蒸發(fā)損耗，使得降水增加與水儲(chǔ)量減少相對(duì)應(yīng)，這在黃淮海地區(qū)尤其明顯[31，47，50]；二是松花江流域冬春冰凍，蒸發(fā)減少，雖然降水不多，但降水很多以冰雪的形式積蓄起來(lái)，反而水儲(chǔ)量在增加（西北諸河2月-4月份不一致也可能是此原因），松花江流域4月-5月為春耕時(shí)期，水資源消耗大，水儲(chǔ)量開(kāi)始減少，9月-10月為水稻和玉米成熟季，耗水量減少，水儲(chǔ)量開(kāi)始增加。

3 結(jié)論

（1）近十年，中國(guó)水儲(chǔ)量呈微弱減少趨勢(shì)，減少速率為-0.65 mm/a。其空間分布格局呈西北-東南的條帶狀分布，北部松花江流域、中部塔里木木盆地、長(zhǎng)江流域大部、東南諸河和珠江流域、南部云南南部小部分地區(qū)呈增加趨勢(shì)；分隔這三個(gè)水儲(chǔ)量增加區(qū)域的是兩個(gè)水儲(chǔ)量減少條帶，北部從新疆北部一直到海河、淮河流域，南部從西藏南部到云貴高原、海南島，西南諸河、海河、淮河、黃河、遼河和西北諸河流域均呈顯著減少趨勢(shì)。

（2）水儲(chǔ)量年波動(dòng)特征上，中國(guó)水儲(chǔ)量變化較為平穩(wěn)，近十年振幅為2.44 cm?；春恿饔颉⒑：恿饔?、珠江流域、松花江流域振幅較大，振幅均超過(guò)10 cm；西北諸河流域振幅最小，為2.52 cm。由水儲(chǔ)量和年降水的線性相關(guān)分析看出，全國(guó)及十大流域除黃河流域、海河流域、西北諸河流域外均為正相關(guān)，并通過(guò)0.05的相關(guān)性顯著性檢驗(yàn)；東南諸河流域、珠江流域、長(zhǎng)江流域的相關(guān)系數(shù)分別達(dá)到0.75、0.72和0.72；降水在我國(guó)大多流域是引起水儲(chǔ)量變化的重要原因，黃河、海河以及西北諸河流域受其他因素影響，因素可能包括農(nóng)業(yè)灌溉、蒸散發(fā)等。

（3）水儲(chǔ)量年內(nèi)分布上，我國(guó)冬春季水儲(chǔ)量較為虧缺，夏秋季水儲(chǔ)量較為盈余。12月和次年1月，全國(guó)除黃淮海區(qū)域水儲(chǔ)量有盈余外，其他區(qū)域基本處于虧缺狀態(tài)，其中長(zhǎng)江及其南部虧缺較為嚴(yán)重；2月-5月，西北和東北區(qū)域水儲(chǔ)量轉(zhuǎn)為微弱盈余，黃淮海區(qū)域則由盈余轉(zhuǎn)為虧缺；6月-9月，水儲(chǔ)量盈余范圍由東南區(qū)域向內(nèi)陸深入，盈余范圍覆蓋全國(guó)大部分區(qū)域，長(zhǎng)江流域、西南諸河流域、珠江流域及東南諸河流域水儲(chǔ)量盈余量大；10月-12月，水儲(chǔ)量開(kāi)始由盈余變?yōu)樘澣?，范圍逐步由西北、東北、東南向全國(guó)覆蓋。與降水的對(duì)比發(fā)現(xiàn)，黃河流域、長(zhǎng)江流域、東南諸河流域、西南諸河流域以及珠江流域水儲(chǔ)量年內(nèi)分布與降水量年內(nèi)分布一致性較好；松花江流域、遼河流域、海河流域和淮河流域水儲(chǔ)量年內(nèi)分布與降水一致性較差；西北諸河流域2月-4月份水儲(chǔ)量與降水一致性較差外，其他月份一致性較好。

參考文獻(xiàn)（References）：

[1] Li B，Rodell M，Zaitchik B F，et al.Assimilation of GRACE terrestrial water storage into a land surface model：Evaluation and potential value for drought monitoring in western and central Europe[J].Journal of Hydrology，2012，446-447：103-115.

[2] Schmidt R，Schwintzer P，F(xiàn)lechtner F，et al.GRACE observations of changes in continental water storage[J].Global and Planetary Change，2006，50（1-2）：112-126.

[3] Ramillien G，F(xiàn)rappart F，Cazenave A，et al.Time variations of land water storage from an inversion of 2 years of GRACE geoids[J].Earth and Planetary Science Letters，2005，235（1-2）：283-301.

[4] Rodell M，F(xiàn)amiglietti J S.The potential for satellite-based monitoring of groundwater storage changes using GRACE：the High Plains aquifer，Central US[J].Journal of Hydrology，2002，263（1-4）：245-256.

[5] 任永強(qiáng)，潘云，宮輝力.海河流域地下水儲(chǔ)量空間變化趨勢(shì)分析[J].首都師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2014（02）：89–98.（REN Yong-qiang，PAN Yun，GONG Hui-li.Haihe Basin groundwater reserves space trend analysis[J].Journal of Capital Normal University （Natural Science Edition），2014（02）：89-98.（in Chinese） ）

[6] 許朋琨，張萬(wàn)昌.GRACE反演近年青藏高原及雅魯藏布江流域陸地水儲(chǔ)量變化[J].水資源與水工程學(xué)報(bào)，2013（01）：23–29.（XU Peng-kun，ZHANG Wan-chang.Inversion of terrestrial water storage changes in recent years for Qinghai - Tibetan plateau and Yarlung Zangbo River basin by GRACE[J].Journal of Water Resources & Water Engineering，2013（01）：23-29.（in Chinese） ）

[7] Xu M，Ye B，Zhao Q，et al.Estimation of water balance in the source region of the Yellow River based on GRACE satellite data[J].Journal of Arid Land，2013（03）：384-395.

[8] 孫倩，塔西甫拉提·特依拜，丁建麗，等.利用GRACE數(shù)據(jù)監(jiān)測(cè)中亞地區(qū)陸地水儲(chǔ)量動(dòng)態(tài)變化的研究[J].天文學(xué)報(bào)，2014（06）：498-511.（SUN Qian，TASHPOLAT Tiyip，DING Jian-li，et al.GRACE data-based estimation of spatial variations in water storage over the Central Asia during 2003-2013[J].Acta Astronomica Sinica，2014（06）：498–511.（in Chinese） ）

[9] Liu W，Miao Y-X.Analysis of Land Water Storage in Southwest China Based on GRACE Data[J].Chinese Astronomy and Astrophysics，2011，35（4）：421-427.

[10] Huang Y，Salama M S，Krol M S，et al.Analysis of long-term terrestrial water storage variations in the Yangtze River basin[J].Hydrol.Earth Syst.Sci.，2013，17（5）：1985-2000.

[11] Llovel W，Becker M，Cazenave A，et al.Global land water storage change from GRACE over 2002–2009；Inference on sea level[J].Comptes Rendus Geoscience，2010，342（3）：179-188.

[12] Xavier L，Becker M，Cazenave A，et al.Interannual variability in water storage over 2003-2008 in the Amazon Basin from GRACE space gravimetry，in situ river level and precipitation data[J].Remote Sensing of Environment，2010，114（8）：1629-1637.

[13] Singh A，Seitz F，Schwatke C.Inter-annual water storage changes in the Aral Sea from multi-mission satellite altimetry，optical remote sensing，and GRACE satellite gravimetry[J].Remote Sensing of Environment，2012，123：187-195.

[14] 許民，葉柏生，趙求東.基于GRACE重力衛(wèi)星數(shù)據(jù)的黃河源區(qū)實(shí)際蒸發(fā)量估算[J].冰川凍土，2013（01）：138–147.（XU Min，YE Bai-sheng，ZHAO Qiu-dong.Estimation over the real evaporation in the source regions of the Yellow River using GRACE satellite data[J].Journal of Glaciology and Geocryology，2013（01）：138-147.（in Chinese） ）

[15] Long D，Longuevergne L，Scanlon B R.Uncertainty in evapotranspiration from land surface modeling，remote sensing，and GRACE satellites[J].Water Resources Research，2014，50（2）：1131-1151.

[16] Billah M M，Goodall J L，Narayan U，et al.A methodology for evaluating evapotranspiration estimates at the watershed-scale using GRACE[J].Journal of Hydrology，2015，523：574-586.

[17] Morishita Y，Heki K.Characteristic precipitation patterns of El Ni o/La Ni a in time-variable gravity fields by GRACE[J].Earth and Planetary Science Letters，2008，272（3-4）：677-682.

[18] Houborg R，Rodell M，Li B，et al.Drought indicators based on model-assimilated Gravity Recovery and Climate Experiment （GRACE） terrestrial water storage observations[J].Water Resources Research，2012，48（7）：W07525.

[19] Long D，Shen Y，Sun A，et al.Drought and flood monitoring for a large karst plateau in Southwest China using extended GRACE data[J].Remote Sensing of Environment，2014，155：145-160.

[20] 曹艷萍，南卓銅，胡興林.利用GRACE重力衛(wèi)星數(shù)據(jù)反演黑河流域地下水變化[J].冰川凍土，2012（03）：680-689.（CAO Yan-ping，NAN Zhuo-tong，HU Xing-lin.Changes of groundwater storage in the Heihe River Basin derived from GRACE gravity satellite data[J].Journal of Glaciology and Geocryology，2012（03）：680-689.（in Chinese） ）

[21] 冉全，潘云，王一如，等.GRACE衛(wèi)星數(shù)據(jù)在海河流域地下水年開(kāi)采量估算中的應(yīng)用[J].水利水電科技進(jìn)展，2013（02）：42-46，67.（RAN Quan，PAN Yun，WANG Yi-ru，et al.Estimation of annual groundwater exploitation in Haihe River Basin by use of GRACE satellite data[J].Advances in Science and Technology of Water Resources，2013（02）：42-46，67.（in Chinese） ）

[22] Strassberg G，Scanlon B R，Chambers D.Evaluation of groundwater storage monitoring with the GRACE satellite：Case study of the High Plains aquifer，central United States[J].Water Resources Research，2009，45（5）：W05410.

[23] Tiwari V M，Wahr J，Swenson S.Dwindling groundwater resources in northern India，from satellite gravity observations[J].Geophysical Research Letters，2009，36（18）：L18401.

[24] Henry C M，Allen D M，Huang J.Groundwater storage variability and annual recharge using well-hydrograph and GRACE satellite data[J].Hydrogeology Journal，2011，19（4）：741-755.

[25] Scanlon B R，Longuevergne L，Long D.Ground referencing GRACE satellite estimates of groundwater storage changes in the California Central Valley，USA[J].Water Resources Research，2012，48（4）：W04520.

[26] Huang J，Halpenny J，van der Wal W，et al.Detectability of groundwater storage change within the Great Lakes Water Basin using GRACE[J].Journal of Geophysical Research：Solid Earth，2012，117（B8）：B08401.

[27] 翟寧，王澤民，伍岳，等.利用GRACE反演長(zhǎng)江流域水儲(chǔ)量變化[J].武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2009（04）：436–439.（ZHAI Ning，WANG Ze-min，WU Yue，et al.Recovery of Yangtze River Basin water storage variations by GRACE observations[J].Geomatics and Information Science of Wuhan University，2009（04）：436-439.（in Chinese） ）

[28] 許民，葉柏生，趙求東.2002-2010年長(zhǎng)江流域GRACE水儲(chǔ)量時(shí)空變化特征[J].地理科學(xué)進(jìn)展，2013（01）：68–77.（XU Min，YE Bai-sheng，ZHAO Qiu-dong.Temporal and spatial pattern of water storage changes over the Yangtz river basin during 2002-2010 based on GRACE satellite data[J].Progress in Geography，2013（01）：68-77.（in Chinese） ）

[29] Zhang Z，Chao B F，Chen J，et al.Terrestrial water storage anomalies of Yangtze River Basin droughts observed by GRACE and connections with ENSO[J].Global and Planetary Change，2015，126：35-45.

[30] 尼勝楠，陳劍利，李進(jìn)，等.利用GRACE衛(wèi)星時(shí)變重力場(chǎng)監(jiān)測(cè)長(zhǎng)江、黃河流域水儲(chǔ)量變化[J].大地測(cè)量與地球動(dòng)力學(xué)，2014（04）：49-55.（NI Sheng-nan，CHEN Jian-li，LI Jin，et al.Terrestrial water storage change in the Yangtze and Yellow River basins from GRACE time-variable gravity measurements[J].Journal of Geodesy and Geodynamics，2014（04）：49-55.（in Chinese） ）

[31] Feng W，Zhong M，Lemoine J-M，et al.Evaluation of groundwater depletion in North China using the Gravity Recovery and Climate Experiment （GRACE） data and ground-based measurements[J].Water Resources Research，2013，49（4）：2110-2118.

[32] Huang Z，Pan Y，Gong H，et al.Subregional-scale groundwater depletion detected by GRACE for both shallow and deep aquifers in North China Plain[J].Geophysical Research Letters，2015，42（6）：2014GL062498.

[33] Wang J，Jiang D，Huang Y，et al.Drought analysis of the Haihe River Basin based on GRACE terrestrial water storage[J].The Scientific World Journal，2014，2014：e578372.

[34] Moiwo J P，Yang Y，Tao F，et al.Water storage change in the Himalayas from the Gravity Recovery and Climate Experiment （GRACE） and an empirical climate model[J].Water Resources Research，2011，47（7）：W07521.

[35] 李瓊，羅志才，鐘波，等.利用GRACE時(shí)變重力場(chǎng)探測(cè)2010年中國(guó)西南干旱陸地水儲(chǔ)量變化[J].地球物理學(xué)報(bào)，2013（06）：1843–1849.（LI Qiong，LUO Zhi-cai，ZHONG Bo，et al.Terrestrial water storage changes of the 2010 southwest China drought detected by GRACE temporal gravity field[J].Chinese Journal of Geophysics，2013（06）：1843-1849.（in Chinese） ）

[36] 段建賓，鐘敏，閆昊明，等.利用重力衛(wèi)星觀測(cè)資料解算中國(guó)大陸水儲(chǔ)量變化[J].大地測(cè)量與地球動(dòng)力學(xué)，2007（03）：68–71.（DUAN Jian-bin，ZHONG Min，YAN Hao-ming，et al.Recovery of land water storage variations in Chinese mainland by use of GRACE data[J].Journal of Geodesy and Geodynamics，2007（03）：68-71.（in Chinese） ）

[37] 邢樂(lè)林，李輝，劉冬至，等.利用GRACE時(shí)變重力場(chǎng)監(jiān)測(cè)中國(guó)及其周邊地區(qū)的水儲(chǔ)量月變化[J].大地測(cè)量與地球動(dòng)力學(xué)，2007（04）：35-37，56.（XING Le-lin，LI Hui，LIU Dong-zhi，et al.Monthly water distribution in China and its adjacent area from time-variable gravity field of GRACE[J].Journal of Geodesy and Geodynamics，2007（04）：35-37，56.（in Chinese） ）

[38] 盧飛，游為，范東明，等.由GRACE RL05數(shù)據(jù)反演近10年中國(guó)大陸水儲(chǔ)量及海水質(zhì)量變化[J].測(cè)繪學(xué)報(bào)，2015（02）：160-167.（LU Fei，YOU Wei，F(xiàn)AN Dong-ming，et al.Chinese continental water storage and ocean water mass variations analysis in recent ten years based on GRACE RL05 data[J].Acta Geodaetica et Cartographica Sinica，2015（02）：160-167.（in Chinese） ）

[39] Zhao Q，Wu W，Wu Y.Variations in China′s terrestrial water storage over the past decade using GRACE data[J].Geodesy and Geodynamics，2015，6（3）：187-193.

[40] 高歌，許崇育.1961-2010年中國(guó)十大流域水分盈虧量時(shí)空變化特征[J].地理學(xué)報(bào)，2015（03）：380-391.（GAO Ge，XU Chong-yu.Characteristics of water surplus and deficit change in 10 major river basins in China during 1961-2010[J].Acta Geographica Sinica，2015（03）：380-391.（in Chinese） ）

[41] 羅子人，鐘敏，邊星，等.地球重力場(chǎng)空間探測(cè)：回顧與展望[J].力學(xué)進(jìn)展，2014（00）：291-337.（LUO Zi-ren，ZHONG Min，BIAN Xing，et al.Mapping earth's gravity in space：review and future perspective[J].Advances in Mechanics，2014（00）：291-337.（in Chinese） ）

[42] Landerer F W，Swenson S C.Accuracy of scaled GRACE terrestrial water storage estimates[J].Water Resources Research，2012，48（4）：W04531.

[43] Swenson S，Wahr J.Post-processing removal of correlated errors in GRACE data[J].Geophysical Research Letters，2006，33（8）：L08402.

[44] Fan H，Hu J，He D.Trends in precipitation over the low latitude highlands of Yunnan，China[J].Journal of Geographical Sciences，2013，23（6）：1107-1122.

[45] Li Y-G，He D，Hu J-M，et al.Variability of extreme precipitation over Yunnan Province，China 1960-2012[J].International Journal of Climatology，2015，35（2）：245-258.

[46] Matsuo K，Heki K.Time-variable ice loss in Asian high mountains from satellite gravimetry[J].Earth and Planetary Science Letters，2010，290（1-2）：30-36.

[47] Tang Q，Zhang X，Tang Y.Anthropogenic impacts on mass change in North China[J].Geophysical Research Letters，2013，40（15）：3924-3928.

[48] Fritz S，See L，Mccallum I，et al.Mapping global cropland and field size[J].Global Change Biology，2015，21（5）：1980-1992.

[49] Siebert S，Burke J，F(xiàn)aures J M，et al.Groundwater use for irrigation-a global inventory[J].Hydrol Earth Syst Sci，2010，14（10）：1863-1880.

[50] Moiwo J P，Tao F，Lu W.Analysis of satellite-based and in situ hydro-climatic data depicts water storage depletion in North China Region[J].Hydrological Processes，2013，27（7）：1011-1020.