近十年中國陸地水儲量變化及其時空分布規(guī)律

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

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

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

中圖分類號：TV213 文獻標(biāo)志碼：A 文章編號：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ū)域水儲量由區(qū)域地表水儲量、冰川積雪水儲量、土壤水儲量和地下水儲量組成[1]。從水量平衡角度看，陸地水儲量的變化是由降水、蒸散發(fā)、徑流等活動過程的綜合反映[2]。水儲量的盈余和虧缺反映了某地水儲量相對于其長期水儲量的正負(fù)偏離，分析水儲量的盈虧時空變化分布規(guī)律，對了解區(qū)域水平衡和干濕狀況具有重要意義。

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

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

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

1.1 數(shù)據(jù)

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

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

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

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

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

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

1.2 研究方法

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

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

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

2 結(jié)果與分析

2.1 近十年中國水儲量變化趨勢

2.1.1 變化趨勢空間分布特點

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

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

2.1.2 全國及十大流域變化趨勢

對2002年4月到2015年1月全國和十大流域線性變化趨勢進行統(tǒng)計（表1）。中國近10年水儲量呈微弱負(fù)增長，平均年減少速率約為0.65 mm等效水高。各流域上，東南諸河、松花江、長江以及珠江呈增加趨勢，平均年增加速率分別為3.22 mm/a、2.74 mm/a、1.96 mm/a和1.72 mm/a，其中東南諸河和松花江增長趨勢顯著；西南諸河、海河、淮河、黃河、遼河以及西北諸河均呈顯著減少趨勢，平均年減少速率分別為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 水儲量年波動特征

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

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

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

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

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

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

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

年內(nèi)振幅上，全國振幅為4.81 cm，全國十大流域依次從高到低分別為珠江流域（18.36 cm）、西南諸河流域（17.27 cm）、長江流域（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)盈余異常，遼河流域無明顯盈虧轉(zhuǎn)換月份；其他流域由虧損轉(zhuǎn)為盈余起始月份從早到晚依次為2月（西北諸河流域）、4月（東南諸河流域）、6月（珠江流域、長江流域）、7月（黃河流域、西南諸河流域）、8月（海河流域、淮河流域）；由盈余轉(zhuǎn)虧損起始月份分別為9月（西北諸河流域）、10月（東南諸河流域）、12月（黃河流域、長江流域、西南諸河流域、珠江流域）、次年4月（海河流域、淮河流域）。

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

3 結(jié)論

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

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

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

參考文獻（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] 任永強，潘云，宮輝力.海河流域地下水儲量空間變化趨勢分析[J].首都師范大學(xué)學(xué)報（自然科學(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] 許朋琨，張萬昌.GRACE反演近年青藏高原及雅魯藏布江流域陸地水儲量變化[J].水資源與水工程學(xué)報，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)測中亞地區(qū)陸地水儲量動態(tài)變化的研究[J].天文學(xué)報，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ū)實際蒸發(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ù)在海河流域地下水年開采量估算中的應(yīng)用[J].水利水電科技進展，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反演長江流域水儲量變化[J].武漢大學(xué)學(xué)報：信息科學(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年長江流域GRACE水儲量時空變化特征[J].地理科學(xué)進展，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] 尼勝楠，陳劍利，李進，等.利用GRACE衛(wèi)星時變重力場監(jiān)測長江、黃河流域水儲量變化[J].大地測量與地球動力學(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時變重力場探測2010年中國西南干旱陸地水儲量變化[J].地球物理學(xué)報，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)星觀測資料解算中國大陸水儲量變化[J].大地測量與地球動力學(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] 邢樂林，李輝，劉冬至，等.利用GRACE時變重力場監(jiān)測中國及其周邊地區(qū)的水儲量月變化[J].大地測量與地球動力學(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年中國大陸水儲量及海水質(zhì)量變化[J].測繪學(xué)報，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年中國十大流域水分盈虧量時空變化特征[J].地理學(xué)報，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] 羅子人，鐘敏，邊星，等.地球重力場空間探測：回顧與展望[J].力學(xué)進展，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.