曹 陽，楊明祥，雷曉輝，張 靜

(1. 城市環(huán)境過程與數(shù)字模擬國家實驗室基地，首都師范大學，北京 100048；2. 三維信息獲取與應用教育部重點實驗室，首都師范大學，北京100048；3. 中國水利水電科學研究院水資源研究所，北京 100038)

0 引 言

水是國民經(jīng)濟的命脈，是人類賴以生存和發(fā)展的基礎，充足的水資源對于維持我國經(jīng)濟長期健康穩(wěn)定發(fā)展尤為重要，因此探究我國不同地區(qū)水資源現(xiàn)狀和變化趨勢，有利于水資源的優(yōu)化利用，對于緩解我國水資源緊缺現(xiàn)狀、提高生態(tài)環(huán)境質(zhì)量、促進社會經(jīng)濟的繁榮，具有非常重要的戰(zhàn)略意義。研究陸地水儲量的傳統(tǒng)方法主要有：①通過遙感衛(wèi)星進行反演：這種方法局限于只能獲取地表十余厘米深度土壤含水量的變化信息；②地面站點觀測：該方法缺點在于人力、物力投入大且觀測范圍小，在地形復雜區(qū)域難以進行工作；③水文建模模擬：由于模型自身的不確定性以及有些地區(qū)觀測數(shù)據(jù)缺失等因素，會致使模型模擬結(jié)果的精度存在較大的偏差。由此可見傳統(tǒng)水資源監(jiān)測方法對區(qū)域陸地水儲量研究的弊端越來越突出。通過重力衛(wèi)星則可以提取區(qū)域水儲量變化信息，該法不僅有效避免了傳統(tǒng)方法的缺點，也為對全球范圍內(nèi)的陸地水儲量的研究提供了嶄新的思路。

GRACE重力位模型通過扣除潮汐(海洋潮汐、固體潮汐和地球自轉(zhuǎn)帶來的極潮)和非潮汐(大氣和海洋)變化造成的影響，剩余的信號則主要表示地球表面非大氣、非海洋的變化信息，在陸地地區(qū)則主要表示陸地水儲量的變化情況[1]。因此，自2004年GRACE衛(wèi)星時變重力場數(shù)據(jù)公布以來，可利用其對區(qū)域陸地水儲量變化和水循環(huán)變化規(guī)律在月尺度上進行研究分析，這促進了地球物質(zhì)遷移和水文領域研究的蓬勃發(fā)展，同時也開辟了陸地水儲量變化研究的新思路。本文綜述了GRACE衛(wèi)星在國內(nèi)陸地水儲量的應用現(xiàn)狀和展望，為GRACE衛(wèi)星在國內(nèi)水文領域的研究提供參考。

1 GRACE重力衛(wèi)星

從Baker[2]在1960年首先提出可以利用衛(wèi)星與衛(wèi)星間的跟蹤模式來觀測地球重力場變化的設想以來，大量學者和研究團體在重力衛(wèi)星探測地球重力場變化的研究中獲得了重要進展(如表1所示)。

表1 衛(wèi)星重力測量理論研究和技術發(fā)展進程[6]Tab.1 Process of theory study and technology developmentof satellite gravity measurement[6]

GRACE是由美國和德國共同研制并于2002年發(fā)射升空的重力衛(wèi)星，其設計目的主要是通過得到高精度的地球時變重力場，為研究地球質(zhì)量變化提供有效觀測數(shù)據(jù)[3,4]。該計劃最初的構(gòu)想是在大地測量學方向得到廣泛運用，但在不斷的探索中，學者發(fā)現(xiàn)GRACE數(shù)據(jù)在氣候?qū)W和水文學等領域也能得到很好的運用。目前GRACE應用領域包括地球物理學、海洋學、冰川學和水文學等相關領域，并在水文學領域?qū)崿F(xiàn)了初次使用空間探測來研究區(qū)域范圍內(nèi)的陸地水儲量(Terrestrial Water Storage, TWS)的突破[5]。

GRACE由兩顆近極圓軌道運行的低空衛(wèi)星(GRACE A和GRACE B)組成，兩顆衛(wèi)星之間的距離保持在170～270 km之間。由于地球系統(tǒng)物質(zhì)在不斷的遷移，導致其質(zhì)量也在不斷的重新分配，這使得地球平均重力場處在變化之中，主要因素則來自于水的質(zhì)量轉(zhuǎn)移。因此，地球重力場的局部變動會導致GRACE兩顆子衛(wèi)星二次過境時的飛行姿態(tài)發(fā)生變化，導致兩顆子衛(wèi)星之間的相對位置發(fā)生變動。這一微小的距離變化會被重力衛(wèi)星上的K波段測距系統(tǒng)所測得并將這種變化表現(xiàn)在它的數(shù)據(jù)序列中。得到的數(shù)據(jù)序列既包含了地球重力場成分，也包含了其他非保守力作用信息。后者則需要通過GPS記錄和非保守力加速度儀數(shù)據(jù)計算扣除，剩余部分可被用作估測地球重力場的變動。

2 GRACE重力衛(wèi)星反演陸地水儲量原理

利用解算好的GRACE月重力場模型數(shù)據(jù)求解，其數(shù)學公式如下：

(1)

由于GRACE重力場系數(shù)的誤差隨l的增大而增大，因此計算過程中還得考慮高階項造成的誤差。為降低高階項誤差對計算結(jié)果造成的影響，Jekeli[7]提出了一種利用平滑核函數(shù)計算表面區(qū)域平均變化的新方法，在公式(2)中引入高斯(Gaussian)平滑核函數(shù)W(a),即：

Wt[ΔCmlcos(mφ)+ΔSmlsin(mφ)]

(2)

3 GRACE在國內(nèi)陸地水儲量中的研究現(xiàn)狀

GRACE通過獲取區(qū)域水儲量的變化信息為系統(tǒng)評價水資源變化提供重要科學依據(jù)。Zhao等[8]利用GRACE衛(wèi)星2003-2013年數(shù)據(jù)對我國陸地水儲總量的變化情況做了研究總結(jié)，并將我國水儲量變化較大的區(qū)域劃分成5個區(qū)(A-E)。山西省和新疆天山地區(qū)水儲量有明顯的下降趨勢，其原因主要是因為人類活動導致了該地區(qū)的水儲量的減少。相比之下，新疆、西藏和青海省的交界地區(qū)由于生態(tài)環(huán)境的提高和缺少人類活動的影響，水儲量則表現(xiàn)出增長的趨勢。同時，重慶、貴州、湖南三省交界區(qū)域和江西、浙江、福建三省交界區(qū)域的陸地水儲量也出現(xiàn)增加態(tài)勢，具體原因是由于這些地區(qū)年降水量豐富，而且沒有明顯的地下水資源量消耗。GRACE重力衛(wèi)星受其空間分辨率(20 萬km2)的限制，其數(shù)據(jù)多應用于大尺度的區(qū)域研究。下文通過敘述GRACE在我國河流流域、湖泊流域和高原平原地區(qū)等區(qū)域水儲量的研究來說明GRACE的應用領域和研究現(xiàn)狀。

3.1 GRACE在河流流域尺度的研究應用

相對于區(qū)域水儲量的研究，流域水儲量涉及多種水量交換而更復雜。目前GRACE衛(wèi)星進行國內(nèi)河流流域的水儲量的研究還主要集中在長江、黃河、海河和黑河等大江大河上，由于GRACE衛(wèi)星估測水儲量存在一定的誤差，因此很多學者將GRACE衛(wèi)星數(shù)據(jù)和水文模型數(shù)據(jù)結(jié)合作對比分析。許民等[9]通過高斯平滑濾波對GRACE觀測到的時變重力場數(shù)據(jù)進行反演后得出長江流域的陸地水儲量變化情況。結(jié)果表明長江流域的水儲量在年際上表現(xiàn)出上升態(tài)勢，并出現(xiàn)顯著的季節(jié)性變化特征。同時對比GRACE反演結(jié)果與GLDAS(Global Land Data Assimilation System)模型結(jié)果，兩者的吻合度較高。Luo等[10]結(jié)合GRACE數(shù)據(jù)估算了珠江流域的水儲量變化，并與WGHM(WaterGAP Global Hydrology Model)和GLDAS水文模型計算結(jié)果比較，最后表明兩個水文模型模擬的珠江水儲量變化證實GRACE反演的珠江流域水儲量的季節(jié)和年際變化特征的準確性。Moiwo等[11]利用GRACE時變重力場數(shù)據(jù)研究了海河流域的水儲量變動情況，并將估算結(jié)果與實測水文數(shù)據(jù)進行對比分析。兩者模擬的海河流域水儲量變化時間序列具有很好的一致性，證實了實測水文測量的土壤水分和地下水資料能有效地反映半干旱地區(qū)河流流域的水儲量變化。尼勝楠等[12]對CSR數(shù)據(jù)中心公布的GRACE level-2數(shù)據(jù)進行了高斯平滑濾波處理，并對水文模型數(shù)據(jù)作球諧系數(shù)展開和高斯平滑，探討了多年長江和黃河流域的水儲量在季節(jié)上和年際上的變化特征，同時通過監(jiān)測黃河流域2002年到2003年的干旱現(xiàn)象，證明了GRACE衛(wèi)星在監(jiān)測干旱事件上的能力。但由于GRACE衛(wèi)星空間分辨率較小，因此限制了其在較小流域尺度中的應用。

3.2 GRACE在湖泊、水庫流域尺度的研究應用

由于湖泊、水庫具有調(diào)蓄水量的作用，水儲量變化明顯，目前GRACE相關研究中國內(nèi)涉及的湖泊、水庫流域主要包括洞庭湖、鄱陽湖和三峽水庫等幾個大尺度流域。國內(nèi)外研究水儲量變化大多數(shù)采用的是解算好的GRACE衛(wèi)星月尺度數(shù)據(jù)。如廖夢思等[13]選用了最新解算的GRACE RL05 2003-2012年的數(shù)據(jù)與30個水文氣象站點的降水量數(shù)據(jù)相結(jié)合，通過空間插值和線性趨勢法反演出洞庭湖流域水儲量變化特征。汪漢勝等[14]結(jié)合2002年-2004年的22個月GRACE時變重力場數(shù)據(jù)通過高斯平滑法估算了三峽水庫流域的水儲量變化情況，并將反演結(jié)果與CPC(Climate Prediction Center)水文模型結(jié)果作對比分析，兩者均反映了三峽水庫流域水儲量變化特征。王富[15]在鄱陽湖的水儲量變化研究中通過15個月的平均值和高斯平滑半徑1 000 km來降低數(shù)據(jù)處理誤差，并在反演結(jié)果與VIC模型和CLM模型數(shù)據(jù)比較的基礎上進行了最小二乘法和樣條插值法擬合處理。對比發(fā)現(xiàn)最小二乘法估算結(jié)果更優(yōu)，模型表現(xiàn)出的鄱陽湖水儲量的季節(jié)和年際變化與GRACE反演結(jié)果吻合度良好。雖然目前GRACE反演的水儲量變化與水文模型模擬結(jié)果大致趨勢一致，但仍然存在一定的差異，主要原因一個是由于GRACE衛(wèi)星獲取數(shù)據(jù)及數(shù)據(jù)處理分析時產(chǎn)生的誤差，另一個原因是目前的水文模型本身的不完善導致了差異性明顯。但可以肯定的是，GRACE衛(wèi)星數(shù)據(jù)能反映出流域水儲量變化特征的總體趨勢。

3.3 GRACE在我國西部和華北等區(qū)域尺度的研究應用

由于水資源是限制我國西部和華北地區(qū)社會經(jīng)濟發(fā)展的重要因素，因此探明這些地區(qū)水儲量變得尤為重要。當前GRACE反演國內(nèi)區(qū)域水儲量的研究大多數(shù)集中在我國的西部和華北等區(qū)域，采用的數(shù)據(jù)包括GRACE衛(wèi)星數(shù)據(jù)和水文模型數(shù)據(jù)，其中GRACE數(shù)據(jù)可從CSR、GFZ、JPL三個機構(gòu)公開的不同數(shù)據(jù)產(chǎn)品中選擇，水文模型數(shù)據(jù)包括全球陸地資料同化系統(tǒng)、美國氣候預報模式和全球水資源分析和預測水文模型公布的水文模擬數(shù)據(jù)。Su等[16]運用CSR公布的GRACE衛(wèi)星八年(2002-2010年)的時變重力場資料評價了華北平原地區(qū)水儲量變化趨勢，研究發(fā)現(xiàn)在這段時間內(nèi)華北地區(qū)的水儲量時序特征表現(xiàn)為以平均每年-11 mm下降，并將反演的水儲量同CPC水文模型和GLDAS數(shù)據(jù)作了比對探討，對比結(jié)果表明GRACE反演的華北平原的水儲量變化與它們均體現(xiàn)出較好一致性，最后分析了過度的地下水開采是導致華北平原水儲量減少的重要原因。蘇子校[17]采用CSR機構(gòu)發(fā)布的GRACE level-2數(shù)據(jù)對柴達木盆地2003-2011年水儲量進行計算，研究表明該區(qū)域水儲量逐年增加，平均增速約為0.27 mm/月，總增長量約為29.46 mm，其中西南部地區(qū)增長最多。同時水儲量的季節(jié)性變化與降水量的季節(jié)性變化趨勢相同。Yao等[18]運用GRACE衛(wèi)星數(shù)據(jù)對四川盆地2001年1月-2015年2月的水儲量變化趨勢進行了評價。研究發(fā)現(xiàn)這段時期內(nèi)該區(qū)域的陸地水儲量大致表現(xiàn)出上升態(tài)勢，并且2003年和2006年的干旱事件也在水儲量變化趨勢中有所體現(xiàn)。另外，孫倩[19]和許民[20]分別探討了新疆和祁連山的陸地水儲量變化情況。

4 國內(nèi)陸地水儲量變化的因子分析

4.1 氣候因子

蘇曉莉等[21]通過GRACE數(shù)據(jù)反演了全球陸地水儲量的季節(jié)性變化趨勢。研究結(jié)果顯示只有我國和其他少數(shù)國家的流域水儲量季節(jié)性變化表現(xiàn)明顯，尤其是我國的長江流域，2005年1-7月份，其水儲量呈現(xiàn)減少態(tài)勢，并在4-6月份下降了10 cm左右，直到7-9月份水儲量才出現(xiàn)上升狀態(tài)，但在12月份和2016年的1-2月份水儲量又表現(xiàn)減少趨勢。我國氣候類型豐富，因此利用GRACE探測的地球重力場研究我國陸地水儲量變化時，降水量成為影響我國陸地水儲量變化的一個很重要的因素。我國氣候類型主要包括熱帶季風氣候、亞熱帶季風氣候、溫帶季風氣候、溫帶大陸性氣候和高原山地氣候，除西北和青藏高原地區(qū)受大陸性氣候和高原氣候控制外，我國大部分地區(qū)均受季風氣候影響，具體表現(xiàn)為夏季炎熱多雨，冬季低溫少雨。因此在探討國內(nèi)陸地水儲量變化與降水量的關系時，氣候類型因素不可忽略。

王杰等[22]將估算到的2002年8月-2010年6月期間的云南省水儲量變幅情況與相應的月降水資料對比，并采用了趨勢分析法分析了云南省水儲量的變化特征，結(jié)果表明水儲量變化與降水變化基本保持一致。李瓊等[23]通過GRACE反演的云貴川三省的水儲量變化結(jié)果與該地區(qū)的TRMM月降雨量數(shù)據(jù)進行比較，同時分析研究區(qū)陸地水儲量變化與月降雨量的時間序列變化特征。顯示該地區(qū)的水儲量變化與降雨量變化具有強相關聯(lián)系。結(jié)果也證明了降雨量的變化是導致我國陸地水儲量變化的首要因素。許民等[24]進行天山地區(qū)的水儲量變化研究時，采用了Wang等[25]的趨勢計算方法，計算出該地區(qū)的水儲量和降水量的月變化趨勢。研究表明每年的5-9月是天山地區(qū)降水主要集中的月份，在這期間該地區(qū)水儲量也比其他月份多，是一年當中水儲量最多的幾個月份。而從12月到次年的3月，降水量和水儲量則表現(xiàn)出相對匱乏狀態(tài)。結(jié)果表明該地區(qū)水儲量變化對降水量變化比較敏感。

4.2 人口社會經(jīng)濟因子

我國人口數(shù)量眾多，第一產(chǎn)業(yè)和第二產(chǎn)業(yè)發(fā)展迅速，人口社會經(jīng)濟因子在人口密度大或工農(nóng)業(yè)較為發(fā)達的地區(qū)，是影響該地區(qū)陸地水儲量變化的一個重要因素。因此在進行國內(nèi)陸地水儲量變化研究時，尤其是在農(nóng)業(yè)、工業(yè)相對發(fā)達地區(qū)，研究分析人口經(jīng)濟因子對水儲量的變化變得尤為重要。張銳[26]在對華北平原水儲量變化成因探討中，綜合考慮到華北平原是我國重要的糧食產(chǎn)地，農(nóng)業(yè)用地面積廣闊，用水量大。因此將人口密度、第一產(chǎn)業(yè)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值數(shù)據(jù)與估算到的華北平原地區(qū)的水儲量變化研究結(jié)果進行關聯(lián)度分析。研究結(jié)果顯示“第一產(chǎn)業(yè)生產(chǎn)總值”與該地區(qū)水儲量變化的灰色綜合關聯(lián)度達到了0.681，“第二產(chǎn)業(yè)生產(chǎn)總值”為0.601，均說明華北平原地區(qū)水儲量變化與人口經(jīng)濟因子聯(lián)系密切。另外，“第一產(chǎn)業(yè)生產(chǎn)總值”與華北平原地表水變化的灰色綜合關聯(lián)度是0.602。在華北平原地下水變化方面，“第一產(chǎn)業(yè)生產(chǎn)總值”的灰色綜合關聯(lián)度達到0.641，在眾多影響因子中關聯(lián)性最大，甚至高于降水的影響。說明了農(nóng)業(yè)用水對華北平原地區(qū)地下水變化影響最大。

5 GRACE在國內(nèi)水儲量研究的主要瓶頸

5.1 空間分辨率低，只能研究大尺度的陸地水儲量

GRACE在進行陸地水儲量的研究時，雖然在探測國內(nèi)幾個大的流域和區(qū)域的水儲量中得到很好的運用，但為了減少中高階位系數(shù)對研究結(jié)果的誤差影響，而只截取低階的重力場模型，與此同時也失去了中高階位帶來的重力場信息，這種方法的缺點是空間分辨率低，只能研究大尺度(350 km以上)的陸地水儲量[27]。如果它的測量面積越小，重力數(shù)據(jù)就越少，計算出水儲量信息就很困難。在廣大研究者不斷地探索下，倪葦[28]將GRACE和GOCE模型進行譜組合計算，表明了GRACE和GOCE模型可以進行互補，彌補了GRACE模型只對中低階重力場信息敏感的不足，也將GRACE在水儲量的應用領域推向更廣提供了可能。在更小尺度區(qū)域水儲量的研究上也取得進展，汪漢勝等[14]通過CPC模型與GRACE重力反演結(jié)果比較，證明了GRACE數(shù)據(jù)在監(jiān)測三峽庫區(qū)水儲量中的合理性。另外，解決GRACE衛(wèi)星空間分辨率低還有另一個方法，在區(qū)域地下水研究的基礎上，基于GRACE衛(wèi)星反演數(shù)據(jù)研究細化GRACE衛(wèi)星空間數(shù)據(jù)(降尺度)的研究方法，即在GRACE衛(wèi)星反演的區(qū)域地下水儲量動態(tài)數(shù)據(jù)基礎上，結(jié)合區(qū)域水文地質(zhì)參數(shù)的分布來研究降尺度的方法。

5.2 GRACE在水儲量研究應用中的誤差

自重力衛(wèi)星GRACE為國內(nèi)水儲量的研究提供了可靠的觀測數(shù)據(jù)，但由于數(shù)據(jù)本身和應用中的誤差使水儲量研究結(jié)果出現(xiàn)一定的偏差。水儲量誤差中主要有信號泄露、球諧系數(shù)導致的誤差[29]。

5.2.1 信號泄露誤差

信號泄露誤差是指研究區(qū)域外的地球物理信號進入了研究區(qū)域內(nèi)，或者研究區(qū)域內(nèi)的信號未考慮完全，尤其在小尺度區(qū)域的研究中誤差表現(xiàn)更為突出[30]，一般情況下，誤差大小與區(qū)域或流域的尺度大小和緯度高低具有較強的相關性[31]。均值函數(shù)[32,33]、最佳平滑核[34]和高斯濾波[35]等濾波處理方法可以通過選取最優(yōu)平滑核能一定程度上降低信號泄露誤差帶來的影響。Guo等[36]結(jié)合各向異性高斯濾波器證實了此方法可以顯著地降低信號泄露誤差。其核心是產(chǎn)品的兩個高斯函數(shù)具有明顯的橫向和縱向平滑半徑，當表示地球表面數(shù)公里區(qū)域時，縱向平滑半徑定義為一個固定的經(jīng)度間隔，這個間隔在赤道較長，在較高的緯度較短。這主要是根據(jù)GRACE數(shù)據(jù)的分辨率，以及在沒有過度平滑的緯向下產(chǎn)生均勻平滑結(jié)果。這種泄漏減少方法也可以通過去除一個模型來研究一個區(qū)域內(nèi)的信號是遠遠大于周圍的區(qū)域，或者周圍區(qū)域的信號可以被減少到非常低的幅度的區(qū)域。

5.2.2 儀器誤差

衛(wèi)星間距離變化率的系統(tǒng)噪聲誤差、超穩(wěn)定振蕩器誤差、加速度誤差以及衛(wèi)星軌道誤差即是儀器誤差。GRACE重力衛(wèi)星研究國內(nèi)水儲量情況的精確度某種程度上是由GRACE本身的誤差決定的。重力場球諧系數(shù)階數(shù)增大會引起誤差增大，因此在進行GRACE反演國內(nèi)水儲量變化的處理中需要為降低誤差的影響而截取到最優(yōu)的階數(shù)，但在截取階數(shù)的處理中又必不可少地會出現(xiàn)相應地截斷誤差。曹艷萍等[37]在監(jiān)測黑河流域水儲量變化中通過引入高斯平滑核函數(shù)W(a)，減少了來自高階項的誤差影響，得出了更為準確的反演公式。

5.2.3 混淆誤差

混淆誤差是指重力場中那些次級的短周期變化混入月平均重力場而導致的誤差。而在大部分國內(nèi)水儲量的研究上，選用的數(shù)據(jù)都是經(jīng)過模型去掉GRACE衛(wèi)星觀測數(shù)據(jù)中的大氣和海洋潮汐等的影響。例如尼勝楠等[2]在進行長江、黃河流域的水儲量變化監(jiān)測中選用的是GRACE level-2數(shù)據(jù)產(chǎn)品，該產(chǎn)品是一種已經(jīng)將非潮汐大氣、高頻海洋信號和各種潮汐的影響去除的數(shù)據(jù)。這些模型在理論上是沒有偏差的，但實際上大氣-海洋去混淆模型(AOD)都避免不了相應的誤差，這些誤差可能會導致更大誤差的潮汐模擬，影響在國內(nèi)反演水儲量的研究結(jié)果的精確度。目前混淆誤差也是制約GRACE在水儲量研究的一個重要因素。國內(nèi)在這方面的研究還比較缺乏，國外研究成果中，Seo等[38]通過GRACE合成數(shù)據(jù)和地球物理模型去除與球諧系數(shù)相關的混淆誤差。

6 結(jié) 語

自重力衛(wèi)星GRACE在國內(nèi)陸地水儲量變化研究領域越來越廣泛，并且取得了一系列重要的研究成果。利用GRACE近月平均重力場模型系列數(shù)據(jù)研究不同流域、區(qū)域的水儲量情況，有助于了解我國干旱缺水區(qū)域的水儲量變化，為國家制定發(fā)展戰(zhàn)略，促進地區(qū)社會經(jīng)濟發(fā)展和生態(tài)環(huán)境保護提供科學依據(jù)。GRACE重力衛(wèi)星在監(jiān)測陸地水儲量變化的同時還可以通過水量平衡方程進一步的研究國內(nèi)地下水儲量的變化情況，為區(qū)域地下水管理提供新的方法。目前已有不少專家學者對我國不同區(qū)域地下水儲量變化情況做了研究總結(jié)[39-41]。此外，GRACE得到的陸地水儲量變化還可用于區(qū)域蒸散發(fā)時變量[42]、計算模擬冰雪徑流和寒區(qū)凍結(jié)地的滲透量變化等領域[43]。GRACE在這些領域的發(fā)展與應用加深了我們對國內(nèi)陸地水儲量理解與認識，也為水資源的空間調(diào)度和旱澇災害預警提供新的技術支持。

GRACE的應用潛力在水文界備受關注，但由于受GRACE衛(wèi)星空間分辨率的限制，目前在國內(nèi)水儲量變化的監(jiān)測還主要集中在長江、黃河和青藏高原、華北平原等大的流域和區(qū)域尺度上，而在小流域和小區(qū)域尺度上，由于誤差較大、水文模型不完善等原因還鮮有重要的研究成果。因此GRACE的潛力還遠沒有得到最大程度的發(fā)掘。后續(xù)的重力衛(wèi)星也將被陸續(xù)研制并發(fā)射，例如具有更高分辨率的GRACE Follow-on(雙星)衛(wèi)星將于2017年升空，其觀測到的地球重力場精度比目前的GRACE高10倍。因此具有更高的時間和空間分辨率的重力衛(wèi)星數(shù)據(jù)的應用，無論是研究我國水儲量變化還是在其他相關領域都是值得期待的方向。此外，目前GRACE在國內(nèi)水儲量的研究流域和區(qū)域?qū)θ斯ひ蛩匾鸬目鐓^(qū)域調(diào)水的水儲量變化還少有研究，如南水北調(diào)工程，因此在這方面還有巨大潛力。通過對這類大的工程區(qū)域的水儲量變化研究或許可以為相關部門提供科學的數(shù)據(jù)參考和決策依據(jù)。另外，當前GRACE在國內(nèi)水儲量的研究還主要集中在傳統(tǒng)的河流和區(qū)域，與國家重大政策接軌程度欠缺。隨著“一帶一路”政策的實施和落實，探究“一帶一路”地區(qū)水儲量的變化變得尤為重要，對新“絲綢之路”沿線地區(qū)的水儲量研究，有利于科學地認識“一帶一路”地區(qū)的水儲量變化情況，深化我國當前經(jīng)濟開放水平，實現(xiàn)我國經(jīng)濟發(fā)展的再平衡。

□

[1] Hu X G, Chen J L, Zhou Y H, et al. Seasonal water storage change of the Yangtze River basin detected by GRACE[J].Science in China: Series D Earth Sciences, 2006,49(5):483-491.

[2] Baker R M L. Orbit determination from range and range-rate data[C]∥The Semi-Annual Meeting of the American Rocket Society, Los Angeles, 1960.

[3] Tapley B D, Reigber C. The GRACE mission: status and future plans[J]. EOS Trans AGU, 2000,81(48):307.

[4] 周旭華,吳 斌,彭碧波,等.全球水儲量變化的GRACE衛(wèi)星監(jiān)測[J].地球物理學報, 2006,49(6):1 644-1 650.

[5] Wahr J, Swenson S, Zlotnicki V, et al. Time-variable gravity from GRACE: first results[J]. Geophy Res Lett, 2004,31:L11501.

[6] 許厚澤,陸 洋,鐘 敏,等.衛(wèi)星重力測量及其在地球物理環(huán)境變化監(jiān)測中的應用[J].中國科學:地球科學，2012,42(6):843-853.

[7] Jekeli C. Alternative Methods to Smooth the Earth's Gravity Field[R].Columbus: Department of Geodetic Science and Surveying, 1981:1-48.

[8] Zhao Q, Wu W W, Wu Y L. Variations in China's terrestrial water storage over the past decade using GRACE data[J].Geodesy and Geodynamics, 2015,6(3):187-193.

[9] 許 民,葉柏生,趙求東.2002-2010年長江流域GRACE水儲量時空變化特征[J].地理科學進展, 2013,32(1):68-77.

[10] Luo Z C, Yao C L, Li Q, et al. Terrestrial water storage changes over the Pearl River Basin from GRACE and connections with pacific climate variability[J].Geodesy and Geodynamics, 2016,7(3):39-49.171-179.

[11] Moiwo J P, Yang Y H, Li H L, et al. Comparison of GRACE with in situ hydrological measurement data shows storage depletion in Hai River basin Northern China[J].Water S A, 2009,35(5):663-670.

[12] 尼勝楠,陳劍利,李 進,等.利用GRACE衛(wèi)星時變重力場監(jiān)測長江、黃河流域水儲量變化[J].大地測量與地球動力學, 2014,34(4):49-50.

[13] 廖夢思,章新平,黃 煌,等.利用GRACE衛(wèi)星監(jiān)測近10年洞庭湖流域水儲量變化[J].地球物理學進展, 2016,31(1):61-68.

[14] 汪漢勝,王志勇,袁旭東,等.基于GRACE時變重力場的三峽水庫補給水系水儲量變化[J].地球物理學報, 2007,50(3):730-736.

[15] 王 富.基于重力衛(wèi)星數(shù)據(jù)反演區(qū)域水儲量[D].成都：電子科技大學, 2014.

[16] Su X L, Ping J S, Ye Q X. Terrestrial Water Variations in the North China Plain Revealed by the GRACE Mission[J].Science China: Earth Sciences, 2012,42(6):917-922.

[17] 蘇子校.利用GRACE衛(wèi)星數(shù)據(jù)估算柴達木盆地水儲量變化[D].北京: 中國地質(zhì)大學, 2015.

[18] Yao C L, Luo Z C, Wang H H, GRACE-Derived Terrestrial Water Storage Changes in the Inter-Basin Region and Its Possible Influencing Factors: A Case Study of the Sichuan Basin, China[J].Remote Sense, 2016,(8):444-458.

[19] 孫 倩.基于GRACE和GLDAS的新疆水資源時空變化研究[D].烏魯木齊:新疆大學, 2015.

[20] 許 民,張世強,王 建,等.利用GRACE重力衛(wèi)星監(jiān)測祁連山水儲量時空變化[J].干旱區(qū)地理, 2014,37(3):458-467.

[21] 蘇曉莉,平勁松,黃 倩,等.重力衛(wèi)星檢測到的全球陸地水儲量變化[J].中國科學院上海天文臺年刊, 2009,30:14-21.

[22] 王 杰,黃 英,曹艷萍,等.利用GRACE重力衛(wèi)星觀測研究近7年云南省水儲量變化[J].節(jié)水灌溉, 2012,(5):1-5.

[23] 李 瓊,羅志才,鐘 波,等.利用GRACE時變重力場探測2010年中國西南干旱陸地水儲量變化[J].地球物理學報, 2013,56(6):1 843-1 849.

[24] 許 民,葉柏生,趙求東,等.利用GRACE重力衛(wèi)星監(jiān)測新疆天山山區(qū)水儲量時空變化[J].干旱區(qū)研究, 2013,30(3):404-411.

[25] Wang J, Ye B S, Liu F J, et al. Variations of NDVI over elevational zones during the past two decades and climatic controls in the Qilian Mountains, Northwestern China[J].Arctic, Antarctic, and Alpine Research, 2011,43(1):127-136.

[26] 張 銳.華北平原水儲量變化及其影響機制分析[D].北京: 首都師范大學, 2013.

[27] 姜永濤,張永志,王 帥,等.GRACE時變重力場濾波方法研究[J].測繪通報, 2014,(11):1-5.

[28] 倪 葦.GRACE和GOCE衛(wèi)星重力場模型的譜組合研究[J].測繪地理信息, 2015,40(3):66-70.

[29] 鄭秋月,陳 石.應用GRACE衛(wèi)星重力數(shù)據(jù)計算陸地水變化的進展評述.地球物理學進展, 2015,30(6):2 603-2 615.

[30] Baur O, Kuhn M, Featherstone W E. GRACE-derived ice-mass variations over Greenland by accounting for leakage effects[J].Journal of Geophysical Research: Solid Earth (1978-2012), 2009,114(B6):B06407.

[31] Seo K W, Wilson C R, Famiglietti J S, et al. Terrestrial water mass load changes from Gravity Recovery and Climate Experiment (GRACE)[J].Water Resources Research, 2006,42(5):W05417.

[32] Willis J K, Chamber D P, Nerem R S. Assessing the Globally Averaged Sea Level Budged on Seasonal to Interannual Timescales[J].Journal of Geophysical Research, 2008,113(C6):C6015.

[33] Wahr J, Molenaar M, Bryan F. Time Variability of the Earth's Gravity Field: Hydrological and Oceanic Effects and their Possible Detection Using GRACE[J].Journal of Geophysical Research, 1998,103(B12):30 205-30 230.

[34] Swenson S, Wahr J. Monitoring Changes in Continental Water Storage with GRACE[J].Space Science Reviews, 2003,108(1):345-354.

[35] Swenson S, Wahr J. Methods for Inferring Regional Surface-mass Anomalies from Gravity Recovery and Climate Experiment(GRACE)Measurements of Time-variable Gravity[J].Journal of Geophysical Research, 2002,107(B9):2 193.

[36] Guo J Y, Duan X J, Shum C K. Non-isotropic Gaussian smoothing and leakage reduction for determining mass changes over land and ocean using GRACE data[J].Geophysical Journal International,2010,181(1):290-302.

[37] 曹艷萍,南桌銅.利用GRACE重力衛(wèi)星監(jiān)測黑河流域水儲量變化, 2011,26(6):719-727.

[38] Seo K W, Wilson C R, Famiglietti J S, et al. Terrestrial water mass load changes from Gravity Recovery and Climate Experiment (GRACE)[J].Water Resources Research, 2006,42(5):W05417.

[39] Huo A D, Peng J B, Chen X H, et al. Groundwater storage and depletion trends in the Loess areas of China[J]. Environ Earth Sci, 2016,75:1 167-1 178.

[40] Xiang L W, Wang H S, Steffen H, et al. Groundwater storage changes in the Tibetan Plateau and adjacent areas revealed from GRACE satellite gravity data[J].Earth and Planetary Science Letters, 2016,449:228-239.

[41] Jiang Q N, Ferreira V G, Chen J S. Monitoring groundwater changes in the Yangtze River basin using satellite and model data[J].Arab J Geosci, 2016,9:500-508.

[42] Ramillien G, Frappart F, Guntner A, et al. Time-variations of the Regional Evapotranspiration Rate from GRACE Satellite Gravimetry[J].Water Resources Research, 2006,42:W10403.

[43] Niu G, Yang Z. Effects of Frozen Soil on Snowmelt Runoff and Soil Water Storage at a Continental Scale[J]. Journal of Hydrometeorology, 2006,7(5):937-952.