何思源，谷延超，范東明，趙鴻彬，鄭 蓉

1. 西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院，四川 成都 611756； 2. 中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司，陜西 西安 710043； 3. 中國(guó)電力工程顧問(wèn)集團(tuán)西北電力設(shè)計(jì)院有限公司，陜西 西安 710075

云南省是一個(gè)嚴(yán)重缺水的省份，近些年云南地區(qū)經(jīng)歷了持續(xù)干旱，陸地水儲(chǔ)量變化直接影響著區(qū)域經(jīng)濟(jì)和局部生態(tài)環(huán)境[1]。陸地水儲(chǔ)量的變化涉及地表水、降水、蒸發(fā)、徑流、土壤水、地下水及其他成分的相互作用。傳統(tǒng)的監(jiān)測(cè)方法(如土壤水分傳感器、探地雷達(dá)、無(wú)線傳感器等)僅對(duì)陸地水儲(chǔ)量中的地表及淺層水體進(jìn)行監(jiān)測(cè)，并不能有效監(jiān)測(cè)深層地下水變化[2]。GRACE(Gravity Recovery And Climate Experiment)衛(wèi)星于2002年3月發(fā)射，能夠以前所未有的高精度和單月甚至雙周的時(shí)間分辨率獲取地球重力場(chǎng)的中長(zhǎng)波部分[3]，用于分析全球或局部區(qū)域水文信號(hào)的季節(jié)性及長(zhǎng)期變化等[1,4,5]。因供電系統(tǒng)老化，GRACE系統(tǒng)自2011年之后出現(xiàn)間斷性數(shù)據(jù)缺失，不利于區(qū)域陸地水儲(chǔ)量的監(jiān)測(cè)，同時(shí)GRACE與其之后的GRACE Follow-on計(jì)劃將會(huì)存在數(shù)據(jù)空白[6]，因此尋找一種可替代的方法對(duì)云南區(qū)域陸地水儲(chǔ)量進(jìn)行監(jiān)測(cè)十分必要。

密集分布的GPS臺(tái)站能夠精確地記錄地表質(zhì)量變遷引起的負(fù)荷形變[7-10]，GPS負(fù)荷形變觀測(cè)量為研究地表質(zhì)量變遷(陸地水儲(chǔ)量變化)提供可選方案[11-14]。文獻(xiàn)[15]首次利用GPS垂直位移觀測(cè)量研究南加州地區(qū)陸地水儲(chǔ)量的空間分布；隨后大量的研究結(jié)果表明利用GPS垂直位移形變能夠較好地反演陸地水以及地下水的季節(jié)性信號(hào)及長(zhǎng)期變化[16-17]。本文首先分析利用云南省及其周邊現(xiàn)有的47個(gè)中國(guó)大陸環(huán)境構(gòu)造監(jiān)測(cè)網(wǎng)絡(luò)(簡(jiǎn)稱(chēng)陸態(tài)網(wǎng))連續(xù)運(yùn)行觀測(cè)站反演云南地區(qū)陸地水儲(chǔ)量的可行性，然后利用實(shí)際觀測(cè)的GPS垂直位移形變數(shù)據(jù)反演2010—2014年云南省陸地水儲(chǔ)量變化，最后將其與GRACE、GLDAS以及降雨數(shù)據(jù)進(jìn)行綜合對(duì)比分析，以驗(yàn)證GPS觀測(cè)數(shù)據(jù)反演云南地區(qū)陸地水儲(chǔ)量變化的可行性。

1 基本原理與方法

1.1 地表垂直負(fù)荷形變

地表水、雪、冰以及大氣作用于固體地球時(shí)，地球表面會(huì)產(chǎn)生負(fù)荷形變，利用格林函數(shù)可建立負(fù)荷質(zhì)量與形變之間的關(guān)系[18]。負(fù)荷形變包含水平方向和垂直方向，垂直位移形變對(duì)負(fù)荷源較為敏感，其負(fù)荷形變量約為水平負(fù)荷形變的2～3倍[19]。負(fù)荷引起的垂直位移的格林函數(shù)為[20]

(1)

式中，Pn為勒讓德多項(xiàng)式；θ為距離圓盤(pán)中心的角半徑；G為牛頓萬(wàn)有引力常數(shù)；a為地球半徑；hn負(fù)荷勒夫數(shù)；g為重力加速度；Γn函數(shù)的推導(dǎo)公式如下

(2)

(3)

圖1(a)為相同質(zhì)量不同半徑和厚度的圓盤(pán)置于地表，產(chǎn)生的負(fù)荷形變量與距離的關(guān)系圖，可知：近場(chǎng)的負(fù)荷響應(yīng)非常明顯，當(dāng)距圓盤(pán)中心的距離為圓盤(pán)半徑的3倍時(shí)，負(fù)荷形變僅為中心負(fù)荷形變量的1/10，當(dāng)距圓盤(pán)中心的距離為圓盤(pán)半徑的10倍時(shí)，這種響應(yīng)可忽略不計(jì)。因此臺(tái)站垂直位移形變主要受周邊有限范圍內(nèi)的負(fù)荷源影響，當(dāng)采用GPS反演局部區(qū)域陸地水儲(chǔ)量時(shí)，邊界附近的臺(tái)站同時(shí)受到區(qū)域內(nèi)和區(qū)域外的負(fù)荷源的影響，若忽略邊界外的質(zhì)量變化，則會(huì)造成邊界外的負(fù)荷源被約束至邊界區(qū)域，使得邊界處反演結(jié)果不正確。因此反演過(guò)程中實(shí)際反演區(qū)域比研究區(qū)域范圍大，應(yīng)顧及負(fù)荷效應(yīng)響應(yīng)的敏感范圍進(jìn)行反演區(qū)域的擴(kuò)展[16]。本文討論了不同擴(kuò)展范圍情況下同一邊界點(diǎn)的反演等效水高值，圖1(b)為當(dāng)空間格網(wǎng)為0.25°時(shí)反演的邊界點(diǎn)等效水高與擴(kuò)展半徑間的關(guān)系，結(jié)果表明：隨著擴(kuò)展范圍的增大，計(jì)算得到的同一邊界點(diǎn)的陸地水(terrestrial water storage,TWS)等效水高值趨于穩(wěn)定，當(dāng)擴(kuò)展范圍2.5°時(shí)反演結(jié)果已基本達(dá)到穩(wěn)定，結(jié)合文獻(xiàn)[16]中對(duì)邊界的擴(kuò)展討論，根據(jù)本研究區(qū)域多個(gè)邊界點(diǎn)的反饋結(jié)果選取2°為擴(kuò)展邊界范圍。

圖1 垂直位移負(fù)荷形變響應(yīng)示意圖及反演邊界擴(kuò)展結(jié)果Fig.1 Diagrammatic diagram showing the vertical response for mass disk and the inversion results influenced by boundary extension

1.2 反演模型

本文沿用文獻(xiàn)[15]的反演模型

(4)

式中，σ為垂直位移觀測(cè)值向量標(biāo)準(zhǔn)差；A為格林函數(shù)系數(shù)矩陣，由式(1)求得；x為格網(wǎng)點(diǎn)對(duì)應(yīng)的等效水高；b為測(cè)站GPS垂直位移觀測(cè)量；β為平滑因子；L為拉普拉斯算子。

由于區(qū)域內(nèi)GPS觀測(cè)臺(tái)站數(shù)量有限導(dǎo)致方程個(gè)數(shù)小于未知數(shù)個(gè)數(shù)，導(dǎo)致法方程系數(shù)矩陣秩虧，因此利用GPS觀測(cè)數(shù)據(jù)反演陸地水儲(chǔ)量是一個(gè)不適定方程問(wèn)題。Hoerl提出的嶺估計(jì)法是一種經(jīng)典的正則化方法，方法關(guān)鍵在于嶺參數(shù)的求取，即式(4)反演模型中β的求解。本文采用L曲線法選取最優(yōu)的正則化參數(shù)，具體方法與實(shí)現(xiàn)過(guò)程請(qǐng)參考文獻(xiàn)[21]，如圖2所示，確定最優(yōu)的嶺參數(shù)β為0.01。

2 GPS觀測(cè)數(shù)據(jù)反演陸地水儲(chǔ)量可行性分析及反演結(jié)果

2.1 GPS觀測(cè)數(shù)據(jù)

本文采用的GPS時(shí)間序列來(lái)源于中國(guó)大陸環(huán)境監(jiān)測(cè)網(wǎng)GAMIT/GLOBK解算結(jié)果(ftp:∥ftp.cgps.ac.cn)，數(shù)據(jù)處理過(guò)程中已扣除了固體潮、海潮和極潮等影響，垂直位移負(fù)荷形變主要反映非潮汐大氣、海洋和水文的影響[7,10]。為了研究陸地水儲(chǔ)量的變化，利用大氣模型Merra和非潮汐海洋OMCT模型[22]扣除大氣負(fù)荷和非潮汐海洋負(fù)荷的影響。GPS臺(tái)站單日解的中誤差較大[8,13]，為了消除隨機(jī)噪聲的影響，通過(guò)月內(nèi)平均得到GPS臺(tái)站月垂直位移時(shí)間序列，進(jìn)而反演得到單月的云南地區(qū)陸地水儲(chǔ)量。本文選取97°E～107°E，20°N～30°N范圍內(nèi)的47個(gè)GPS臺(tái)站，其點(diǎn)位分布如圖3(a)所示，臺(tái)站的觀測(cè)數(shù)據(jù)時(shí)段為2010—2014年。圖3(a)為離散點(diǎn)通過(guò)克里金內(nèi)插得到的云南地區(qū)的垂直位移周年振幅，圖3(b)為部分臺(tái)站(YNML、YNLC和YNGM)垂直位移時(shí)間序列?？煽闯觯篩NML、YNLC臺(tái)站均有明顯的季節(jié)性信號(hào)，然而YNGM臺(tái)站時(shí)間序列的季節(jié)性信號(hào)不明顯，數(shù)據(jù)異常[23]，數(shù)據(jù)分析顯示除YNGM外，其他臺(tái)站數(shù)據(jù)正常，另外本次所用數(shù)據(jù)中部分臺(tái)站個(gè)別月份數(shù)據(jù)缺失，因此采用主成分分析補(bǔ)齊缺失數(shù)據(jù)并消除數(shù)據(jù)異常；云南省西南部地區(qū)周年振幅偏大，約為11 mm，東北部地區(qū)偏小，周年振幅約為5 mm。

2.2 利用當(dāng)前GPS臺(tái)站反演陸地水可行性分析

GPS觀測(cè)精度受儀器設(shè)備和相關(guān)模型缺陷的影響，目前GPS垂直方向的精度較低，利用當(dāng)前的GPS垂直位移觀測(cè)量能否給出穩(wěn)健的陸地水儲(chǔ)量結(jié)果仍需深入的研究和分析。為了探討當(dāng)前GPS臺(tái)站反演陸地水的可行性，本文選取2010—2014年期間全球陸地?cái)?shù)據(jù)同化模型GLDAS Noah數(shù)據(jù)(http:∥daac.gsfc.nasa.gov)，計(jì)算得到目標(biāo)區(qū)域內(nèi)水儲(chǔ)量周年振幅，以其周年振幅為云南區(qū)域內(nèi)陸地水分布真值，根據(jù)式(2)計(jì)算各臺(tái)站的模擬垂直位移Sup，模擬隨機(jī)分布的觀測(cè)噪聲Δ，通過(guò)包含誤差的垂直位移Sup+Δ反演陸地水儲(chǔ)量。根據(jù)數(shù)據(jù)預(yù)處理得到GPS數(shù)據(jù)周年振幅平均誤差為0.8 mm，因此對(duì)模擬負(fù)荷形變量施加Δnorm(0，0.8)隨機(jī)誤差。圖4(a)為水文模型GLDAS陸地水周年振幅，(b)為反演得到的陸地水儲(chǔ)量周年振幅。反演結(jié)果與真實(shí)的陸地水分布的相關(guān)系數(shù)為0.9，具有很強(qiáng)的一致性：云南省境內(nèi)反演的等效水高在數(shù)值以及空間分布上與實(shí)際結(jié)果相符；但對(duì)于云南省周邊觀測(cè)點(diǎn)較少的地區(qū)(A、B、C、D區(qū)域)，反演結(jié)果與實(shí)際分布差異較大。研究表明反演結(jié)果受到點(diǎn)位密集程度的影響[16]，點(diǎn)位越密集反演結(jié)果越穩(wěn)健，而對(duì)于點(diǎn)位稀疏的區(qū)域，反演結(jié)果相對(duì)較差。本文選取的GPS臺(tái)站主要位于云南省內(nèi)及北部周邊區(qū)域，云南省內(nèi)以及圖中E區(qū)反演結(jié)果與真值具有很好的一致性，而對(duì)于云南省區(qū)域外如圖中A、B、C、D區(qū)域反演結(jié)果較差。為了測(cè)試隨機(jī)誤差對(duì)反演結(jié)果的影響，本文在等條件下進(jìn)行1000次上述試驗(yàn)，試驗(yàn)統(tǒng)計(jì)分析結(jié)果表明：反演得到的陸地水儲(chǔ)量與真值的平均相關(guān)系數(shù)為0.8，空間分布上具有很好的一致性。模擬試驗(yàn)表明當(dāng)前云南地區(qū)及其周邊的47個(gè)GPS臺(tái)站垂直位移觀測(cè)數(shù)據(jù)能夠有效地反映云南地區(qū)陸地水儲(chǔ)量變化信息。

圖2 利用L-曲線圖進(jìn)行正則化參數(shù)選取Fig.2 Optimized regularization parameter based on L-curve between solution roughness and observation misfits

圖3 云南地區(qū)垂直位移周年信號(hào)及觀測(cè)臺(tái)站垂直分量的時(shí)間序列Fig.3 Annual signal result and GPS time series of vertical component in Yunnan province and adjacent areas

2.3 陸地水儲(chǔ)量變化計(jì)算結(jié)果

采用扣除非潮汐海洋和大氣負(fù)荷后的GPS月垂直位移時(shí)間序列，選取上述參數(shù)分析的結(jié)果(取擴(kuò)展范圍2°，平滑因子β=0.01，預(yù)設(shè)圓盤(pán)半徑16 km)，根據(jù)式(5)計(jì)算云南省2010—2014年0.25×0.25°陸地水變化，云南省陸地水儲(chǔ)量變化具有明顯的地域性和季節(jié)性，限于篇幅所限，文中僅給出2011年云南省陸地水儲(chǔ)量變化時(shí)間序列(如圖5所示)。

3 綜合對(duì)比分析

為了分析GPS反演的陸地水儲(chǔ)量變化的可靠性，本文利用GRACE數(shù)據(jù)、GLDAS數(shù)據(jù)和TRMM降雨數(shù)據(jù)進(jìn)行驗(yàn)證。文中采用CSR(Center For Space Research)發(fā)布的2010—2014年的GRACE RL05 Mascon數(shù)據(jù)[24]，降雨數(shù)據(jù)為全球0.25°×0.25°分辨率的TRMM(Tropical Rainfall Measurement Mission，https:∥pmm.nasa.gov)格網(wǎng)數(shù)據(jù)。圖6為云南省(a)GPS、(b)GRACE和(c)GLDAS陸地水儲(chǔ)量周年振幅以及(d)TRMM降雨量周年振幅；圖7為GPS、GRACE、GLDAS和TRMM數(shù)據(jù)的時(shí)間序列變化圖。

圖4 利用當(dāng)前GPS臺(tái)站反演陸地水儲(chǔ)量變化的模擬試驗(yàn)Fig.4 Simulation experiments using the GPS data to infer the variation of water storage

圖5 云南省2011年GPS陸地水儲(chǔ)量變化Fig.5 The terrestrial water storage of Yunnan province inferred from GPS vertical displacements in 2011

圖6 云南省陸地水儲(chǔ)量或降雨數(shù)據(jù)周年振幅Fig.6 Annual amplitude of TWS or precipitation data in Yunnan province

圖7 云南省GPS、GRACE、GLDAS的陸地水儲(chǔ)量及降雨數(shù)據(jù)時(shí)間序列Fig.7 TWS time series from GPS，GRACE and GLDAS and TRMM data

(1) GPS陸地水儲(chǔ)量周年振幅呈明顯東西分化趨勢(shì)，西部山地區(qū)域遠(yuǎn)大于東部盆地地區(qū)，陸地水周年振幅最大值約為350 mm；最小值為約90 mm，與降雨數(shù)據(jù)空間分布基本吻合。云南省大部分地區(qū)水儲(chǔ)量季節(jié)性變化在10月(夏季末)達(dá)到最大值，在4月(冬季末)達(dá)到最小值。GPS陸地水儲(chǔ)量周年振幅在空間分布上與GRACE和GLDAS陸地水儲(chǔ)量具有較高一致性，計(jì)算結(jié)果顯示GPS陸地水周年振幅與GRACE、GLDAS的相關(guān)系數(shù)分別為0.79和0.71，但從圖6和圖7中可以發(fā)現(xiàn)：GPS計(jì)算得到的云南省陸地水周年振幅約為90～350 mm，而GRACE和GLDAS的結(jié)果約為50～250 mm；GPS的周年振幅在數(shù)值上明顯大于GRACE和GLDAS，云南省GRACE陸地水周年振幅僅為GPS反演結(jié)果的79%左右。

(2) GPS陸地水儲(chǔ)量長(zhǎng)期變化趨勢(shì)(如圖7所示)表明2010—2014年期間云南省陸地水儲(chǔ)量以20 mm/a的速度緩慢增長(zhǎng)，該結(jié)果與GRACE時(shí)間序列的長(zhǎng)期趨勢(shì)一致。通過(guò)GRACE與GLDAS數(shù)據(jù)結(jié)果得到的云南省地下水變化的時(shí)間序列，結(jié)果表明云南省2010—2012年期間地下水處于快速增長(zhǎng)的趨勢(shì)，速度大約為50 mm/a，2013年后無(wú)顯著變化。

(3) GPS垂直觀測(cè)量受多種誤差源影響，主要包含衛(wèi)星及接收機(jī)鐘差、對(duì)流層及電離層延遲、多路徑影響等，均影響著陸地水儲(chǔ)量的反演精度。圖7中2012年的1月GPS陸地水儲(chǔ)量顯著低于相鄰月份和2014年的3月GPS陸地水儲(chǔ)量顯著高于相鄰月份。通過(guò)對(duì)比GRACE、GLDAS以及TRMM數(shù)據(jù)發(fā)現(xiàn)該時(shí)間并無(wú)異常極端氣候，說(shuō)明2012年的1月和2014年的3月GPS陸地水儲(chǔ)量異常值是由GPS垂直觀測(cè)量誤差引起的，并非真實(shí)的陸地水儲(chǔ)量變化信號(hào)。

上述結(jié)果表明：利用GPS反演的陸地水儲(chǔ)量變化與GRACE、GLDAS以及TRMM具有較好的時(shí)空一致性，GPS和GRACE陸地水儲(chǔ)量長(zhǎng)期變化趨勢(shì)一致。但除個(gè)別月份GPS陸地水儲(chǔ)量變化因隨機(jī)噪聲的影響與其他數(shù)據(jù)存在差異外，GPS反演結(jié)果周年振幅明顯大于其他技術(shù)手段。因GPS和GRACE均能夠反映陸地水儲(chǔ)量所有成分的變化，針對(duì)這種結(jié)果上的差異，圖8為選取的YNML和YNGM臺(tái)站處的對(duì)應(yīng)的GPS和GRACE單月等效水高散點(diǎn)分析圖，橫軸和縱軸分別表示GPS和GRACE得到的等效水高，若兩者完全一致，則離散點(diǎn)應(yīng)位于“k=1”標(biāo)準(zhǔn)直線附近。

由圖8可看出：利用直線擬合離散點(diǎn)分布時(shí)，離散點(diǎn)比較集中在擬合直線附近，且實(shí)際擬合直線的斜率明顯小于標(biāo)準(zhǔn)直線，說(shuō)明GPS陸地水儲(chǔ)量周年振幅明顯大于GRACE陸地水儲(chǔ)量的周年振幅。本文通過(guò)計(jì)算得到所用47個(gè)臺(tái)站的均方根誤差為0.8 mm，與前期陸態(tài)網(wǎng)測(cè)站垂直方向觀測(cè)精度一致(平均單日重復(fù)度為3.9 mm，臺(tái)站數(shù)據(jù)的月重復(fù)度在0.5～1.2 mm，平均為0.7 mm[13])，故利用該數(shù)值進(jìn)行可行性分析，結(jié)果表明在該誤差水平的影響下利用GPS得到的陸地水儲(chǔ)量是穩(wěn)健可靠的。因此GPS陸地水儲(chǔ)量周年振幅偏大并非由于白噪聲產(chǎn)生，而是由其他周年地球物理信號(hào)或者虛假信號(hào)引起的。

圖8 由GPS和GRACE得到的YNML和YNGM兩站點(diǎn)的單月水高等效水高對(duì)比圖Fig.8 Comparison of monthly water thickness inferred from GPS and GRACE at sites YNML and YNGM

已扣除固體潮、海潮、極潮和非潮汐大氣海洋效應(yīng)的GPS觀測(cè)值僅包含水文信號(hào)或其他虛假信號(hào)，就目前對(duì)GPS周年信號(hào)的分析而言，部分原因(如：交年項(xiàng)誤差[6,13]、海潮建模不準(zhǔn)確引起的周年信號(hào)[25])可產(chǎn)生虛假的周年信號(hào)。云南地區(qū)并非沿海地區(qū)，其海潮模型不準(zhǔn)確而引起的周年信號(hào)較小[26]。文獻(xiàn)[16]通過(guò)對(duì)GPS臺(tái)站長(zhǎng)時(shí)間序列的分析得到交年項(xiàng)振幅占周年項(xiàng)振幅的平均比例為21.24%±13.56%，并且指出交年項(xiàng)誤差作為GPS時(shí)間序列的潛在系統(tǒng)誤差，會(huì)產(chǎn)生明顯虛假的物理信號(hào)[16]。本文采用GRACE數(shù)據(jù)為Mascon求解結(jié)果，泄漏誤差較小，且并未采用任何空間濾波，故能夠較真實(shí)地反映陸地水儲(chǔ)量信號(hào)。因本文所用數(shù)據(jù)觀測(cè)時(shí)間較短，故難以對(duì)交年項(xiàng)誤差進(jìn)行準(zhǔn)確估計(jì)，但根據(jù)文獻(xiàn)[16]的研究結(jié)果可知GPS周年振幅大于GRACE結(jié)果的原因可能是因交年項(xiàng)誤差引起的虛假信號(hào)，但該誤差對(duì)陸地水儲(chǔ)量長(zhǎng)期變化趨勢(shì)影響不大。

4 結(jié) 論

本文在負(fù)荷形變的基礎(chǔ)上，通過(guò)模擬試驗(yàn)驗(yàn)證了當(dāng)前云南地區(qū)的GPS垂直位移觀測(cè)數(shù)據(jù)可作為獨(dú)立觀測(cè)量反演陸地水儲(chǔ)量變化，為陸地水儲(chǔ)量監(jiān)測(cè)提供另外一種可行的方案。GPS反演結(jié)果的準(zhǔn)確性與觀測(cè)數(shù)據(jù)精度和點(diǎn)的密集程度相關(guān)，點(diǎn)位越密集反演結(jié)果的空間分辨率越高，觀測(cè)數(shù)據(jù)的精度越高反演結(jié)果越準(zhǔn)確。利用陸態(tài)網(wǎng)2010—2014年的觀測(cè)數(shù)據(jù)反演云南省陸地水儲(chǔ)量變化，綜合對(duì)比分析表明：GPS陸地水儲(chǔ)量與GRACE、GLDAS結(jié)果具有較好的一致性；云南省水儲(chǔ)量周年振幅在空間分布上西部高山高原地區(qū)大于東部盆地地區(qū)；云南省陸地水總體變化趨勢(shì)以20 mm/a速度增長(zhǎng)，地下水2010—2012年以50 mm/a速度增長(zhǎng)。GPS反演陸地水的方法可為監(jiān)測(cè)云南地區(qū)陸地水儲(chǔ)量提供參考，反演結(jié)果可用于當(dāng)前GRACE缺失月份數(shù)據(jù)的補(bǔ)充，同時(shí)也可用于GRACE和GRACE Follow-on期間的云南地區(qū)陸地水儲(chǔ)量的監(jiān)測(cè)和研究。

[1] 李瓊, 羅志才, 鐘波, 等. 利用GRACE時(shí)變重力場(chǎng)探測(cè)2010年中國(guó)西南干旱陸地水儲(chǔ)量變化[J]. 地球物理學(xué)報(bào), 2013, 56(6): 1843-1849.

LI Qiong, LUO Zhicai, 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, 56(6): 1843-1849.

[2] CHEW C C, SMALL E E. Terrestrial Water Storage Response to the 2012 Drought Estimated from GPS Vertical Position Anomalies[J]. Geophysical Research Letters, 2014, 41(17): 6145-6151.

[3] TAPLEY B D, BETTADPUR S, WATKINS M, et al. The Gravity Recovery and Climate Experiment: Mission Overview and Early Results[J]. Geophysical Research Letters, 2004, 31(9): L09607.

[4] 鐘敏, 段建賓, 許厚澤, 等. 利用衛(wèi)星重力觀測(cè)研究近5年中國(guó)陸地水量中長(zhǎng)空間尺度的變化趨勢(shì)[J]. 科學(xué)通報(bào), 2009, 54(9): 1290-1294.

ZHONG Min, DUAN Jianbin, XU Houze, et al. Trend of China Land Water Storage Redistribution at Medi-and Large-spatial Scales in Recent Five Years by Satellite Gravity Observations[J]. Chinese Science Bulletin, 2009, 54(9): 1290-1294.

[5] 盧飛, 游為, 范東明, 等. 由GRACE RL05數(shù)據(jù)反演近10年中國(guó)大陸水儲(chǔ)量及海水質(zhì)量變化[J]. 測(cè)繪學(xué)報(bào), 2015, 44(2): 160-167. DOI: 10.11947/j.AGCS.2015.20130753.

LU Fei, YOU Wei, FAN Dongming, et al. Chinese Continental Water Storage and Ocean Water Mass Variations Analysis in Recent Ten Years Based on GEACE RL05 Data[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(2): 160-167. DOI: 10.11947/j.AGCS.2015.20130753.

[6] FAMIGLIETTI J S, RODELL M. Water in the Balance[J]. Science, 2013, 340(6138): 1300-1301.

[7] DAVIS J L, ELSEGUI P, MITROVICA J X, et al. Climate-driven Deformation of the Solid Earth from GRACE and GPS[J]. Geophysical Research Letters, 2004, 31(24): L24605.

[8] NAHMANI S, BOCK O, BOUIN M N, et al. Hydrological Deformation Induced by the West African Monsoon: Comparison of GPS, GRACE and Loading Models[J]. Journal of Geophysical Research: Solid Earth, 2012, 117(B5): B05409.

[9] 姜衛(wèi)平, 夏傳義, 李昭, 等. 環(huán)境負(fù)載對(duì)區(qū)域GPS基準(zhǔn)站時(shí)間序列的影響分析[J]. 測(cè)繪學(xué)報(bào), 2014, 43(12): 1217-1223. DOI: 10.13485/j.cnki.11-2089.2014.0149.

JIANG Weiping, XIA Chuanyi, LI Zhao, et al. Analysis of Environmental Loading Effects on Regional GPS Coordinate Time Series[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(12): 1217-1223. DOI: 10.13485/j.cnki.11-2089.2014.0149.

[10] 盛傳貞, 甘衛(wèi)軍, 梁詩(shī)明, 等. 滇西地區(qū)GPS時(shí)間序列中陸地水載荷形變干擾的GRACE分辨與剔除[J]. 地球物理學(xué)報(bào), 2014, 57(1): 42-52.

SHENG Chuanzhen, GAN Weijun, LIANG Shiming, et al. Identification and Elimination of Non-tectonic Crustal Deformation Caused by Land Water from GPS Time Series in the Western Yunnan Province based on GRACE Observations[J]. Chinese Journal of Geophysics, 2014, 57(1): 42-52.

[11] BORSA A A, AGNEW D C, CAYAN D R. Ongoing Drought-induced Uplift in the Western United States[J]. Science, 2014, 345(6204): 1587-1590.

[12] BIRHANU Y, BENDICK R. Monsoonal Loading in Ethiopia and Eritrea from Vertical GPS Displacement Time Series[J]. Journal of Geophysical Research: Solid Earth, 2015, 120(10): 7231-7238.

[13] GU Yanchao, YUAN Linguo, FAN Dongming, et al. Seasonal Crustal Vertical Deformation Induced by Environmental Mass Loading in Mainland China Derived from GPS, GRACE and Surface Loading Models[J]. Advances in Space Research, 2017, 59(1): 88-102.

[14] 王林松, 陳超, 鄒蓉, 等. 利用GPS與GRACE監(jiān)測(cè)陸地水負(fù)荷導(dǎo)致的季節(jié)性水平形變: 以喜馬拉雅山地區(qū)為例[J]. 地球物理學(xué)報(bào), 2014, 57(6): 1792-1804.

WANG Linsong, CHEN Chao, ZOU Rong, et al. Using GPS and GRACE to Detect Seasonal Horizontal Deformation Caused by Loading of Terrestrial Water: A Case Study in the Himalayas[J]. Chinese Journal of Geophysics, 2014, 57(6): 1972-1804.

[15] ARGUS D F, FU Yuning, LANDERER F W. Seasonal Variation in Total Water Storage in California Inferred from GPS Observations of Vertical Land Motion[J]. Geophysical Research Letters, 2014, 41(6): 1971-1980.

[16] FU Yuning, ARGUS D F, LANDERER F W. GPS as An Independent Measurement to Estimate Terrestrial Water Storage Variations in Washington and Oregon[J]. Journal of Geophysical Research: Solid Earth, 2015, 120(1): 552-566.

[17] JIN Shuanggen, ZHANG Tengyu. Terrestrial Water Storage Anomalies Associated with Drought in Southwestern USA from GPS Observations[J]. Surveys in Geophysics, 2016, 37(6): 1139-1156.

[18] FARRELL W E. Deformation of the Earth by Surface Loads[J]. Reviews of Geophysics, 1972, 10(3): 761-797.

[19] WAHR J, KHAN S A, VAN DAM T, et al. The Use of GPS Horizontals for Loading Studies, with Applications to Northern California and Southeast Greenland[J]. Journal of Geophysical Research: Solid Earth, 2013, 118(4): 1795-1806.

[20] WANG Hansheng, XIANG Longwei, JIA Lulu, et al. Load Love Numbers and Green’s Functions for Elastic Earth Models PREM, iasp91, ak135, and Modified Models with Refined Crustal Structure from Crust 2.0[J]. Computers & Geosciences, 2012, 49: 190-199.

[21] 王振杰, 歐吉坤. 用L-曲線法確定嶺估計(jì)中的嶺參數(shù)[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2004, 29(3): 235-238.

WANG Zhenjie, OU Jikun. Determining the Ridge Parameter in a Ridge Estimation Using L-curve Method[J]. Geomatics and Information Science of Wuhan University, 2004, 29(3): 235-238.

[22] PETROV L. The International Mass Loading Service[EB/OL]. [2015-02-28]. http:∥arxiv.org/abs/1503.00191.

[23] ZHAN Wei, LI Fei, HAO Weifeng, et al. Regional Characteristics and Influencing Factors of Seasonal Vertical Crustal Motions in Yunnan, China[J]. Geophysical Journal International, 2017, 210(3): 1295-1304.

[24] SAVE H, BETTADPUR S, TAPLEY B D. High-Resolution CSR GRACE RL05 Mascons[J]. Journal of Geophysical Research: Solid Earth, 2016, 121(10): 7547-7569.

[25] PENNA N T, KING M A, STEWART M P. GPS Height Time Series: Short-Period Origins of Spurious Long-Period Signals[J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B2): B02402.

[26] 廖海華, 鐘敏, 周旭華. 利用GRACE衛(wèi)星重力資料解算氣候驅(qū)動(dòng)的地表周年垂直形變[J]. 地球物理學(xué)報(bào), 2010, 53(5): 1091-1098.

LIAO Haihua, ZHONG Min, ZHOU Xuhua. Climate-Driven Annual Vertical Deformation of the Solid Earth Calculated from GRACE[J]. Chinese Journal of Geophysics, 2010, 53(5): 1091-1098.