利用GRACE模型反演長江流域陸地水儲(chǔ)量的變化

王 微,李偉偉,王奉偉

(1.同濟(jì)大學(xué) 測繪與地理信息學(xué)院，上海 200092； 2.山東科技大學(xué) 測繪科學(xué)與工程學(xué)院，山東 青島 266590)

1 研究背景

長江是我國最長的河流，其流域面積約1.8×106km2，橫跨19個(gè)省、市、自治區(qū)。長江流域呈多級階梯地形，上中游以湖北省宜昌為界，中下游以江西省湖口為界。受季風(fēng)影響，流域內(nèi)降雨量年內(nèi)分布不均，而陸地水儲(chǔ)量變化能綜合反映降雨、蒸散、徑流和地下水的變化情況。因此，研究長江流域的陸地水儲(chǔ)量變化，對了解其水資源變化以及合理安排水資源的調(diào)配等均具有重要意義。

GRACE重力衛(wèi)星計(jì)劃由NASA和DLR聯(lián)合實(shí)施，能夠監(jiān)測空間中長尺度的時(shí)變信息，同時(shí)還可以揭示厘米級的等效水高變化[1-2]，目前已被廣泛應(yīng)用于全球與區(qū)域的陸地水儲(chǔ)量變化研究中。陳劍利和馮偉等[3-5]，利用月時(shí)變重力場模型反演了亞馬遜流域陸地水儲(chǔ)量的變化，并成功監(jiān)測到了該流域2005年和2010年發(fā)生的特大干旱以及2009年的洪水災(zāi)害。許多學(xué)者利用GRACE監(jiān)測技術(shù)，在長江流域展開了研究。比如：許民等定量研究了水儲(chǔ)量變化速率及時(shí)空分布特性[6]，結(jié)果發(fā)現(xiàn)2002～2010年全流域水儲(chǔ)量的月增長速率為0.43 mm，并分別給出了流域上中下游的水儲(chǔ)量變化速率。超能芳等基于CSR RL05模型數(shù)據(jù)，驗(yàn)證了不同濾波方法在長江流域水儲(chǔ)量變化監(jiān)測中的適用性[7]，結(jié)果發(fā)現(xiàn)各向異性濾波更適用于區(qū)域水儲(chǔ)量變化研究。李杰等研究了2003年1月至2012年12月長江中下游平原陸地水儲(chǔ)量變化的趨勢[8]，并對其干旱的成因進(jìn)行了分析。孫麗等利用CSR RL05模型求得的2003～2013年長江流域陸地水儲(chǔ)量變化的增速約為0.6 cm/a[9]。田小娟等基于CSR RL05模型，研究了2002～2016年溫度、降水、人類活動(dòng)以及ENSO等對長江流域水儲(chǔ)量變化的影響[10]。不同的研究時(shí)間段，陸地水儲(chǔ)量的增速不同；不同的時(shí)變重力場模型，其反演結(jié)果及其精度也不盡相同[11]。

目前，具體到流域內(nèi)各省水儲(chǔ)量變化的研究較少，而且大多是基于CSR RL05模型進(jìn)行研究。相比于CSR RL05模型，CSR RL06采用了新的背景模型并強(qiáng)化了對條帶誤差的抑制，精度有了一定的提高[12]。而同濟(jì)大學(xué)采用優(yōu)化的短弧法解算的Tongji-Grace2018時(shí)變重力場模型[11]，于2019年發(fā)布在ICGEM官方網(wǎng)站(http:∥icgem.gfz-potsdam.de/series)，其噪聲比CSR RL05和RL06模型更小，反演結(jié)果表明：其信噪比更大。因此，本文采用Tongji-Grace2018與CSR RL06模型進(jìn)行實(shí)驗(yàn)分析，定量估計(jì)長江流域2003年1月至2016年8月的陸地水儲(chǔ)量變化，對將近14 a的陸地水儲(chǔ)量變化趨勢及周期進(jìn)行了分析，主要研究其月平均和季節(jié)平均的陸地水儲(chǔ)量時(shí)空分布特性，并定量分析流域內(nèi)各省的陸地水儲(chǔ)量變化情況。

2 數(shù)據(jù)與反演方法

本文選用96階Tongji-Grace2018和CSR RL06模型，時(shí)間跨度為2003年1月至2016年8月。由于衛(wèi)星傳感器等方面的問題，導(dǎo)致某些月份數(shù)據(jù)缺失，研究時(shí)間段內(nèi)，Tongji-Grace2018模型有148個(gè)月的數(shù)據(jù)，CSR RL06有149個(gè)月的數(shù)據(jù)。此外，重力衛(wèi)星軌道的幾何構(gòu)型對重力場低階項(xiàng)不敏感，其C20項(xiàng)精度較差，因此，替換為SLR解算的C20項(xiàng)[13]。直接利用GRACE模型反演的質(zhì)量變化存在著明顯的南北條帶誤差和高頻誤差，而且模型誤差隨球諧系數(shù)階數(shù)的增加而增加，所以將球諧系數(shù)截?cái)嗟?0階[14]，采用扇形濾波抑制高階項(xiàng)誤差影響，結(jié)合去相關(guān)濾波來減少條帶誤差[15-16]。本文取研究時(shí)間段內(nèi)的平均值作為背景場，將其作為衡量陸地水儲(chǔ)量變化的標(biāo)準(zhǔn)，當(dāng)某月陸地水儲(chǔ)量變化為負(fù)值時(shí)，說明該月份的陸地水儲(chǔ)量低于研究時(shí)間段內(nèi)的陸地水儲(chǔ)量的平均值，以“虧損”描述，反之，以“盈余”表示。

GRACE時(shí)變重力場反演的陸地水儲(chǔ)量變化通常用以cm為單位的等效水高表示，數(shù)值上為質(zhì)量變化面密度與水密度的比值，具體公式如下[17]：

(1)

Wl,m=Wl×Wm

(2)

Wl根據(jù)式(3)遞推得到

(3)

將l替換為m可得Wm。

由于GRACE時(shí)變重力場模型系數(shù)同一次的奇(偶)數(shù)階之間存在著相關(guān)性，采用多項(xiàng)式擬合高階系數(shù)并從原系數(shù)中扣除擬合項(xiàng)，以消除位系數(shù)之間的相關(guān)性[17]。本文利用P4M6進(jìn)行去相關(guān)濾波[18]，即對于某一次m(m≥6)，分別對奇數(shù)和偶數(shù)階系數(shù)進(jìn)行4次多項(xiàng)式擬合，利用最小二乘法確定多項(xiàng)式各擬合系數(shù)如下：

(4)

式中:l為階數(shù)，m為次數(shù)，ai、bi(i=0,1,2,3,4)為待求解的多項(xiàng)式未知參數(shù)。去除相關(guān)誤差后的球諧系數(shù)為

(5)

重力場模型系數(shù)去除相關(guān)誤差后，可將式(1)表示如下：

(6)

陸地水儲(chǔ)量變化時(shí)間序列主要包含趨勢信號、季節(jié)性信號(周年、半年信號)和周期為161 d的S2分潮，用式(7)對其進(jìn)行擬合[10]。

(7)

3 反演結(jié)果與分析

3.1 長江流域水儲(chǔ)量年際和年內(nèi)變化

利用式(6)求得長江流域0.25°×0.25°陸地水儲(chǔ)量變化等效水高，再將區(qū)域內(nèi)每個(gè)格網(wǎng)點(diǎn)的等效水高按面積加權(quán)求得流域內(nèi)平均水儲(chǔ)量變化時(shí)間序列。根據(jù)式(7)，對圖1中的平均陸地水儲(chǔ)量變化序列進(jìn)行擬合，具體研究時(shí)段內(nèi)水儲(chǔ)量的變化趨勢與振幅情況如表1所示。反演結(jié)果顯示：兩個(gè)模型反演的水儲(chǔ)量變化均呈增加的趨勢，分別為0.58 cm/a(Tongji-Grace2018)和0.63 cm/a(CSR RL06)，周年振幅為4.23 cm(Tongji-Grace2018)和4.04 cm(CSR RL06)。另外，本文的趨勢計(jì)算結(jié)果與其他研究結(jié)果基本一致[6-9]，存在差異的原因是研究時(shí)間段、模型及數(shù)據(jù)預(yù)處理方式和具體的研究范圍不完全一致。在本文的反演結(jié)果中，Tongji-Grace2018估計(jì)的周年振幅比CSR RL06的估計(jì)周年振幅要大，趨勢項(xiàng)和半年振幅略小。此外，由圖1可知：2006年和2011年的陸地水儲(chǔ)量變化偏小，2016年的陸地水儲(chǔ)量增長幅度最大。由此可以初步判定長江流域在2006年和2011年以及2016年發(fā)生了水文事件。該發(fā)現(xiàn)與其他研究結(jié)果相符：2006年和2011年是長江流域極端少雨事件的發(fā)生年份，2016年是多雨極值出現(xiàn)的年份，說明陸地水儲(chǔ)量的變化與ENSO以及降雨量的增減密切相關(guān)[8,13,19]。

圖2給出了兩個(gè)模型反演的長江流域水儲(chǔ)量年內(nèi)分布情況，兩個(gè)模型計(jì)算的年內(nèi)分布狀況基本一致：研究時(shí)間段內(nèi)，陸地水儲(chǔ)量變化年內(nèi)分配不均，1月最小，5～6月開始由虧轉(zhuǎn)盈，7月達(dá)到峰值，11月開始虧損。其中，7～9月是水儲(chǔ)量盈余的高峰階段，期間要做好防洪工作，1～4月水儲(chǔ)量虧損嚴(yán)重，要注意防旱。

圖2 2003～2016年長江流域陸地水儲(chǔ)量變化年內(nèi)分布Fig.2 Monthly average water storage change in the Yangtze River Basin from 2003 to 2016

3.2 長江流域水儲(chǔ)量變化空間分布

由于流域面積較大，平均水儲(chǔ)量變化值不能反映流域內(nèi)水儲(chǔ)量變化的空間分布信息，因此給出了陸地水儲(chǔ)量變化月平均值的空間分布(見圖3)。兩個(gè)模型反演的水儲(chǔ)量變化空間分布基本一致。1月份流域中下游以及四川省西部邊緣處于虧損狀態(tài)，流域內(nèi)以江西、湖南、湖北三省的水儲(chǔ)量虧損最為嚴(yán)重；3月份流域水儲(chǔ)量虧損程度減輕，源區(qū)水儲(chǔ)量變化很?。?月份流域全區(qū)的水儲(chǔ)量開始由虧轉(zhuǎn)盈。7月份長江中下游水儲(chǔ)量增長可達(dá)6 cm，這是由于6月份的中下旬和7月份是長江中下游地區(qū)的梅雨季，東南季風(fēng)帶來的太平洋暖濕氣流使流域的中下游持續(xù)陰雨天氣，在降雨增加的同時(shí)減少了陸地水的蒸散，使得中下游的陸地水儲(chǔ)量大幅增加。進(jìn)入9月份，源區(qū)以及湖北省北部和四川省的西南部保持輕微的盈余，其他地區(qū)的水儲(chǔ)量變化量級很小，處于平衡狀態(tài)。隨著降水的減少，11月份長江中下游地區(qū)的水儲(chǔ)量虧損，源區(qū)和上游的水儲(chǔ)量處于輕度虧損狀態(tài)，但程度低于1～4月份。

圖3 2003～2016年長江流域月平均水儲(chǔ)量變化空間分布Fig.3 Spatial distribution of monthly average water storage in the Yangtze River Basin from 2003 to 2016

根據(jù)流域的氣候特征，按季節(jié)劃分(3～5月為春季，6～8月為夏季，9～11月為秋季，12至次年2月為冬季)，給出了季節(jié)平均水儲(chǔ)量變化的空間分布(見圖4)。兩種模型的反演結(jié)果吻合度較高，長江流域陸地水儲(chǔ)量變化的季節(jié)性特征明顯。整體來看，水儲(chǔ)量夏季盈余，春冬兩季虧損(冬季虧損程度高于春季)，秋季處于基本平衡狀態(tài)，這與長江流域亞熱帶季風(fēng)性氣候夏季高溫多雨、冬季溫和少雨的特點(diǎn)相一致。

圖4 2003～2016年長江流域季節(jié)平均水儲(chǔ)量變化空間分布Fig.4 Spatial distribution of seasonal average water storage in the Yangtze River Basin from 2003 to 2016

3.3 長江流域各省陸地水儲(chǔ)量變化趨勢及振幅

為更清楚地了解長江流域內(nèi)各省陸地水儲(chǔ)量的變化情況，為預(yù)防以及應(yīng)對洪澇災(zāi)害提供保障，本節(jié)利用Tongji-Grace2018和CSR RL06模型分別分析了各省水儲(chǔ)量變化趨勢及周期等特性。結(jié)果顯示，兩個(gè)模型對各省的反演結(jié)果吻合較好。整體而言，除湖北省北部、安徽省北部以及四川省西部邊緣外，各省的陸地水儲(chǔ)量均呈現(xiàn)出逐年增加的特點(diǎn)。其中，江西和湖南兩省的增速最大，約為1 cm/a。而湖北省北部與安徽省北部的水儲(chǔ)量呈減少的趨勢，這與華北平原逐年嚴(yán)重的深層地下水虧損有關(guān)[20]。此外，四川省西部以及青海省南部的水儲(chǔ)量也逐年減少，主要是由全球氣候變暖導(dǎo)致青藏高原冰川融化所致[7]。這與利用面積加權(quán)計(jì)算的各省份的陸地水儲(chǔ)量變化趨勢中安徽省、湖北省和四川省的陸地水儲(chǔ)量的變化趨勢值較小的結(jié)果基本一致(見表2和圖5)。

圖5 2003～2016年長江流域陸地水儲(chǔ)量變化趨勢及周年、半年振幅Fig.5 Liner trend and annual,semiannual amplitude of TWSC in the Yangtze River Basin from 2003 to 2016

表2 長江流域各省陸地水儲(chǔ)量變化趨勢及振幅Tab.2 Liner trends and amplitudes of TWSC of the provinces in the Yangtze River Basin

流域內(nèi)陸地水儲(chǔ)量變化周年、半年振幅的分布有較強(qiáng)的空間異質(zhì)性。進(jìn)一步定量分析表明：湖南和江西兩省的周年振幅最大，約為6cm。除安徽省和湖北省之外，Tongji-Grace2018計(jì)算的其余省份的周年振幅均略大于CSR RL06；但對于半年振幅，結(jié)果相反；兩個(gè)模型計(jì)算的趨勢基本一致。

4 結(jié)果與討論

本研究利用Tongji-Grace2018月時(shí)變重力場模型反演了長江流域2003～2016年約14 a的陸地水儲(chǔ)量變化情況，清晰地展示了長江流域月平均以及季節(jié)平均的陸地水儲(chǔ)量變化，定量估計(jì)了流域內(nèi)各省的陸地水儲(chǔ)量變化趨勢與振幅，并將Tongji-Grace 2018的反演結(jié)果與CSR RL06的反演結(jié)果進(jìn)行了對比。研究結(jié)果對監(jiān)測長江流域的旱澇災(zāi)害，進(jìn)一步地分析蒸散、徑流以及地下水等的變化，完善長江流域水循環(huán)體系等具有重要的參考價(jià)值，具體結(jié)果如下。

(1) 長江流域水儲(chǔ)量變化呈上升的趨勢，而且季節(jié)特征明顯，兩個(gè)模型反演的水儲(chǔ)量變化趨勢基本一致，分別為0.58 cm/a(Tongji-Grace2018)和0.63 cm/a(CSR RL06)，周年振幅為4.23 cm(Tongji-Grace2018)和4.04 cm(CSR RL06)?？梢?，在長江流域Tongji-Grace2018反演的周年信號更明顯，而CSR RL06模型估計(jì)的趨勢項(xiàng)略大。反演結(jié)果的差異與模型解算策略有關(guān)。

(2) 月平均水儲(chǔ)量給出了不同月份陸地水儲(chǔ)量的變化情況：6～10月水儲(chǔ)量盈余，11月至次年5月水儲(chǔ)量虧損。其中，水儲(chǔ)量變化最小值出現(xiàn)在1月份，5～6月水儲(chǔ)量由虧轉(zhuǎn)盈，峰值出現(xiàn)在7月份，與長江流域的梅雨季節(jié)吻合。

(3) 長江流域水儲(chǔ)量變化空間異質(zhì)性較大，中下游地區(qū)水儲(chǔ)量的變化更為明顯，上游地區(qū)水儲(chǔ)量變化相對較小。兩個(gè)模型反演的流域季節(jié)平均水儲(chǔ)量的變化空間分布表現(xiàn)為：夏季盈余，春、冬季虧損，秋季變化平衡。

(4) 除了與華北平原、青藏高原接壤的流域邊緣地帶以外，流域內(nèi)各省的陸地水儲(chǔ)量變化均呈增加的趨勢。其中，湖南、江西兩省的陸地水儲(chǔ)量增速較大，約為1 cm/a。湖南省和江西省的周年振幅最大，約為6 cm。

GRACE時(shí)變重力場模型能夠揭示厘米級等效水高變化，是很好的質(zhì)量變化反演手段。但反演結(jié)果受誤差影響較重，誤差主要來自于GRACE衛(wèi)星的觀測誤差及軌道設(shè)計(jì)等。因此，針對GRACE模型進(jìn)行濾波處理，抑制模型噪聲并獲得更多水文信號的問題，一直是研究的熱點(diǎn)。雖然本研究選取了信噪比較高的Tongji-Grace2018模型進(jìn)行反演，同時(shí)進(jìn)行扇形濾波和去相關(guān)處理，反演結(jié)果中仍不可避免地存在噪聲。相信隨著GRACE Follow-On計(jì)劃的實(shí)施，陸地水儲(chǔ)量變化的空間分辨率以及監(jiān)測精度會(huì)有進(jìn)一步的提高。